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-- -. Abstract 

^ The fast progress in improving the sensitivity of the gravitational-wave detectors, we all 

^O have witnessed in the recent years, has propelled the scientific community to the point at which 

^^ quantum behavior of such immense measurement devices as kilometer-long interferometers 

1^ starts to matter. The time when their sensitivity will be mainly limited by the quantum noise 

. of light is around the corner, and finding ways to reduce it will become a necessity. Therefore, 

C^ the primary goal we pursued in this review was to familiarize a broad spectrum of readers 

^^ with the theory of quantum measurements in the very form it finds application in the area 

Cn of gravitational-wave detection. We focus on how quantum noise arises in gravitational-wave 

" interferometers and what limitations it imposes on the achievable sensitivity. We start from 

^ the very basic concepts and gradually advance to the general linear quantum measurement 

theory and its application to the calculation of quantum noise in the contemporary and planned 
interferometric detectors of gravitational radiation of the first and second generation. Special 
attention is paid to the concept of the Standard Quantum Limit and the methods of its 
surmounting. 



1 Introduction 

The more-than-ten- years-long history of the large-scale laser gravitation- wave (GW) detectors (the 
first one, TAMA j6j started to operate in 1999, and the most powerful pair, the two detectors of 
the LIGO project f^, in 2001, not to forget about the two European members of the international 
interferometric GW detectors network, also having a pretty long history, namely, the German- 
British interferometer GEO 600 [3 located near Hannover, Germany, and the joint European large- 
scale detector Virgo [7], operating near Pisa, Italy) can be considered both as a great success and 
a complete failure, depending on the point of view. On the one hand, virtually all technical 
requirements for these detectors have been met, and the planned sensitivity levels have been 
achieved. On the other hand, no GWs have been detected thus far. 

The possibility of this result had been envisaged by the community, and during the same last ten 
years, plans for the second-generation detectors were developed [144, 69, 11, 168, 1, 4 . Currently 
(2012), both LIGO detectors are shut down, and their upgrade to the Advanced LIGO, which 
should take about three years, is underway. The goal of this upgrade is to increase the detectors' 
sensitivity by about one order of magnitude p,,3.9j , and therefore the rate of the detectable events by 
three orders of magnitude, from some 'half per year' (by the optimistic astrophysical predictions) 
of the second generation detectors to, probably, hundreds per year. 

This goal will be achieved, mostly, by means of quantitative improvements (higher optical 
power, heavier mirrors, better seismic isolation, lower loss, both optical and mechanical) and 
evolutionary changes of the interferometer configurations, most notably, by introduction of the 
signal recycling mirror. As a result, the second-generation detectors will be quantum noise limited. 
At higher GW frequencies, the main sensitivity limitation will be due to phase fluctuations of 
light inside the interferometer (shot noise). At lower frequencies, the random force created by the 
amplitude fluctuations (radiation-pressure noise) will be the main or among the major contributors 
to the sum noise. 

It is important that these noise sources both have the same quantum origin, stemming from 
the fundamental quantum uncertainties of the electromagnetic field, and thus that they obey the 
Heisenberg uncertainty principle and can not be reduced simultaneously 03]. In particular, the 
shot noise can (and will, in the second generation detectors) be reduced by means of the optical 
power increase. However, as a result, the radiation-pressure noise will increase. In the 'naively' 
designed measurement schemes, built on the basis of a Michelson interferometer, kin to the first 
and the second generation GW detectors, but with sensitivity chiefly limited by quantum noise, 
the best strategy for reaching a maximal sensitivity at a given spectral frequency would be to 
make these noise source contributions (at this frequency) in the total noise budget equal. The 
corresponding sensitivity point is known as the Standard Quantum Limit (SQL) [51] [37] . 

This limitation is by no means an absolute one, and can be evaded using more sophisticated 
measurement schemes. Starting from the first pioneering works oriented on solid-state GW detec- 



tors [33 l 134 1 [145] ■ many methods of overcoming the SQL were proposed, including the ones suitable 
for practical implementation in laser-interferometer GW detectors. The primary goal of this review 
is to give a comprehensive introduction of these methods, as well as into the underlying theory of 
linear quantum measurements, such that it remains comprehensible to a broad audience. 

The paper is organized as follows. In Section ^ we give a classical (that is, non-quantum) 
treatment of the problem, with the goal to familiarize the reader with the main components of 
laser GW detectors. In Section [3] we provide the necessary basics of quantum optics. In Section [4] 
we demonstrate the main principles of linear quantum measurement theory, using simplified toy 
examples of the quantum optical position meters. In Section [5] we provide the full-scale quantum 
treatment of the standard Fabry-Perot-Michelson topology of the modern optical GW detectors. 
At last, in Section [6] we consider three methods of overcoming the SQL, which are viewed now 
as the most probable candidates for implementation in future laser GW detectors. Concluding 



remarks are presented in Section [7] Throughout the review we use the notations and conventions 
presented in Table [T] below. 

Table 1: Notations and conventions, used in this review, given in 
alphabetical order for both, greek (first) and latin (after greek) 
symbols. 
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SQL beating factor 
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annihilation and creation operators of photons with fre- 
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two-photon amplitude quadrature operator 

two-photon phase quadrature operator 

Symmetrised (cross) correlation of the field quadrature 
operators {i, j — c, s) 

light beam cross section area 
speed of light 

light quantization normalization constant 

Resonance denominator of the optical cavity transfer func- 
tion, defining its characteristic conjugate frequencies ("cav- 
ity poles") 

electric field strength 

classical complex amplitude of the light 

classical quadrature amplitudes of the light 

vector of classical quadrature amplitudes 

back-action force of the meter 

signal force 

dimensionless GW signal (a.k.a. metrics variation) 

homodyne vector 

Hamiltonian of a quantum system 

Planck's constant 

identity matrix 

optical power 

circulating optical power in a cavity 

circulating optical power per interferometer arm cavity 

normalized circulating power 
optical pumping wave number 
rigidity, including optical rigidity 
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Kimble's optomechanical coupling factor 
optomechanical coupling factor of the Sagnac speedmeter 

cavity length 

probe-body mass 

general linear meter readout observable 

matrix of counterclockwise rotation (pivoting) by angle a 

amplitude squeezing factor (e'') 

power squeezing factor in decibels 

power reflectivity of a mirror 

reflection matrix of the Fabry-Perot cavity 

noise power spectral density (double-sided) 

measurement noise power spectral density (double-sided) 

back-action noise power spectral density (double-sided) 

cross-correlation power spectral density (double-sided) 

vacuum quantum state power spectral density matrix 

squeezed quantum state power spectral density matrix 

[—0] squeezing matrix 

power transmissivity of a mirror 

transmissivity matrix of the Fabry-Perot cavity 

test-mass velocity 

optical energy 

Wigncr function of the quantum state ji/") 

test-mass position 

dimensionless oscillator (mode) displacement operator 

dimensionless oscillator (mode) momentum operator 



2 Interferometry for GW Detectors: Classical Theory 

2.1 Interferometer as a weak force probe 

In order to have a firm basis for understanding liow quantum noise influences the sensitivity of a 
GW detector it would be illuminating to give a brief description of the interferometers as weak 
force/tiny displacement meters. It is by no means our intention to give a comprehensive survey 
of this ample field that is certainly worthy of a good book, which there are in abundance, but 
rather to provide the reader with the wherewithal for grasping the very principles of the GW 
interferometers operation as well as of other similar ultrasensitive optomechanical gauges. The 
reader interested in a more detailed description of the interferometric techniques being used in the 
field of GW detectors might enjoy reading this book |T7] or the comprehensive Living Reviews on 
the subject by Freise and Strain [65] and by Pitkin et al. |125j . 

2.1.1 Light phase as indicator of a weak force 



50 = 47rg 




Figure 1: Scheme of a simple weak force measurement: an external signal force G pulls the mirror 
from its equilibrium position a: = 0, causing displacement 5x. The signal displacement is measured 
by monitoring the phase shift of the light beam, reflected from the mirror. 

Let us, for the time being, imagine that we are capable of measuring an electromagnetic (e.g., 
light) wave-phase shift 54> with respect to some coherent reference of the same frequency. Having 
such a hypothetical tool, what would be the right way to use it, if one had a task to measure 
some tiny classical force? The simplest device one immediately conjures up is the one drawn in 
Figure [l] It consists of a movable totally-reflective mirror with mass M and a coherent paraxial 
light beam, that impinges on the mirror and then gets reflected towards our hypothetical phase- 
sensitive device. The mirror acts as a probe for an external force G that one seeks to measure. 
The response of the mirror on the external force G depends upon the details of its dynamics. For 
definiteness, let the mirror be a harmonic oscillator with mechanical eigenfrequency f^m = 27rfm- 
Then the mechanical equation of motion gives a connection between the mirror displacement x 
and the external force G in the very familiar form of the harmonic oscillator equation of motion: 

Mi + M^lx = Git) , =^ x{t) = xo{t) + f dt' xUt - t')G{t'), (1) 

Jo 

where a;o(t) = x(0) cosr^^i + p(0)/(Mrim) sinri^i is the free motion of the mirror defined by its 



initial displacement a;(0) and momentum p{0) at t = and 

is the oscillator Green's function. It is easy to see that the reflected light beam carries in its 
phase the information about the displacement Sx{t) — x{t) — x{0) induced by the external force G. 
Indeed, there is a phase shift between the incident and reflected beams, that matches the additional 
distance the light must propagate to the new position of the mirror and back, i.e., 

2ujqSx 6x 

6(j) = = 47r— , (3) 

c Ao 

with Wo = 2-Kc/Xo the incident light frequency, c the speed of light and Aq the light wavelength. 
Here we implicitly assume mirror displacement to be much smaller than the light wavelength. 

Apparently, the information about the signal force G{t) can be obtained from the measured 
phase shift by post-processing of the measurement data record 5(j){t) ex 5x(t) by substituting it into 
Eq. (fTl) instead of x. Thus, the estimate of the signal force G reads: 
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(4) 



This kind of post-processing pursues an evident goal of getting rid of any information about the 

eigenmotion of the test object while keeping only the signal-induced part of the total motion. 

The above time-domain expression can be further simplified by transforming it into a Fourier 

domain, since it does not depend anymore on the initial values of the mirror displacement x{Qi) 

and momentum p(0): 

Mr 
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denotes a Fourier transform of an arbitrary time-domain function A{t). If the expected signal 
spectrum occupies a frequency range that is much higher than the mirror-oscillation frequency ^m 
as is the case for ground based interferometric GW detectors, the oscillator behaves as a free mass 
and the term proportional to fJ,^ in the equation of motion can be omitted yielding: 

G'--{n) = -^p^Ha. (7) 

2a;o 

2.1.2 Michelson interferometer 

Above, we assumed a direct light phase measurement with a hypothetical device in order to detect 
a weak external force, possibly created by a GW. However, in reality, direct phase measurement 
are not so easy to realize at optical frequencies. At the same time, physicists know well how to 
measure light intensity (amplitude) with very high precision using different kinds of photodetectors 
ranging from ancient -yet-die-hard reliable photographic plates to superconductive photodetectors 
capable of registering individual photons (71j . How can one transform the signal, residing in 
the outgoing light phase, into amplitude or intensity variation? This question is rhetorical for 
physicists, for interference of light as well as the multitude of interferometers of various design and 
purpose have become common knowledge since a couple of centuries ago. Indeed, the amplitude 
of the superposition of two coherent waves depends on the relative phase of these two waves, thus 
transforming phase variation into the variation of the light amplitude. 
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Figure 2: Scheme of a Michelson interferometer. When the end mirrors of the interferometer arms 
Mn^e are at rest the length of the arms L is such that the hght from the laser gets reflected back 
entirely (bright port), while at the dark port the reflected waves suffer destructive interference 
keeping it really dark. If, due to some reason, e.g., because of GWs, the lengths of the arms 
changed in such a way that their difference was 5L, the photodetector at the dark port should 
measure light intensity /dark(i5-^) = "2^(1 ~ cos47r^). 



For the detection of GWs, the most popular design is the Michelson interferometer PUIITTII^ , 
which schematic view is presented in Figure [2J Let us briefly discuss how it works. Here, the light 
wave from a laser source gets split by a semi-transparent mirror, called a beamsplitter, into two 
waves with equal amplitudes, travelling towards two highly-reflective mirrors A/„ jMto get reflected 
off them, and then recombine at the beamsplitter. The readout is performed by a photodetector, 
placed in the signal port. The interferometer is usually tuned in such a way as to operate at a dark 
fringe, which means that by default the lengths of the arms are taken so that the optical paths for 
light, propagating back and forth in both arms, are equal to each other, and when they recombine 
at the signal port, they interfere destructively, leaving the photodetector unilluminated. On the 
opposite, the two waves coming back towards the laser, interfere constructively. The situation 
changes if the end mirrors get displaced by some external force in a differential manner, i.e., such 
that the difference of the arms lengths is non-zero: SL = L^ — Ln ^ 0. Let a laser send to the 
interferometer a monochromatic wave that, at the beamsplitter, can be written as 

Euscrit) = Eocos{ujot)- 
Hence, the waves reflected off the interferometer arms at the beamsplitter (before interacting with 



^Here, we adopt the system of labeling parts of the interferometer by the cardinal directions, they are located 
with respect to the interferometer central station, e.g., M„ and A/e in Figure p] stand for 'northern' and 'eastern' 
end mirrors, respectively. 
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And the intensity of the outgoing hght in both ports can be found using a relation I oc E^ with 
overline meaning time-average over many osciUation periods: 

Idark(<5i/Ao) = y ( 1 ~ cos47r— j , and Ibright('5i/Ao) = y ( 1 + cos47r— j . (8) 

Apparently, for small differential displacements 6L ^ Aq, the Michelson interferometer tuned 
to operate at the dark fringe has a sensitivity to ^ {5L/Xo)'^ that yields extremely weak light 
power on the photodetector and therefore very high levels of dark current noise. In practice, the 
interferometer, in the majority of cases, is slightly detuned from the dark fringe condition that can 
be viewed as an introduction of some constant small bias SLq between the arms lengths. By this 
simple trick experimentalists get linear response to the signal nonstationary displacement Sx{t): 

Idark(oa;/Ao) = — l-cos47r 



2 V Ao 



Comparison of this formula with Eq. ([3]) should immediately conjure up the striking similarity 
between the response of the Michelson interferometer and the single moving mirror. The nonsta- 
tionary phase difference of light beams in two interferometer arms S4>{t) = ATT6x{t) / Xq is absolutely 
the same as in the case of a single moving mirror (cf. Eq. ([3])). It is no coincidence, though, but 
a manifestation of the internal symmetry that all Michelson-type interferometers possess with re- 
spect to coupling between mechanical displacements of their arm mirrors and the optical modes of 



the outgoing fields. In the next section 2.1.3 we show how this symmetry displays itself in GW 
interferometers . 

2.1.3 Gravitational waves' interaction with interferometer 

Let us see how a Michelson interferometer interacts with the GW. For this purpose we need to 
understand, on a very basic level, what a GW is. Following the poetic, yet precise, definition by Kip 
Thorne, 'gravitational waves are ripples in the curvature of spacetime that are emitted by violent 
astrophysical events, and that propagate out from their source with the speed of light' [I8 llll2j . 
A weak GW far away from its birthplace can be most easily understood from analyzing its action 
on the probe bodies motion in some region of spacetime. Usually, the deformation of a circular 



■^Here and below we keep to a definition of the reflectivity coefliicient of the mirrors that implies that the reflected 
wave acquires a phase shift equal to it with respect to the incident wave if the latter impinged the reflective surface 
from the less optically dense medium (air or vacuum). In the opposite case, when the incident wave encounters 
reflective surface from inside the mirror, i.e., goes from the optically more dense medium (glass), it is assumed to 
acquire no phase shift upon reflection. 



ring of free test particles is considered (see Chapter 26: Section 26.3.2 of [THj for more rigorous 
treatment) when a GW impinges it along the z-direction, perpendicular to the plane where the 
test particles are located. Each particle, having plane coordinates (x, y) with respect to the center 
of the ring, undergoes displacement Sr = (Sx, Sy) from its position at rest, induced by GWs: 



Sx ■ 
5x 



/l_|_X , 



Sy 



-h+y, 



1. 



-/ixy, Sy = -h-^x. 



(10) 
(11) 



Here, /i+ = hj^-lt — z/c) and hx = hx{t — z/c) stand for two independent polarizations of a GW that 
creates an acceleration field resulting in the above deformations. The above expressions comprise 
a solution to the equation of motion for free particles in the tidal acceleration field created by a 
GW: 

Sr = - (h+x+'hxy)e^x + {-'h+y+'hxx)ey , 

with e^ — {1, 0}^ and By = {0, 1}^ the unit vectors pointing in the x and y direction, respectively. 



a) /i+-polarized GW 
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b) /ix -polarized GW 
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Figure 3: Action of the GW on a Michelson interferometer: (a) /i+-polarized GW periodically 
stretch and squeeze the interferometer arms in the x- and j/-directions, (b) /ix -polarized GW though 
have no impact on the interferometer, yet produce stretching and squeezing of the imaginary test 
particle ring, but along the directions, rotated by 45° with respect to the x and y directions of the 
frame. The lower pictures feature field lines of the corresponding tidal acceleration fields ex /i+,x- 

For our Michelson interferometer, one can consider the end mirrors to be those test particles 
that lie on a circular ring with beamsplitter located in its center. One can choose arms directions 
to coincide with the frame x and y axes, then the mirrors will have coordinates (0, L„) and (Lg, 0), 
correspondingly. For this case, the action of the GW field on the mirrors is featured in Figure [3] 
It is evident from this picture and from the above formulas that an /i x -polarized component of the 
GW does not change the relative lengths of the Michelson interferometer arms and thus does not 
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contribute to its output signal; at the same time, /i^_-polarized GWs act on the end masses of the 
interferometer as a pair of tidal forces of the same value but opposite in direction: 

Gn = —-Mnhj^Ln , Gg — -Mg^h-\.Lf. . 

Assuming Ge — — G„ = G, M„ = Me = M, and L^ ~ i„ = L, one can write down the 
equations of motion for the interferometer end mirrors that are now considered free (f2m <C J^gw) 
as: 

Mi^G, My = -G , 

and for the differential displacement of the mirrors SL = Lg — Ln = x — y, which, we have shown 
above, the Michelson interferometer is sensitive to, one gets the following equation of motion: 

M5L^2G(t)^Mh+{t)L (12) 

that is absolutely analogous to Eq. (fTl) for a single free mirror with mass M. Therefore, we have 
proven that a Michelson interferometer has the same dynamical behavior with respect to the tidal 
force G{t) — Mh^{t)L/2 created by GWs, as the single movable mirror with mass M to some 
external generic force G{t). 

The foregoing conclusion can be understood in the following way: for GWs are inherently 
quadruple and, when the detector's plane is orthogonal to the wave propagation direction, can only 
excite a differential mechanical motion of its mirrors, one can reduce a complicated dynamics of 
the interferometer probe masses to the dynamics of a single effective particle that is the differential 
motion of the mirrors in the arms. This useful observation appears to be invaluably helpful for 
calculation of the real complicated interferometer responses to GWs and also for estimation of its 
optical quantum noise, that comprises the rest of this review. 

2.2 Prom incident wave to outgoing light: light transformation in the 
GW interferometers 

To proceed with the analysis of quantum noise in GW interferometers we first need to familiarize 
ourselves with how a light field is transformed by an interferometer and how the ability of its 
mirrors to move modifies the outgoing field. In the following paragraphs, we endeavor to give a 
step-by-step introduction to the mathematical description of light in the interferometer and the 
interaction with its movable mirrors. 

2.2.1 Light propagation 

We first consider how the light wave is described and how its characteristics transform, when 
it propagates from one point of free space to another. Yet the real light beams in the large 
scale interferometers have a rather complicated inhomogeneous transverse spatial structure, the 
approximation of a plane monochromatic wave should suffice for our purposes, since it comprises 
all the necessary physics and leads to right results. Inquisitive readers could find abundant material 
on the field structure of light in real optical resonators in particular, in the introductory book |171| 
and in the Living Review by Vinet |155) . 

So, consider a plane monochromatic linearly polarized light wave propagating in vacuo in the 
positive direction of the a;-axis. This field can be fully characterized by the strength of its electric 
component E{t — x/c) that should be a sinusoidal function of its argument C ~ t^ ^/c and can be 
written in three equivalent ways: 

E{() = £o cos [woC - (f>o] = £c cos woC + £s sinwoC = tjf ■, (13) 
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where £o and (f>Q are called amplitude and phase, E^ and Eg take names of cosine and sine quadra- 
ture amplitudes, and complex number £ — \£\e'-^^^^ is known as the complex amplitude of the 
electromagnetic wave. Here, we see that our wave needs two real or one complex parameter to be 
fully characterized in the given location x at a given time t. The 'amplitude-phase' description is 
traditional for oscillations but is not very convenient since all the transformations are nonlinear 
in phase. Therefore, in optics, either quadrature amplitudes or complex amplitude description 
is applied to the analysis of wave propagation. All three descriptions are related by means of 
straightforward transformations: 



£o = v/^F+f| 



/2\£\ 



£c = ^^ = \/2Rc [£] = £o cos ( 



V2 



V2^ 



tancio — £s/£c = 8.Tg£, (f>ohi[0, 2it\ 
£s = ^-^ = \/2Im [£] = £o sin(/)o , 

£* = ^■=-'^- 



(14) 



V2 



— £a_p — i4>a 



The aforesaid means that for complete understanding of how the light field transforms in the 
optical device, knowing the rules of transformation for only two characteristic real numbers - real 
and imaginary parts of the complex amplitude suffice. Note also that the electric field of a plane 
wave is, in essence, a function of a single argument C^ — t — x/c (for a forward propagating wave) 
and thus can be, without loss of generality, substituted by a time dependence of electric field in 
some fixed point, say with x = Q, thus yielding £'(C) ^ E{t). We will keep to this convention 
throughout our review. 

Now let us elaborate the way to establish a link between the wave electric field strength values 
taken in two spatially separated points, xi = Q and X2 = L. Obviously, if nothing obscures 
light propagation between these two points, the value of the electric field in the second point at 
time t is just the same as the one in the first point, but at earlier time, i.e., at t' = t — L/c: 
£'(^) (i) = iijC^) {t — L/c) . This allows us to introduce a transformation that propagates EM- wave 
from one spatial point to another. For complex amplitude £, the transformation is very simple: 



g{L) ^ gic^oL/c£(0) 



(15) 



Basically, this transformation is just a counterclockwise rotation of a wave complex amplitude 
vector on a complex plane by an angle (J)l = [^^^] j o • This fact becomes even more evident if 
we look at the transformation for a 2-dimensional vector of quadrature amplitudes £ = {£c,£s\ , 
that are: 



£(L)^ 



COS<pL 

sin (t>L 



- sm (t>L 

COS(f)L 



where 
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cost! 
sin 9 



— sin 6 
cos 6* 
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(16) 



(17) 



stands for a standard counterclockwise rotation (pivoting) matrix on a 2D plane. In the special 
case when the propagation distance is much smaller than the light wavelength L <C A, the above 
two expressions can be expanded into Taylor's series in (f>L = 2TrL/X <C 1 up to the first order: 



S(^«^) = (1 + z,/.l)^(°) 



and 

£:(i«A) 
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(18) 



(I + ,5P[0i])f("), (19) 
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where I stands for an identity matrix and i5P[(?!)i] is an infinitesimal increment matrix that generate 
the difference between the field quadrature amphtudes vector £ after and before the propagation, 
respectively. 

It is worthwhile to note that the quadrature amplitudes representation is used more frequently 
in literature devoted to quantum noise calculation in GW interferometers than the complex am- 
plitudes formalism and there is a historical reason for this. Notwithstanding the fact that these 
two descriptions are absolutely equivalent, the quadrature amplitudes representation was chosen 
by Caves and Schumaker as a basis for their two-photon formalism for the description of quantum 
fluctuations of light [HlllS] that became from then on the workhorse of quantum noise calculation. 



More details about this extremely useful technique are given in the sections 3.1 and 3.2 of this 
review. Unless otherwise specified, we predominantly keep ourselves to this formalism and give all 
results in terms of it. 

2.2.2 Modulation of light 

Above, we have seen that a GW signal displays itself in the modulation of the phase of light, passing 
through the interferometer. Therefore, it is illuminating to see how the modulation of the light 
phase and/or amplitude manifests itself in a transformation of the field complex amplitude and 
quadrature amplitudes. Throughout this section we assume our carrier field is a monochromatic 
light wave with frequency wq, amplitude Sq and initial phase (f>o — 0: 

^car(i) = So COS woi = Re [foe"'"°*] . 

Amplitude modulation. The modulation of light amplitude is straightforward to analyze. Let 
us do it for pedagogical sake: imagine one managed to modulate the carrier field amplitude slow 
enough compared to the carrier oscillation period, i.e., fi <C wo, then: 

EAuit) == £o(l + emC08{flt + (f)m)) COS UJQt, 

where em ^ 1 and 4>m are some constants called modulation depth and relative phase, respectively. 
The complex amplitude of the modulated wave equals to 

c 

£AM{t) = ^ (1 + e™ cos(rJi + 0,„)) , 

and the carrier quadrature amplitudes are, apparently, transformed as follows: 
ScAuit) ^£oil + erncos{flt + (l)m)) and Ss.AMit) ^0. 

The fact that the amplitude modulation shows up only in the quadrature that is in phase with the 
carrier field sets forth why this quadrature is usually named amplitude quadrature in the literature. 
In our review, we shall also keep to this terminology and refer to cosine quadrature as amplitude 
one. 

Illuminating also is the calculation of the modulated light spectrum, that in our simple case of 
single frequency modulation is straightforward: 



^AM(t) = Re 



£^g-iLUnt _|_ ^oemg-i0™g-j(LJo+n)t _|_ ^^l£!Zlgi</.™g-i(cun-f2)t 



Apparently, the spectrum is discrete and comprises three components, i.e., the harmonic at carrier 
frequency ujq with amplitude Ai^^ = £q and two satellites at frequencies ujq ± f2, also referred 
to as modulation sidebands, with (complex) amplitudes A^g±Q = emi5oe^*'^'"/2. The graphical 
interpretation of the above considerations is given in the left panel of Figured Here, carrier fields 
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as well as sidebands are represented by rotating vectors on a complex plane. The carrier field vector 
has length Eq and rotates clockwise with the rate ujq, while sideband components participate in 
two rotations at a time. The sum of these three vectors yields a complex vector, whose length 
oscillates with time, and its projection on the real axis represents the amplitude-modulated light 
field. 

The above can be generalized to an arbitrary periodic modulation function A(t) — 
X]fe°=i ^fc cos{kil. + (j)k), with -Eam(^) ~ ^o(l + £m^(^)) coswoi. Then the spectrum of the modulated 
light consists again of a carrier harmonic at wq and an infinite discrete set of sideband harmonics 
at frequencies ujq ± kfl [k — 1, cx)): 

^ oo 

EAuit) = focoswot+ "^ y^^fc {cos[(i:jo - kn)t - (pk] + cos[(wo + kil.)t + (j)k]} ■ (20) 

fe=i 



Further generalization to an arbitrary (real) non-periodic modulation function A{t) 

duj 
oo 27r ^ 



J^^^A{n)e-'^* is apparent: 



foe-'""* + e^Soe'"^"' / —A{n)e-'''^* 



EAuit) - Re 

So cos^oi + ^ / ^ {A{co - ^o) + A{lj + LJo)} e"^"* . (21) 

2 J-oo 27r 

From the above expression, one readily sees the general structure of the modulated light spectrum, 
i.e., the central carrier peaks at frequencies iwo and the modulation sidebands around it, whose 
shape retraces the modulation function spectrum A{uj) shifted by the carrier frequency ±ujq. 

Phase modulation. The general feature of the modulated signal that we pursued to demonstrate 
by this simple example is the creation of the modulation sidebands in the spectrum of the modulated 
light. Let us now see how it goes with a phase modulation that is more related to the topic of the 
current review. The simplest single-frequency phase modulation is given by the expression: 

Epuit) — So cos{ujot + Sjn cos{Qt + 0m)) , 



where fi <C Wq, and the phase deviation (5^ is assumed to be much smaller than 1. Using Eqs. (14 1, 
one can write the complex amplitude of the phase-modulated light as: 

p f^^^ — ^0 i5„cos(nt+0„) 
opM[t) — -y=e ^ ' , 

and quadrature amplitudes as: 

ScPuit) = So cos [Sm cos(51i + 4>m)] and Ss^puit) = So sin [6m cos{Qt + 0^)] . 

Note that in the weak modulation limit {dm <C 1), the above equations can be approximated as: 

ScPuit) — So and Ss^puit) ^ S^So cos{Qt + 0„i) . 

This testifies that for a weak modulation only the sine quadrature, which is tt/2 out-of-phase with 
respect to the carrier field, contains the modulation signal. That is why this sine quadrature is 
usually referred to as phase quadrature. It is also what we will call this quadrature throughout the 
rest of this review. 
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a) Amplitude modulation phasor diagram 



Im / 
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b) Phase modulation phasor diagram 
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Figure 4: Phasor diagrams for amplitude {Left panel) and phase {Right panel) modulated light. 
Carrier field is given by a brown vector rotating clockwise with the rate wq around the origin of the 
complex plane frame. Sideband fields are depicted as blue vectors. The lower (wq — Vi) sideband 
vector origin rotates with the tip of the carrier vector, while its own tip also rotates with respect 
to its origin counterclockwise with the rate f2. The upper {ljq + J7) sideband vector origin rotates 
with the tip of the upper sideband vector, while its own tip also rotates with respect to its origin 
counterclockwise with the rate fJ. Modulated oscillation is a sum of these three vectors an is given 
by the red vector. In the case of amplitude modulation (AM), the modulated oscillation vector 
is always in phase with the carrier field while its length oscillates with the modulation frequency 
i7. The time dependence of its projection onto the real axis that gives the AM-light electric field 
strength is drawn to the right of the corresponding phasor diagram. In the case of phase modulation 
(PM), sideband fields have a 7r/2 constant phase shift with respect to the carrier field (note factor 
i in front of the corresponding terms in Eq. (22 1; therefore its sum is always orthogonal to the 



carrier field vector, and the resulting modulated oscillation vector (red arrow) has approximately 
the same length as the carrier field vector but outruns or lags behind the latter periodically with 
the modulation frequency J7. The resulting oscillation of the PM light electric field strength is 
drawn to the right of the PM phasor diagram and is the projection of the PM oscillation vector 
on the real axis of the complex plane. 



In order to get the spectrum of the phase-modulated light it is necessary to refer to the theory 
of Bessel functions that provides us with the following useful relation (known as the Jacobi-Anger 
expansion) : 



IQm COS 



(f2t+0,„) 



Y. i''MSrn)e'' 



k{nt + <prri) 



/i.' — — oo 



where Jk{Sm) stands for the fc-th Bessel function of the first kind. This looks a bit intimidating, 
yet for (5m <C 1 these expressions simplify dramatically, since near zero Bessel functions can be 
approximated as: 



J0{Sm)^l-^-f+O{St,), MS„,)^^f+0{5i), MSrn) = l(^-J 



+ o{5^;:'){k^2) 



Thus, for sufficiently small (5„i, we can limit ourselves only to the terms of order 5^-^ and (5^ 
which yields: 



£;pM(i)^Re 



.Om^O 



fpe*""* + i ° |'gi[(wo+n)t+0™] _|_ ^i[(u:o-n)t~4>m]\ 



(22) 
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and we again face the situation in which modulation creates a pair of sidebands around the carrier 
frequency. The difference from the amphtude modulation case is in the way these sidebands behave 
on the complex plane. The corresponding phasor diagram for phase modulated light is drawn in 
Figure l4] In the case of PM, sideband fields have 7r/2 constant phase shift with respect to the 
carrier field (note factor i in front of the corresponding terms in Eq. (22)); therefore its sum is 



always orthogonal to the carrier field vector, and the resulting modulated oscillation vector has 
approximately the same length as the carrier field vector but outruns or lags behind the latter 
periodically with the modulation frequency O. The resulting oscillation of the PM light electric 
field strength is drawn to the right of the PM phasor diagram and is the projection of the PM 
oscillation vector on the real axis of the complex plane. 

Let us now generalize the obtained results to an arbitrary modulation function $(i): 

Epuit) = £0 cos(wot + (5,„$(i)) . 

In the most general case of arbitrary modulation index 5m, the corresponding formulas are very 
cumbersome and do not give much insight. Therefore, we again consider a simplified situation of 
sufficiently small (5^ ^ 1- Then one can approximate the phase- modulated oscillation as follows: 



EpM{t)=Re 



Eqc 



-iwotgi'5m*(t) 



Re[£oe-'"''*{l + i^™$(i)}] 



When <3?(i) is a periodic function: <I>(i) = J^kLi ^k cos kil + (p^, and in weak modulation limit 
(5to ^ 1, the spectrum of the PM light is apparent from the following expression: 



Epuit) 



Sm£o 



£0 cosujot ^ — 2_. ^k {sin [(wo — kQ)t — cjjk] + sin [{loq + kQ)t + 0^]} 



fc=i 



Re 



£oe' 



-iujQt 



iSm^O 



00 



-il{uio~kn)t~it>k] I ^-i[{(^o+kn)t+4>k] 



k=l 



(23) 



while for the real non-periodic modulation function $(i) = /-^i^*^(^)s *^^* ^^e spectrum can be 
obtained from the following relation: 



EpM (t) 



Re 



S^e-'^ot ^ i^„^£pe— 0* 






^{n)e 



-iflt 



£0 cos Wo i 



5m£a 



duj 

■ 2^ 



{i^{uj - uq) - i^{uj + Wo)} e 



(24) 



And again we get the same general structure of the spectrum with carrier peaks at zLujq and 
shifted modulation spectra i$(a;±wo) as sidebands around the carrier peaks. The difference with 
the amplitude modulation is an additional ±7r/2 phase shifts added to the sidebands. 

2.2.3 Laser noise 

Thus far we have assumed the carrier field to be perfectly monochromatic having a single spectral 
component at carrier frequency wo fully characterized by a pair of classical quadrature amplitudes 
represented by a 2-vector £. In reality, this picture is no good at all; indeed, a real laser emits not 
a monochromatic light but rather some spectral line of finite width with its central frequency and 
intensity fiuctuating. These fluctuations arc usually divided into two categories: (i) quantum noise 
that is associated with the spontaneous emission of photons in the gain medium, and (ii) technical 
noise arising, e.g., from excess noise of the pump source, from vibrations of the laser resonator, or 
from temperature fluctuations and so on. It is beyond the goals of this review to discuss the details 
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of the laser noise origin and methods of its suppression, since there is an abundance of literature 
on the subject that a curious reader might find interesting, e.g., the following works [1211 11221 11231 

nm[mm\. 

For our purposes, the very existence of the laser noise is important as it makes us to reconsider 
the way we represent the carrier field. Apparently, the proper account for laser noise prescribes 
us to add a random time-dependent modulation of the amplitude (for intensity fluctuations) and 
phase (for phase and frequency fluctuations) of the carrier field (13 1: 



E{t) = {£q + e„(i)) cos woi + 0o + (f>n{t) 



where we placed hats above the noise terms on purpose, to emphasize that quantum noise is a part 
of laser noise and its quantum nature has to be taken into account, and that the major part of this 
review will be devoted to the consequences these hats lead to. However, for now, let us consider 
hats as some nice decoration. 

Apparently, the corrections to the amplitude and phase of the carrier light due to the laser noise 
are small enough to enable us to use the weak modulation approximation as prescribed above. In 
this case one can introduce a more handy amplitude and phase quadrature description for the laser 
noise contribution in the following manner: 



E{t) — [£c + ec{t)) cosLOot + [£s + es(i)) sinwo^ , 



(25) 



where e^s are related to e„ and 0„ in the same manner as prescribed by Eqs. (14 1. It is convenient 
to represent a noisy laser field in the Fourier domain: 



,{t) 






,(fi)e 



-iQ.t 



Worth noting is the fact that Cc^si^) is a spectral representation of a real quantity and thus 
satisfies an evident equality e\ s(— f2) ~ ec,s{^) (by f we denote the Hcrmitian conjugate that for 
classical functions corresponds to taking the complex conjugate of this function). What happens if 
we want to know the light field of our laser with noise at some distance L from our initial reference 
point X = 0? For the carrier field component at wq, nothing changes and the corresponding 



transform is given by Eq. (16), yet for the noise component 



5Enoisc{t) = e.c{t) COS ujQt + es{t) sinujot 

there is a slight modification. Since the field continuity relation holds for the noise field to the 
same extent as for the carrier field: 



SEl!^J^^{t)= SEZs.it -L/c) 



the following modification applies: 



i(L] 



m 



e^^ 



m 



JlVlLJc 



UJqL 
COS 

. ^oL 

sm 

c 



— sm 



LO^L 



Wo 



I 






^iQL/cv. 



LOqL 



j("H^)- (26) 



Therefore, for sideband field components the propagation rule shall be modified by adding a 
frequency-dependent phase factor e*^^^/'^ that describes an extra phase shift acquired by a sideband 
field relative to the carrier field because of the frequency difference fl = uj — luq. 
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2.2.4 Light reflection from optical elements 



So, we are one step closer to understanding how to calculate the quantum noise of the light 
coming out of the GW interferometer. It is necessary to understand what happens with light 
when it is reflected from such optical elements as mirrors and beamsplitters. Let us first consider 
these elements of the interferometer fixed at their positions. The impact of mirror motion will be 
considered in the next subsection |2.2.5[ One can also refer to Section 2 of the Living Review by 
Freise and Strain f65' for a more detailed treatment of this topic. 
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Figure 5: Scheme of light reflection off the coated mirror. 

Mirrors of the modern interferometers are rather complicated optical systems usually consisting 
of a dielectric slab with both surfaces covered with multilayer dielectric coatings. These coatings 
are thoroughly constructed in such a way as to make one surface of the mirror highly reflective, 
while the other one is anti-reflective. We will not touch on the aspects of coating technology in 
this review and would like to refer the interested reader to an abundant literature on this topic, 
e.g., to the following book [7S] and reviews and articles [Hi [73 [Ml \IE [HI [Ml [IS Ell ^■ 
For our purposes, assuming the reflective surface of the mirror is flat and lossless should suffice. 
Thus, we represent a mirror by a reflective plane with (generally speaking, complex) coefficients of 
reflection r and r' and transmission t and t' as drawn in Figure^ Let us now see how the ingoing 
and outgoing light beams couple on the mirrors in the interferometer. 



Mirrors: From the general point of view, the mirror is a linear system with 2 input and 2 output 
ports. The way how it transforms input signals into output ones is featured by a 2 x 2 matrix that 
is known as the transfer matrix of the mirror M: 



(27) 



Since we assume no absorption in the mirror, reflection and transmission coefficients should satisfy 
Stokes' relations [IlTl[2n] (see also Section 12.12 of [101^ ): 



'E1''\tj 
ET\i) 


= M- 


'eW) 

Ef{t)_ 


= 


r t 
t' r' 




'ET{t) 

E^\t)_ 



t 



\t'\ 



W^i 



r*t' + r't* = , 



= 0, 



(28) 



that is simply a consequence of the conservation of energy. This conservation of energy yields that 
the optical transfer matrix M must be unitary: M^ = M~^. Stokes' relations leave some freedom 
in defining complex reflectivity and transmissivity coefficients. Two of the most popular variants 
are given by the following matrices: 



Jisym 



R 

is/T 



WT 

Vr 



and 



Ureal — 



R 
f 



(29) 
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where we rewrote transfer matrices in terms of real power reflectivity and transmissivity coefficients 
R = |rp and T — |ip that will find extensive use throughout the rest of this review. 

The transformation rule, or putting it another way, coupling relations for the quadrature am- 
plitudes can easily be obtained from Eq. (27 1. Now, we have two input and two output fields. 



Therefore, one has to deal with 4-dimensional vectors comprising of quadrature amplitudes of 
both input and output fields, and the transformation matrix become 4 x 4-dimensional, which can 
be expressed in terms of the outer product of a 2 x 2 matrix Mjcai by a 2 x 2 identity matrix I: 
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The same rules apply to the sidebands of each carrier field: 
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(30) 
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(31) 
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Figure 6: Scheme of a beamsplitter. 
In future, for the sake of brevity, we reduce the notation for matrices like M^eai I to simply 

^rcal- 



Beam splitters: Another optical element ubiquitous in the interferometers is a beamsplitter 
(see Figure pi). In fact, it is the very same mirror considered above, but the angle of input light 
beams incidence is different from (if measured from the normal to the mirror surface). The 
corresponding scheme is given in Figure p] In most cases, symmetric 50%/50% beamsplitter are 
used, which imply R = T = 1/2 and the coupling matrix M50/50 then reads: 



•150/50 



-l/\/2 l/\/2 

-l/\/2 lly/2 

l/\/2 l/%/2 

l|^/2 1/^/2 



(32) 
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Losses in optical elements: Above, we have made one assumption that is a bit idealistic. 
Namely, we assumed our mirrors and beamsplitters to be lossless, but it could never come true 
in real experiments; therefore, we need some way to describe losses within the framework of our 
formalism. Optical loss is a term that comprises a very wide spectrum of physical processes, 
including scattering on defects of the coating, absorption of light photons in the mirror bulk and 
coating that yields heating and so on. A full description of loss processes is rather complicated. 
However, the most important features that influence the light fields, coming off the lossy optical 
element, can be summarized in the following two simple statements: 

1. Optical loss of an optical element can be characterized by a single number (possibly, frequency 
dependent) e (usually, |e| <C 1) that is called the absorption coefficient, e is the fraction of 
light power being lost in the optical element: 

£;°"*(i) ^ \/n^£;°"*(t). 



2. Due to the fundamental law of nature summarized by the Fluctuation Dissipation Theorem 
(FDT) [?51 [^S], optical loss is always accompanied by additional noise injected into the 
system. It means that additional noise field h uncorrelated with the original light is mixed 
into the outgoing light field in the proportion of y^e governed by the absorption coefficient. 

These two rules conjure up a picture of an effective system comprising of a lossless mirror and 
two imaginary non-symmetric beamsplitters with reflectivity \/l — e and transmissivity t/e that 
models optical loss for both input fields, as drawn in Figure [7J 
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Figure 7: Model of lossy mirror. 
Using the above model, it is possible to show that for a lossy mirror the transformation of 



carrier fields given by Eq. ( 30 ) should be modified by simply multiplying the output fields vector 
by a factor 1 — e: 



£OUt 

Cout 
^2 



= (1- 



real 



Cin 

On 



real 



Si? 



(33) 



where we used the fact that for low loss optics in use in GW interferometers, the absorption 
coefficient might be as small as e ~ 10~^-10~''. Therefore, the impact of optical loss on classical 
carrier amplitudes is negligible. Where the noise sidebands are concerned, the transformation rule 
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given by Eq. (31 1 changes a bit more: 



'm 

'm 



(1 - e)M„ai 

(1 - e)M,eal 











_n2m_ 




1 - e)M,eal • 



Here, we again used the smallness of e <C 1 and also the fact that matrix Mjcai is unitary, i.e., we 
redefined the noise that enters outgoing fields due to loss as {"ii'^, n'2} = M • {fii, n2} , which 
keeps the new noise sources •h[{t) and n'2{t) uncorrelated: {fii{t)h2{t')) = {n'i{t)h'2{t')) = 0. 

2.2.5 Light modulation by mirror motion 

For full characterization of the light transformation in the GW interferometers, one significant 
aspect remains untouched, i.e., the field transformation upon reflection off the movable mirror. 



Above (see Section 2.1.1), we have seen that motion of the mirror yields phase modulation of the 



reflected wave. Let us now consider this process in more detail. 
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Figure 8: Reflection of light from the movable mirror. 

Consider the mirror described by the matrix Mroah introduced above. Let us set the convention 
that the relations of input and output fields is written for the initial position of the movable mirror 
reflective surface, namely for the position where its displacement is a; = as drawn in Figure [8] We 
assume the sway of the mirror motion to be much smaller than the optical wavelength: x/Aq ^ 1. 
The effect of the mirror displacement x{t) on the outgoing field E°^2{t) can be straightforwardly 
obtained from the propagation formalism. Indeed, considering the light field at a fixed spatial point, 
the reflected light field at any instance of time t is just the result of propagation of the incident 
lightby twice the mirror displacement taken at time of reflection and multiplied by reflectivity 



Er\t) = - 
'it)- 



RE'l\t - 2x{t)/c) + VTE'^\t) , 



rpout ( 



TE'^{t) + yJREl\i + 2x{t)lc) 



(34) 



^In fact, the argument of x should be written as t», that is the moment when the actual reflection takes place 
and is the solution to the equation: c(t — t,) = x(tf)^ but since the mechanical motion is much slower than that of 
light one has lixjc <S 1. This fact implies t ^ tt- 
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Remember now our assumption that x <C Aq; 



according to Eq. ( 19 1 the mirror motion modifies 

It 



the quadrature amphtudes in a way that aUows one to separate this effect from the reflection 
means that the result of the light reflection from the moving mirror can be represented as a sum 
of two independently calculable effects, i.e., the reflection off the fixed mirror, as described above 



in Section 2.2.4 and the response to the mirror displacement (see Section 2.2.1), i.e., the signal 
presentable as a sideband vector {si(i7), S2(^)} . The latter is convenient to describe in terms of 
the response vector {i?i , i?2}^ that is defined as: 



si{n) = Rix{t) = -VRSV[2uJoxit)/c]-£'^ = 



S2in) = R2x{t) = -VRdV [2ujox{t)/c] ■ f ^" = - 
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(35) 



Note that we did not include sideband fields e'"2(^) in the definition of the response vector. In 
principle, sideband fields also feel the motion induced phase shift; however, as far as it depends on 
the product of one very small value of 2u)Qx{t)/c = 47rx(i)/Ao ^ 1 by a small sideband amplitude 
|e']"2(^)| ^ l'^i'2|j the resulting contribution to the final response will be dwarfed by that of the 



classical fields. Moreover, the mirror motion induced contribution ( 35 ) is itself a quantity of the 



same order of magnitude as the noise sidebands, and therefore we can claim that the classical 



amplitudes of the carrier fields are not affected by the mirror motion and that the relations ( 30 ) 
hold for a moving mirror too. However, the relations for sideband amplitudes must be modified. 



In the case of a lossless mirror, relations (31 1 turn 



e°"'(f7) 



= Mrcal 






Ri 
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x{Vt). 



(36) 



where x(r^) is the Fourier transform of the mirror displacement x{t) 

poo 

x{n) = / dtx{ty^K 



It is important to understand that signal sidebands characterized by a vector {si(ri), S2(r2)} , 
on the one hand, and the noise sidebands {ei(ri), 62(^2)} , on the other hand, have the same 
order of magnitude in the real GW interferometers, and the main role of the advanced quantum 
measurement techniques we are talking about here is to either increase the former, or decrease the 
latter as much as possible in order to make the ratio of them, known as the signal-to-noise ratio 
(SNR), as high as possible in as wide as possible a frequency range. 

2.2.6 Simple example: the reflection of light from a perfect moving mirror 

All the formulas we have derived here, though being very simple in essence, look cumbersome and 
not very transparent in general. In most situations, these expressions can be simplified significantly 
in real schemes. Let us consider a simple example for demonstration purposes, i.e., consider the 
reflection of a single light beam from a perfectly reflecting [R = 1) moving mirror as drawn in 
Figure |9] The initial phase 0o of the incident wave does not matter and can be taken as zero. 
Then £™ = Eq and f " = 0. Putting these values into Eq. (30) and accounting for 5" — Oi quite 



reasonably results in the amplitude of the carrier wave not changing upon reflection off the mirror, 
while the phase changes by tt: 
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Since we do not have control over the laser noise, the input light has laser fluctuations in both 
quadratures e™ = {ei",e\"} that are transformed in full accordance with Eq. (31)): 



Again, nothing surprising. Let us see what happens with a mechanical motion induced component 



of the reflected wave: according to Eq. (36), the reflected light will contain a motion-induced signal 
in the s-quadrature: 



s{n) 



2^0 
c 



x{n). 



This fact, i.e., that the mirror displacement that just causes phase modulation of the reflected 
field, enters only the s-quadrature, once again justifies why this quadrature is usually referred to 



as phase quadrature (cf. section 2.2.2) 
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Figure 9: Schematic view of light modulation by perfectly reflecting mirror motion. An initially 
monochromatic laser field i?'"(t) with frequency wq = 2itc/Xo gets reflected from the mirror that 
commits slow (compared to optical oscillations) motion x{t) (blue line) under the action of external 
force G. Reflected the light wave phase is modulated by the mechanical motion so that the 
spectrum of the outgoing field E°^^{uj) contains two sidebands carrying all the information about 
the mirror motion. The left panel shows the spectral representation of the initial monochromatic 
incident light wave (upper plot), the mirror mechanical motion amplitude spectrum (middle plot) 
and the spectrum of the phase- modulated by the mirror motion, reflected light wave (lower plot). 

It is instructive to see the spectrum of the outgoing light in the above considered situation. It 
is, expectedly, the spectrum of a phase modulated monochromatic wave that has a central peak at 
the carrier wave frequency and the two sideband peaks on either sides of the central one, whose 
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shape follows the spectrum of the modulation signal, in our case, the spectrum of the mechanical 
displacement of the mirror x{t). The left part of Figurelolillustrates the aforesaid. As for laser noise, 
it enters the outgoing light in an additive manner and the typical (though simplified) amplitude 
spectrum of a noisy light reflected from a moving mirror is given in Figure |10| 
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Figure 10: The typical spectrum (amplitude spectral density) of the light leaving the interferometer 
with movable mirrors. The central peak corresponds to the carrier light with frequency luq, two 
smaller peaks on either side of the carrier represent the signal sidebands with the shape defined by 
the mechanical motion spectrum x{fl); the noisy background represents laser noise. 



2.3 Basics of Detection: Heterodyne and homodyne readout techniques 

Let us now address the question of how one can detect a GW signal imprinted onto the parameters 
of the light wave passing through the interferometer. The simple case of a Michelson interferometer 
considered in Section [2. 1 . 2| where the GW signal was encoded in the phase quadrature of the light 
leaking out of the signal(dark) port, does not exhaust all the possibilities. In more sophisticated 
interferometer setups that are covered in sections [5] and |6) a signal component might be present 
in both quadratures of the outgoing light and, actually, to different extent at different frequencies; 
therefore, a detection method that allows measurement of an arbitrary combination of amplitude 
and phase quadrature is required. The two main methods are in use in contemporary GW detectors: 
these are homodyne and heterodyne detection [HI I140L II 641 IM] . Both are common in radio- 
frequency technology as methods of detection of phase- and frequency-modulated signals. The 
basic idea is to mix a faint signal wave with a strong local oscillator wave, e.g., by means of a 
beamsplitter, and then send it to a detector with a quadratic non-linearity that shifts the spectrum 
of the signal to lower frequencies together with amplification by an amplitude of the local oscillator. 
This topic is also discussed in Section 4 of the Living Review by Freise and Strain [65 with more 
details relevant to experimental implementation. 



2.3.1 Homodyne and DC readout 

Homodyne readout. Homodyne detection uses local oscillator light with the same carrier fre- 
quency as the signal. Write down the signal wave as: 



S{t) — Sc{t) cos wo^ + Ss{t) sinwoi 



and the local oscillator wave as: 



L(t) — Lc(t) cos wo^ + Ls{t) sinwoi ■ 



(37) 



(38) 
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Signal light quadrature amplitudes Sc^sit) might contain GW signal Gc^sit) as well as quantum 
noise nc^s{t) in both quadratures: 

Sc,s{t) = Gc^s{t) +nc^s{t), 

while the local oscillator is a laser light with classical amplitudes: 

(4°), 4°)) = (ioCos^LO, Lo Sin (bLo), 

where we introduced a homodyne angle 4>lo, and laser noise lc,s(t): 

Note that the local oscillator classical amplitude Lq is much larger than all other signals: 

Lo > max [/c,s, Gc^s, ric^s] ■ 



Let mix these two lights at the beamsplitter as drawn in the left panel of Figure 11 and detect 
the two resulting outgoing waves with two photodetectors. The two photocurrents ii^2 will be 
proportional to the intensities I12 of these two lights: 



«i oc /i ex 



12 oc /2 oc 



(L + sy 



a 



f + Lo{Gc + lc + n^) cos(f>Lo + io(G, + ls + n^) sin^^o + O [GI.JI 



{L - 5)2 _ l; 



= ^ - LoiGc -lc + ric) cos 0LO - Lo{Gs - h + ris) sin0LO + O [Gl.Jl^nl,] 



where A stands for time averaging of A{t) over many optical oscillation periods, which reflects the 
inability of photodetectors to respond at optical frequencies and thus providing natural low-pass 
filtering for our signal. The last terms in both expressions gather all the terms quadratic in GW 
signal and both noise sources that are of the second order of smallness compared to the local 
oscillator amplitude Lq and thus are omitted in further consideration. 
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Figure 11: Schematic view of homodyne readout {left panel) and DC readout (right panel) principle 
implemented by a simple Michelson interferometer. 



In a classic homodyne balanced scheme, the difference current is read out that contains only a 
GW signal and quantum noise of the dark port: 



«ho 



= ii - Z2 oc 2Lo [{Gc + nc) cos 4>LO + [Gs + Us) sin (pLo] 



(39) 
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Whatever quadrature the GW signal is in, by proper choice of the homodyne angle ipLO oug can 
recover it with minimum additional noise. That is how homodyne detection works in theory. 

However, in real interferometers, the implementation of a homodyne readout appears to be 
fraught with serious technical difficulties. In particular, the local oscillator frequency has to be 
kept extremely stable, which means its optical path length and alignment need to be actively 
stabilized by a low-noise control system [53] . This inflicts a significant increase in the cost of the 
detector, not to mention the difficulties in taming the noise of stabilising control loops, as the 
experience of the implementation of such stabilization in a Garching prototype interferometer has 
shown [B1IZ3IZS]- 

DC readout. These factors provide a strong motivation to look for another way to implement 
homodyning. Fortunately, the search was not too long, since the suitable technique has already 
been used by Michelson and Morley in their seminal experiment |lllj . The technique is known 
as DC-readout and implies an introduction of a constant arms length difference, thus pulling the 



interferometer out of the dark fringe condition as was mentioned in Section 2.1.2 The advantage 
of this method is that the local oscillator is furnished by a part of the pumping carrier light 
that leaks into the signal port due to arms imbalance and thus shares the optical path with the 
signal sidebands. It automatically solves the problem of phase- locking the local oscillator and 
signal lights, yet is not completely free of drawbacks. The first suggestion to use DC readout 
in GW interferometers belongs to Fritschel [55] and then got comprehensive study by the GW 
community [US] dMl ISl • 

Let us discuss how it works in a bit more detail. The schematic view of a Michelson interfer- 
ometer with DC readout is drawn in the right panel of Figure [TTJ As already mentioned, the local 
oscillator light is produced by a deliberately-introduced constant difference 5L of the lengths of the 
interferometer arms. It is also worth noting that the component of this local oscillator created by 
the asymmetry in the reflectivity of the arms that is always the case in a real interferometer and 
attributable mostly to the difference in the absorption of the 'northern' and 'eastern' end mirrors 
as well as asymmetry of the beamsplitter. All these factors can be taken into account if one writes 
the carrier fields at the beamsplitter after reflection off the arms in the following symmetric form: 



;,out^^N -^0 ,-, _ ^ _^_ _ ,_, ^r /_N ^ 



^r(i) - --^(1 -e„)coswo(i- 2i„/c) = --=(1 - Ae)coswo(? + AL/c) , 

rp c 

E°''\t) = ^{1 - Ce) cosojoit - 2Le/c) = ^ (1 + Ae) COS wqC? - AL/c) , 

v2 V 2 

where e„^e and L„^e stand for optical loss and arm lengths of the corresponding interferometer 
arms, Ae = 27i^g) ^^ ^^^ optical loss relative asymmetry with e — (e„ -I- ee)/2, £ — Eo{l — e) is the 
mean pumping carrier amplitude at the beamsplitter, AL — Ln~ L^ and t ^ t+ ^"2/^° ■ Then the 
classical part of the local oscillator light in the signal (dark) port will be given by the following 
expression: 

L'^1. (t) = " ' ^ = £ Ae cos — cos ujQt + £ sin — sin uot , (40) 

V2 ^ ^ ^ 

where one can define the local oscillator phase and amplitude through the apparent relations: 



Ae c 



tan0nc = ^tan^^^^ <i - fV(A.)^ + ( ^^^ ) ^ "^^^^^^^ , (41) 
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where we have taken into account that wqAX/c <C 1 and the rather smaU absolute value of the 
optical loss coefficient max [e„, eg] ^ f 0^^ ^ f available in contemporary interferometers. One sees 
that were there no asymmetry in the arms optical loss, there would be no opportunity to change 
the local oscillator phase. At the same time, the GW signal in the considered scheme is confined 
to the phase quadrature since it comprises the time-dependent part of AL and thus the resulting 
photocurrent will be proportional to: 

ZDC oc {L + 5)2 ~ (4°c)' + 24°c('r' + ric) cos(/.dc + 2l'^1{Gs + C' + ".) sin^oc, (42) 
where /°"' denote the part of the input laser noise that leaked into the output port: 

lT\t) c, Z^"A6C0S 5^^^ - IT sin ^^^^ c, l^Ae - ^^^^^ , (43) 

c c c 

c c c 

and ric^s stand for the quantum noise associated with the signal sidebands and entering the inter- 
ferometer from the signal port. 

In the case of a small offset of the interferometer from the dark fringe condition, i.e., for 
wqAL/c = 27rAL/Ao ^ 1, the readout signal scales as local oscillator classical amplitude, which 
is directly proportional to the offset itself: L\^q — 2-1: Eq^. The laser noise associated with the 
pumping carrier also leaks to the signal port in the same proportion, which might be considered 
as the main disadvantage of the DC readout as it sets rather tough requirements on the stability 
of the laser source, which is not necessary for the homodyne readout. However, this problem, is 
partly solved in more sophisticated detectors by implementing power recycling and/or Fabry-Perot 
cavities in the arms. These additional elements turn the Michelson interferometer into a resonant 
narrow-band cavity for a pumping carrier with effective bandwidth determined by transmissivities 
of the power recycling mirror (PRM) and/or input test masses (ITMs) of the arm cavities divided 
by the corresponding cavity length, which yields the values of bandwidths as low as ~ 10 Hz. Since 
the target GW signal occupies higher frequencies, the laser noise of the local oscillator around signal 
frequencies turns out to be passively filtered out by the interferometer itself. 

DC readout has already been successfully tested at the LIGO 40-meter interferometer in Cal- 
tech [IM] and implemented in GEO 600 [gIl[51ISD] and in Enhanced LIGO [5SII1]. It has proven 
a rather promising substitution for the previously ubiquitous heterodyne readout (to be considered 
below) and has become a baseline readout technique for future GW detectors [M] . 

2.3.2 Heterodyne readout 

Up until recently, the only readout method used in terrestrial GW detectors has been the het- 
erodyne readout. Yet with more and more stable lasers being available for the GW community, 
this technique gradually gives ground to a more promising DC readout method considered above. 
However, it is instructive to consider briefly how heterodyne readout works and learn some of the 
reasons, that it has finally given way to its homodyne adversary. 

The idea behind the heterodyne readout principle is the generalization of the homodyne read- 
out, i.e., again, the use of strong local oscillator light to be mixed up with the faint signal light 
leaking out the dark port of the GW interferometer save the fact that local oscillator light frequency 
is shifted from the signal light carrier frequency by JIrf ~ several megahertz. In GW interfer- 
ometers with heterodyne readout, local oscillator light of different than wq frequency is produced 
via phase-modulation of the pumping carrier light by means of electro-optical modulator (EOM) 



before it enters the interferometer as drawn in Figure 12 The interferometer is tuned so that the 



readout port is dark for the pumping carrier. At the same time, by introducing a macroscopic 
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Figure 12: Schematic view of heterodyne readout principle implemented by a simple Michelson 
interferometer. Green lines represent modulation sidebands at radio frequency JIrf and blue dotted 
lines feature signal sidebands 

(several centimeters) offset AL of the two arms, which is known as Schnupp asymmetry |136) , 
the modulation sidebands at radio frequency SIrf appear to be optimally transferred from the 
pumping port to the readout one. Therefore, the local oscillator at the readout port comprises 
two modulation sidebands, L\^ct{t) = L^[t) + L^{t), at frequencies wo + ^rf and ujq — ^KF^ re- 
spectively. These two are detected together with the signal sidebands at the photodetector, and 
then the resulting photocurrent is demodulated with the RF-frequency reference signal yielding an 
output current proportional to GW-signal. 

This method was proposed and studied in great detail in the following works [THl 11361 [Ml [7^ 
FfSl 11061 I118J where the heterodyne technique for GW interferometers tuned in resonance with 
pumping carrier field was considered and, therefore, the focus was made on the detection of only 
phase quadrature of the outgoing GW signal light. This analysis was further generalized to detuned 
interferometer configurations in [411 1140] where the full analysis of quantum noise in GW dual- 
recycled interferometers with heterodyne readout was done. 

Let us see in a bit more detail how the heterodyne readout works as exemplified by a simple 
Michelson interferometer drawn in Figure |12| The equation of motion at the input port of the 
interferometer creates two phase-modulation sideband fields {L^(t) and L_(t)) at frequencies ujq^ 
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n 



RF- 



L+{t) = 
L-{t) = 






cos(wo + f^Rp)^ + 
cos(iujo — f^Rp)^ + 



ii°J + Is 



sin(wo + ^Rp)^ , 
sin ( Wo - ^RF)t , 



(0) 



where L\J^., stand for classical quadrature amplitudes of the modulation (upper and lower) side- 
bands^ and l(c,s)±it) represent laser noise around the corresponding modulation frequency. 

Umike the homodyne readout schemes, in the heterodyne ones, not only the quantum noise 
components n'^'l falling into the GW frequency band around the carrier frequency ujq has to be 
accounted for but also those rallying around twice the RF modulation frequencies ujq ± 2r2RF : 



^hot (t) 



(Gc + <°) coscjot + {Gs + <°) sin Wot 



-n 



c 

wq — 2nRF 



cos(a;o + 2ilRF)t + n 



a;o+20r{p 
s 

cjQ — 20RP 



sin(ajo + 2QYiF)t 



-n~" -""'"' cos(wo - 2rjRF)t + n^«~^^'«'' sin(wo - 2n^Y)t ■ 
The analysis of the expression for the heterodyne photocurrent 

ihct oc (Lhct + Sy = S^ + 1^ + 1J + 2(L++L_)5 + 2L^LZ 



(45) 
(46) 

(47) 



gives a clue to why these additional noise components emerge in the outgoing signal. It is easier 
to perform this kind of analysis if we represent the above trigonometric expressions in terms of 
scalar products of the vectors of the corresponding quadrature amplitudes and a special unit-length 
vector H[(f>] — {coscj), sincj)} , e.g.: 



S-hct = {G + n^j"^ • H[ojot] + nl^ 



'cJn+20R,pi 



H[{ujo + r!R,F)t] + n^ , . H[{ljo - 2nnF)t] 



where G = {Gc, Gg} and n^a ~ {^c , n'^} ■ Another useful observation, provided that coq ^ 
max[rii, fl2\, gives us the following relation: 



H[{ujo + ni)t]H^[iLJo + n2)t] = - 



cos(i7i — ^2)^ — sin(ili — n2)t 
sin(rii — fl2)t cos(rii — 5l2)t 



= -v[{n,-n2)t] 



Using this relation it is straightforward to see that the first three terms in Eq. (47) give DC 



components of the photocurrent, while the fifth term oscillates at double modulation frequency 
2nRF- It is only the term 2{L^ + L^)S that is linear in GW signal and thus contains useful 
information: 



2(L+ + L_)S' ~ Ic{t) cos f^Rpi + Is{t) sin fiRpt -I- {oscillations at frequency SriRp} 



where 



Icit) 
Isit) 



'(0) 



-(0) 



{G + n^„ + n^o+2n„-p) ' • i+' + (G + n^„ -I- n„„-2nRF)"^ ' -^- 



and ay is the 2nd Pauli matrix: 







*In the resonance-tuned case, the phase modulation of the input carrier field creates equal magnitude sideband 
fields as discussed in Section |2.2.2| and these sideband fields are transmitted to the output port thanks to Schnupp 
asymmetry in the same state, i.e., they remain equal in magnitude and reside in the phase quadrature. In detuned 
configurations of GW interferometers, the upper and lower RF-sideband fields are transformed differently, which 
influences both their amplitudes and phases at the readout port. 
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In order to extract the desired signal quadrature the photodetector readout current Zhet is mixed 
with (multiplied by) a demodulation function D{t) = Dq cos(51rf^ + 0_d) with the resulting signal 
filtered by a low-pass filter with upper cut-off frequency A <^ ^rf so that only components 



oscillating at GW frequencies f^cw ^ ^rf appear in the output signal (see Figure 12). 



It is instructive to see what the above procedure yields in the simple case of the Michelson 

(0) 
c+ 



interferometer tuned in resonance with RF-sidebands produced by pure phase modulation: L),_/ = 



lI.2 ~ and i^^ = L^_ = Lq. The foregoing expressions simplify significantly to the following: 

/,(<) = 2Lo [Gs + rC + -^ + -^ j and Lit) = -Lq «"-2^2hf _ „^o+20hp^ _ 

Apparently, in this simple case of equal sideband amplitudes ( balanced heterodyne detection) , only 
single phase quadrature of the GW signal can be extracted from the output photocurrent, which 
is fine, because the Michelson interferometer, being equivalent to a simple movable mirror with 
respect to a GW tidal force as shown in Section |2.1.1| and |2.1.2[ is sensitive to a GW signal only 
in phase quadrature. Another important feature of heterodyne detection conspicuous in the above 
equations is the presence of additional noise from the frequency bands that are twice the RF- 
modulation frequency away from the carrier. As shown in [ITj this noise contributes to the total 
quantum shot noise of the interferometer and makes the high frequency sensitivity of the GW 
detectors with heterodyne readout 1.5 times worse compared to the ones with homodyne or DC 
readout. 

For more realistic and thus more sophisticated optical configurations, including Fabry-Perot 
cavities in the arms and additional recycling mirrors in the pumping and readout ports, the anal- 
ysis of the interferometer sensitivity becomes rather complicated. Nevertheless, it is worthwhile 
to note that with proper optimization of the modulation sidebands and demodulation function 
shapes the same sensitivity as for frequency-independent homodyne readout schemes can be ob- 
tained |41) . However, inherent additional frequency- independent quantum shot noise brought by 
the heterodyning process into the detection band hampers the simultaneous use of advanced quan- 
tum non-demolition (QND) techniques and heterodyne readout schemes significantly. 
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3 Quantum Nature of Light and Quantum Noise 

Now is the time to remind ourselves of the word 'quantum' in the title of our review. Thus far, 
the quantum nature of laser light being used in the GW interferometers has not been accounted 
for in any way. Nevertheless, quantum mechanics predicts striking differences for the variances of 
laser light amplitude and phase fluctuations, depending on which quantum state it is in. Squeezed 
vacuum |163( 11011 11381 l43l l94] injection that has been recently implemented in the GEO 600 
detector and has pushed the high-frequency part of the total noise down by 3.5 dB |152l [5] serves 
as a perfect example of this. In this section, we provide a brief introduction into the quantization 
of light and the typical quantum states thereof that are common for the GW interferometers. 

3.1 Quantization of light: Two-photon formahsm 

From the point of view of quantum field theory, a freely propagating electromagnetic wave can 
be characterized in each spatial point with location vector r = {x,y,z) at time i by a Heisenberg 
operator of an electric field strength E{r, t)r\ The electric field Heisenberg operator of a light 
wave traveling along the positive direction of the ^-axis can be written as a sum of a positive- and 
negative-frequency parts: 

E{x,y,z;t)^u{x,y,z)[E'^+\t)+E^-\t)} , (48) 

where u{x,y,z) is the spatial mode shape, slowly changing on a wavelength A scale, and 

^(+)(i)=^°°^/^«.e— S and ^(-) (i) = [^W(i)] ' . (49) 

Here, A is the effective cross-section area of the light beam, and a^: (dj^) is the single photon 



annihilation (creation) operator in the mode of the field with frequency uj. The meaning of Eq. (49 ) 
is that the travelling light wave can be represented by an expansion over the continuum of harmonic 
oscillators - modes of the electromagnetic field, - that are, essentially, independent degrees of 
freedom. The latter implies the commutation relations for the operators a^j and dj^: 



dc., al, 



2Tr 6 {lu - Lj') , and [d^, d„/] = d]^, d^, =0. (50) 



In GW detectors, one deals normally with a close to monochromatic laser light with carrier 
frequency ujq, and a pair of modulation sidebands created by a GW signal around its frequency 
in the course of parametric modulation of the interferometer arm lengths. The light field coming 
out of the interferometer cannot be considered as the continuum of independent modes anymore. 
The very fact that sidebands appear in pairs implies the two-photon nature of the processes taking 
place in the GW interferometers, which means the modes of light at frequencies a;i_2 = wpirj have 
correlated complex amplitudes and thus the new two-photon operators and related formalism is 
necessary to describe quantum light field transformations in GW interferometers. This was realized 
in the early 1980s by Caves and Schumaker who developed the two-photon formalism [HI US], which 
is widely used in GW detectors as well as in quantum optics and optomechanics. 

One starts by defining modulation sideband amplitudes as 



^Insofar as the light beams in the interferometer can be well approximated as paraxial beams, and the polarization 
of the light wave does not matter in most of the considered interferometers, we will omit the vector nature of the 
electric field and treat it as a scalar field with strength defined by a scalar operator-valued function E{x, y, z, t). 
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and factoring out the oscillation at carrier frequency loq in Eqs. (48), which yields: 
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(51) 



where we denote Cq 



^^ and define functions A±(f7) following 



as 



\±m = 



Wo 



and use the fact that wq 3> J^gw enables us to expand the limits of integrals to wq — > 00. The 
operator expressions in front of e*"^°* in the foregoing Eqs. (51) are quantum analogues to the 



complex amplitude £ and its complex conjugate £* defined in Eqs. (14): 
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Again using Eqs. (14), we can define two-photon quadrature amplitudes as: 
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(52) 



Note that so-introduced operators of two-photon quadrature amplitudes dc^s{t) are Hermitian and 
thus their frequency domain counterparts satisfy the relations for the spectra of Hermitian operator: 

alAt) ^ o,c^s{t) =^ al.,{n)=dc,si-^)- 
Now we are able to write down commutation relations for the quadrature operators, which can be 



derived from Eq. ( 50 ) 



[ac(f^), al{n')] = [a,{n), dl{n')] = 2TTS{n - n') 

[adn), dl{n')] = [al{n), a,{n')] = 2Tri6{n - n') , 



(53) 
(54) 



The commutation relations represented by Eqs. (53) indicate that quadrature amplitudes do 



not commute at different times, i.e., [dc(t), dc{t')] — [ds{t), as{t')] 7^ 0, which imply they could 
not be considered for proper output observables of the detector, for a nonzero commutator, as we 
would see later, means an additional quantum noise inevitably contributes to the final measurement 
result. The detailed explanation of why it is so can be found in many works devoted to continuous 
linear quantum measurement theory, in particular, in Chapter 6 of |27j . Appendix 2.7 of |48| or 
in [21]. Where GW detection is concerned, all the authors are agreed on the point that the values of 
GW frequencies 51 (1 Hz ^ 51/27r ^ 10'^ Hz), being much smaller than optical frequencies a;o/27r ^ 



10 Hz, allow one to neglect such weak commutators as those of Eqs. ( 53 ) in all calculations related 



to GW detectors output quantum noise. This statement has gotten an additional ground in the 
calculation conducted in Appendix 2.7 of [48^ where the value of the additional quantum noise 
arising due to the nonzero value of commutators ( 53 ) has been derived and its extreme minuteness 



compared to other quantum noise sources has been proven. Braginsky et al. argued in '53' that 
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the two-photon quadrature aniphtudes defined by Eqs. (52) are not the real measured observables 



at the output of the interferometer, since the photodetectors actually measure not the energy flux 

'•°° dujduj' ^ 
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but rather the photon number flux: 

N{t) 







{2^f 



' UJUJ 



'at,a„,e*("-')* 



dujdw' 



(27r)2 """" 






(55) 



(56) 



The former does not commute with itself: [i{t), i{t')] ^ 0, while the latter apparently does 
[M{t), A/'(t')] — and therefore is the right observable for a self-consistent quantum description of 
the GW interferometer output signal. 

In the course of our review, we shall adhere to the approximate quadrature amplitude operators 



that can be obtained from the exact ones given by Eqs. (52) by setting A±(ri) — > 1, i.e.. 



ac{t) = 



dfl 0.4 



-oo2^ \/2 
°° dna+- dl 
-oo 27r iy^ 



„-iQ.t 



-Hit 



dc{D,) = 



V2 



dsi^) 



a+ — d_ 



iV2 



(57) 



The new approximate two-photon quadrature operators satisfy the following commutation re- 
lations in the frequency domain: 

[ac(r2), ds{n')] = 27:iS{n + n') , and [dc{n), dc{n')] = [d,{n), a,(f7')] = > (58) 

and in the time domain: 

[d,{t), a,(i')] = iS{t - t') , and [a,{t), d,{t')] = [ds{t), a,{t')] = (59) 



Then the electric field strength operator ( 48 ) can be rewritten in terms of the two-photon 
quadrature operators as: 



E{x, y, z; t) — u{x, y, z)Co [dc{t) cosluqI + ds{t) sinojo^] 



(60) 



Hereafter, we will omit the spatial mode factor u(a;, y, z) since it does not influence the final result 
for quantum noise spectral densities. Moreover, in order to comply with the already introduced 
division of the optical field into classical carrier field and to the 1st order corrections to it comprising 
of laser noise and signal induced sidebands, we adopt the same division for the quantum fields, i.e., 
we detach the mean values of the corresponding quadrature operators via the following redefinition 



iold 



^ An 



with A, 



c,s — \"c 



Here, by (a° g) we denote an ensemble average over the 
quantum state \ip) of the light wave: (A) = {ip\A\tp). Thus, the electric field strength operator for 
a plain electromagnetic wave will have the following form: 



E{x, y; t) = Co [{Ac + dc{t)) cos wq* + {A^ + as{t)) sinwoi] 



(61) 



Further, we combine the two-photon quadratures into vectors in the same manner as we used to 
do for classical fields: 

'Ac 



and d ; 
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Now, when we have defined a quantum Heisenberg operator of the electric field of a light wave, 
and introduced quantum operators of two-photon quadratures, the last obstacle on our way towards 
the description of quantum noise in GW interferometers is that we do not know the quantum state 
the light field finds itself in. Since it is the quantum state that defines the magnitude and mutual 
correlations of the amplitude and phase fluctuations of the outgoing light, and through it the total 
level of quantum noise setting the limit on the future GW detectors' sensitivity. In what follows, we 
shall consider vacuum and coherent states of the light, and also squeezed states, for they comprise 
the vast majority of possible states one could encounter in GW interferometers. 

3.2 Quantum states of light 

3.2.1 Vacuum state 

The quantum state of the travelling wave is a subtle structure, for the system it describes comprises 
a continuum of modes. However, each of these modes can be viewed at as a quantum oscillator 
with its own generalized coordinate X^ = {a^ +aj,)/v2 and momentum Y^ = (a^^ — aJ,)/z\/2. The 
ground state of this system, known as a vacuum state |vac), is straightforward and is simply the 
direct product of the ground states \0)^ of all modes over all frequencies w: 

I vac) = (g) |0)^ . (62) 

By definition, the ground state of a mode with frequency uj is the state with minimum energy 
Evac = fiuj/2 and no excitation: 

a„|0)^ = and (0|„al,=0. (63) 

Consider now statistical properties of the vacuum state. The mean values of annihilation and 
creation operators as well as any linear combination thereof that includes quadrature amplitudes, 
are zero: 

(vac|a<^|vac) = (a„) = (4) = ^ {ad^)) = (as(^)) =0. 

Apparently, this also holds for time domain operators: 

(vac|a(t)|vac) = (a(i)) = {a\t)) = ^ (a,(i)) = (a,(t)) =0. 
That the ground state of the oscillator is Gaussian is evident from its q-representation [138) , namely 

V'vac(^^) = (^^|0)^ - ^ exp |-^ 



It means that knowing the second moments of quadrature amplitudes suffices for full character- 
ization of the state |vac). For this purpose, let us introduce a quadrature amplitudes matrix of 
spectral densitie3j§(ri) according to the rule: 



{ac{n) o ac{n')) {a,{n) o as{n')) 

{asin)oac{n')) {a^in) o dsin')) 



= 2nS{n + n') 






= 2TT^5{n + n') . (64) 



"Herein, we make use of a double-sided power spectral density defined on a whole range of frequencies, both 
negative and positive, that yields the following connection to the variance of an arbitrary observable d{t): 



Var [o{i)] = (62{t)>-(6{t)>^ 






It is worth noting that in the GW community, the sensitivity of GW detectors as well as the individual noise 
sousces are usually characterized by a single-sided power spectral density 5^ (C), that is simply defined on positive 
frequencies H ^ 0. The connection between these two is straightforward: So (f2) = 2So{fi) for f2 ^ and otherwise. 
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where Sij{fl) {i,j — c,s) denote (cross) power spectral densities of the corresponding quadrature 
amphtudes {ai{fl) o dj{il')) standing for the symmetrized product of the corresponding quadrature 
operators, i.e.: 

For a vacuum state, this matrix of spectral densities can easily be obtained from the commutation 

"1/2 



relations (581 and equals to: 



iSivac(' 'j 









1/2 



(65) 



which implies that the (double-sided) power spectral densities of the quadrature amplitudes as well 
as their cross-spectral density are equal to: 



Sccrn = Sssin) 



and 



Scs{n)^0. 



In time domain, the corresponding matrix of second moments, known as a covariance matrix with 
elements defined as YijS{t — t') — {ai{t) o aj{t')), is absolutely the same as Svac(^) : 



1/2 
1/2 



(66) 



It is instructive to discuss the meaning of these matrices, S and V, and of the values they 
comprise. To do so, let us think of the light wave as a sequence of very short square-wave light 
pulses with infinitesimally small duration e -> 0. The delta function of time in Eq. (66 1 tells us that 



the noise levels at different times, i.e., the amplitudes of the different square waves, are statistically 
independent. To put it another way, this noise is Markovian. It is also evident from Eq. ( |65[ ) that 
quadrature amplitudes' fluctuations are stationary, and it is this stationarity, as noted in @J that 
makes quadrature amplitudes such a convenient language for describing the quantum noise of light 
in parametric systems exemplified by GW interferometers. 

It is instructive to pay some attention to a pictorial representation of the quantum noise de- 
scribed by the covariance and spectral density matrices V and S. With this end in view let us 
introduce quadrature operators for each short light pulse as follows: 



Xe{t) 



1 



t+e/2 



t-e/2 



dTadr) , 



and 



Ut) ^ 



1 



t+e/2 



t-6/2 



dTCLsi^r) . 



(67) 



These operators X^ (i) and Y^ (i) arc nothing else than dimensionless displacement and momen- 
tum of the corresponding mode (called quadratures in quantum optics), normalized by zero point 
fluctuation amplitudes Xo and Po- X^{t) = £e/Xo and X^{t) = p^jP^. This fact is also justified 
by the value of their commutator: 

X,(i),n(t)l =i. 



There is no difficulty in showing that diagonal elements of the covariance matrix 
the variances of the corresponding mode displacement X^ and momentum Y^ : 



are equal to 



(X2(t)) = l/2, and 



(n'W> = l/2, 



while non-diagonal terms represent correlations between these operators (zero in our case): N cs = 
Vsc = (-^e(0°^e(^)) = 0- ^t the same time, we see that there is no correlation between the pulses, 
justifying the Markovianity of the quantum noise of light in vacuum state: 

{x,{t)x,(s!)) = {%{t)%{^)) = {x,{t) o n(i')) = , t ^ t' . 
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An attempt to measure the light field amplitude as a function of time will give the result depicted 
in Figure [TS] 



(a) 



vacuum state 



(b) 



nj 




- . o ■. ' 
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Figure 13: Light field in a vacuum quantum state |vac). Left panel (a) features a typical result 
one could get measuring the (normalized) electric field strength of the light wave in a vacuum 
state as a function of time. Right panel (b) represents a phase space picture of the results of 
measurement. A red dashed circle displays the error ellipse for the state |vac) that encircles 
the area of single standard deviation for a two-dimensional random vector a of measured light 
quadrature amplitudes. The principal radii of the error ellipse (equal in vacuum state case) are 
equal to square roots of the covariance matrix Vvac eigenvalues, i.e., to l/v2. 



The measurement outcome at each instance of time will be a random variable with zero mean 
and variance defined by a covariance matrix Vvac of Eq. ( 66 ) : 



Var[ij(i)] = {coswo^i sinajo0^vac{coswoi, sinwo^} — ^ • 

In quantum mechanics, it is convenient to describe a quantum state in terms of a Wigner 
Junction, a quantum version of joint (quasi) probability distribution for particle displacement and 
momentum (X^ and Y^ in our case): 



W^|vac)(^s,n) = 



1 



27^^/de^ 



exp {-i^Y,} {X, + e/2|vac) (vac|X, - ^/2) 



: exp \ -i{X„ Y,yY-,\{X,, yA = 1 exp {-{X^ + y/)|Q8) 



where ^ is simply the variable of integration. The above Wigner function describes a Gaussian state, 
which is simply the ground state of a harmonic oscillator represented by a mode with dis placement 

Gaussian 



X^ and momentum Fg . The corresponding plot is given in the left panel of Figure 

states are traditionally pictured by error ellipses on a phase plane, as drawn in the right panel 



14 



of Figure 14 (cf. right panel of Figure 13). Here as well as in Figure 13 a red line in both plots 



circumscribes all the values of X^ and Y^ that fall inside the standard deviation region of the Wigner 
function, i.e., the region where all pertinent points are within 1 standard deviation from the center 
of the distribution. For a vacuum state, such a region is a circle with radius ^/Ycc — V^^Ts = 1/ v2. 
The area of this circle, equal to 1/2 in dimensionless units and to h/2 in case of dimensional 
displacement and momentum, is the smallest area a physical quantum state can occupy in a phase 
space. This fact yields from a very general physical principle, the Heisenberg uncertainty relation, 
that limits the minimal uncertainty product for canonically conjugate observables (displacement 
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Figure 14: Wigner function W\vac){Xs, Ye) of a ground state of harmonic oscillator {left panel) and 
its representation in terms of the noise ellipse (right panel). 

X^ and momentum Y^, in our case) to be less than 1/2 in /i- units: 



1/2 



(Var[l,]) "(varK 



1/2 1 

> - 

2 



The fact that for a ground state this area is exactly equal to 1/2 is due to the fact that it is a 
pure quantum state, i.e., the state of the particle that can be described by a wave function |-0), 
rather than by a density operator p. For more sophisticated Gaussian states with a non-diagonal 
covariance matrix V, the Heisenberg uncertainty relation reads: 



detV> 



1 

4' 



(69) 



and noise ellipse major semi-axes are given by the square root of the matrix V eigenvalues. 

Note the difference between Figures [13] and [M] the former features the result of measurement of 
an ensemble of oscillators (subsequent light pulses with infinitesimally short duration e), while the 
latter gives the probability density function for a single oscillator displacement and momentum. 



3.2.2 Coherent state 



Another important state of light is a coherent state (see, e.g., [1631 11381 11011 [TM] ). It is straight- 
forward to introduce a coherent state \a) of a single mode or a harmonic oscillator as a result of 
its ground state |0) shift on a complex plane by the distance and in the direction governed by a 
complex number a — |a|e*^''s("). This can be caused, e.g., by the action of a classical effective 
force on the oscillator. Such a shift can be described by a unitary operator called a displacement 
operator, since its action on a ground state |0) inflicts its shift in a phase plane yielding a state 
that is called a coherent state: 



\a) ^D[a]\Q) ee e° 



'|0), 



or, more vividly, in q-representation of a corresponding mode of the field [134] 



^,oh{Xu) = {X^\a) 



: exp ■ 



(X^ - V2a)^ 
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The shift described by D[a\ is even more apparent if one writes down its action on an annihilation 
(creation) operator: 

D^a]aD[a] ^ a + a , fl)^ [a]a^ £>[«]= a^ + a*^ . 

Moreover, a coherent state is an eigenstate of the annihilation operator: 

a\a) = a\a) . 

Using the definitions of the mode quadrature operators X^ = X and Y^^ = Y (dimensionless 
oscillator displacement and momentum normalized by zero-point oscillations amplitude) given 
above, one immediately obtains for their mean values in a coherent state: 

{a\X\a) = (0|Z)^HXi)[a]|0) = V2Re[a] , {a\Y\a) = {0\D''[a]YD[a]\0) = V2lm[a] . 

Further calculation shows that quadratures variances: 

Var[X] = {a\X^\a) - ((a|X|a))' = ^ , Var[y] = {a\Y^\a) - ((a|r|a))' = ^ 

have the same values as those for a ground state. These two facts unequivocally testify in favour 
of the statement that a coherent state is just the ground state shifted from the origin of the phase 
plane to the point with coordinates {{X)a, {y)a) — V2(Re[a] ,Re[a]). It is instructive to calculate 
a Wigner function for the coherent state using a definition of Eq. ([68| : 



VF|„) (X, Y)^ f ^ cxp {-iCY} {X + e/2|a) {a\X - ^/2) = 

^^^^^^^ exp |-^{^ - V2Re[a], Y - V2lm[a]yY-^AX - \/2Re[a], Y - V2lm[a]}\ = 

- exp |- {X - V2Re[a]f + {Y - V2Im[a]f] \ , 

which once again demonstrates the correctness of the former statement. 

Generalization to the case of continuum of modes comprising a light wave is straightforward [T9] 



and goes along the same lines as the definition of the field vacuum state, namely (see Eq. (62)): 
|a(c.)) ^ (g) \a)^ = (^D[a{u;)]\0)^ = exp j T ^(a(c.)ai - a*{u)aA |vac) , (70) 

where \a)i^ is the coherent state that the mode of the field with frequency to is in, and a(w) is 
the distribution of complex amplitudes a over frequencies w. Basically, a{cu) is the spectrum of 
normalized complex amplitudes of the field, i.e., a{Ld) oc £{10). For example, the state of a free light 
wave emitted by a perfectly monochromatic laser with emission frequency ujp and mean optical 

power Xq will be defined by a{uj) — tta/j-j-S{uj—uJp), which implies that only the mode at frequency 
ujp will be in a coherent state, while all other modes of the field will be in their ground states. 

Operator D[a] is unitary, i.e., £''''[a]£>[a] = I)[a]-D^[a] = / with / the identity operator, while 
the physical meaning is in the translation and rotation of the Hilbert space that keeps all the 
physical processes unchanged. Therefore, one can simply use vacuum states instead of coherent 
states and subtract the mean values from the corresponding operators in the same way we have 



done previously for the light wave classical amplitudes, just below Eq. (60 1. The covariance matrix 
and the matrix of power spectral densities for the quantum noise of light in a coherent state is thus 
the same as that of a vacuum state case. 

The typical result one can get measuring the electric field strength of light emitted by the 



aforementioned ideal laser is drawn in the left panel of Figure 15 
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Figure 15: Light field in a coherent quantum state \a{uj)). Left panel a) features a typical result one 
could get measuring the (normalized) electric field strength of the light wave in a coherent state as 
a function of time. Right panel b) represents a phase space picture of the results of measurement. 
The red dashed line in the left panel marks the mean value {E{t)). The red arrow in the right panel 
features the vector of the mean values of quadrature amplitudes, i.e., A, while the red dashed circle 
displays the error ellipse for the state |a(a;)) that encircles the area of single standard deviation 
for a two-dimensional random vector a of quadrature amplitudes. The principle radii of the error 
ellipse (equal in the coherent state case) are equal to square roots of the covariance matrix Vcohj 
i.e., to 1/^2. 

3.2.3 Squeezed state 

One more quantum state of light that is worth consideration is a squeezed state. To put it in simple 
words, it is a state where one of the oscillator quadratures variance appears decreased by some factor 
compared to that in a vacuum or coherent state, while the conjugate quadrature variance finds itself 
swollen by the same factor, so that their product still remains Heisenberg-limited. Squeezed states 
of light are usually obtained as a result of a parametric down conversion (PDC) process [5S1 1172] 
in optically nonlinear crystals. This is the most robust and experimentally elaborated way of 
generating squeezed states of light for various apphcations, e.g., for GW detectors |150|.ll53|ll43j . 
or for quantum communications and computation purposes 36J . However, there is another way to 
generate squeezed light by means of a ponderomotive nonlinearity inherent in such optomechanical 
devices as GW detectors. This method, first proposed by Corbitt et al. [53^, utilizes the parametric 
coupling between the resonance frequencies of the optical modes in the Fabry-Perot cavity and the 
mechanical motion of its mirrors arising from the quantum radiation pressure fluctuations inflicting 
random mechanical motion on the cavity mirrors. Further, we will see that the light leaving the 
signal port of a GW interferometer finds itself in a ponderomotively squeezed state (see, e.g., [M] 
for details). A dedicated reader might find it illuminating to read the following review articles on 
this topic [raSlfTM] . 

Worth noting is the fact that generation of squeezed states of light is the process that inherently 
invokes two modes of the field and thus naturally calls for usage of the two-photon formalism 
contrived by Caves and Schumaker |331 US] . To demonstrate this let us consider the physics of 
a squeezed state generation in a nonlinear crystal. Here photons of a pump light with frequency 
ujp — 2ujQ give birth to pairs of correlated photons with frequencies wi and uj2 (traditionally called 
signal and idler) by means of the nonlinear dependence of polarization in a birefringent crystal 
on electric field. Such a process can be described by the following Hamiltonian, provided that 
the pump field is in a coherent state \a)ujp with strong classical amplitude \ap\ ^ 1 (see, e.g.. 
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Section 5.2 of [163] for details): 



Hpuc 

h 



uJia[ai + UJ20'2'^2 + i I X(H^2^ 



-2iujQt 



X oiflae 



2iuj()t 



)■ 



(71) 



where ai^2 describe annihilation operators for the photons of the signal and idler modes and 
X = yoe^"^ is the complex coupling constant that is proportional to the second-order susceptibility 
of the crystal and to the pump complex amplitude. Worth noting is the meaning of t in this 
Hamiltonian: it is a parameter that describes the duration of a pump light interaction with the 
nonlinear crystal, which, in the simplest situation, is either the length of the crystal divided by the 
speed of light c, or, if the crystal is placed between the mirrors of the optical cavity, the same as 
the above but multiplied by an average number of bounces of the photon inside this cavity, which 
is, in turn, proportional to the cavity finesse T. It is straightforward to obtain the evolution of 
the two modes in the interaction picture (leaving apart the obvious free evolution time dependence 



*) solving the Heisenberg equations: 



ai{t) 



di cosh pt + ale^^'^ sinh pt , a2{t) — 0,2 cosh pt + a\e'^^'^ sinh pt . 



-^ J2i 



(72) 



Let us then assume the signal and idler mode frequencies symmetric with respect to the half of 
pump frequency wq = a;p/2: uji — > lo^ ~ ujq -\- Q. and 102 — > w_ = ujq — Q. (di — > d+ and 02 — > d_). 
Then the electric field of a two-mode state going out of the nonlinear crystal will be written as (we 
did not include the pump field here assuming it can be ruled out by an appropriate filter): 

E{t) = Co \x{t) cos Wot -I- Y{t) sinwotj , 



where two-mode quadrature amplitudes X{t) and Y{t) are defined along the lines of Eqs. (57 1 
keeping in mind that only idler and signal components at the frequencies u}± — wq 
be kept in the integral, which yields: 



±VL should 



X{t) 

Y{t) 

where d^'^^ 



1 

V2 
1 

iV2 



d+(t)e'"*-d^f)e-*"* 
(d+(i) + al{t))/V2 and d|'5^ 



d_(t)e-''"*- 
-d_(t)e-'"* 



-dl(t)e*"* 



iQt I p-sqziifit 
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sqz —iQ,t 



^sqztg^nt 



(d+(t) — a_{t))/{i\/2) are the spectral quadrature 
amplitudes of the two-mode field at sideband frequency fl (cf. Eqs. ([52])) after it leaves the nonlinear 
crystal. Substituting Eqs. ( |72[ ) into the above expressions yields transformation rules for quadrature 
amplitudes: 



^sqz 



dc (cosh pt + cos 2(f> sinh pt) + ds sin 20 sinh pt , 
dc sin 2(j) sinh pt + d^ (cosh pt — cos 2(f> sinh pt) , 



(73) 



where dc = (d+ + a_)/\/2 and d^ = (04- — dl)/(iv2) stand for initial values of spectral quadra- 
ture amplitudes of the two-mode light wave created in the PDC process. A close look at these 
transformations written in the matrix form reveals that it can be represented as the following 
sequence: 



~sqz 



= §sqz[pt, 0]a - P[0]§sqz[pt, 0]P[~4>]a 



(74) 



where 



cosh pt + cos 20 sinh pt sin 20 sinh pt 

sin 20 sinh pt cosh pt — cos 20 sinh pt 



J'sqz 



[pt,0] = 



r,pt 



'pt 



(75) 
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are squeezing matrices in general and in special (0 = 0) case, while P[0] stands for a counterclock- 



wise 2D-rotation matrix by angle 4> defined by (17 1. Therefore, the evolution of a two-mode light 



quadrature amplitude vector a in a PDC process described by the Hamiltonian ( 71 ) consists of a 



clockwise rotation by an angle (f) followed by a deformation along the main axes (stretching along 
the flc-axis and proportional squeezing along the a^-axis) and rotation back by the same angle. It 

is straightforward to show that vector X^^"^ ~ \^{t), Y{t) \ = a^^^e^^^^ + asqz*gint transforms 

similarly (here a^^qz* __ | 



sqzf Asqzf 



'}')■ 



This geometric representation is rather useful, particularly for the characterization of a squeezed 
state. If the initial state of the two-mode field is a vacuum state then the outgoing field will be in 
a squeezed vacuum state. One can define it as a result of action of a special squeezing operator 
S[pt, (j)] on the vacuum state 

|sqzo(pt, (J))) = S[pt, (I)] |vac) . (76) 

This operator is no more and no less than the evolution operator for the PDC process in the 
interaction picture, i.e.. 



S[pt,4>] = exp|pi(a+a_e-2*'^ - a^a^.e^^"^)} . (77) 

Action of this operator on the two-photon quadrature amplitudes is fully described by Eqs. ( |76[ ): 

d'^i' = S^pt,<j>]aS[pt,cf>] = P[0]§,q,[pi,O]P[-0]d, 
while annihilation operators of the corresponding modes a± are transformed in accordance with 



Eqs. (72) 



The linearity of the squeezing transformations implies that the squeezed vacuum state is Gaus- 
sian since it is obtained from the Gaussian vacuum state and therefore can be fully characterized 
by the expectation values of operators X and Y and their covariance matrix Vgqz- Let us calculate 
these values: 

(X)3q, = (sqzo(r,,/))|X|sqzo(r,</))) = (vac|^t [r, 0]l^[r, (/.]|vac) = (vac|X(t)|vac) = 0, (y)sqz = 0, 

and for a covariance matrix one can get the following expression: 

V,q, = (sqZo(r, <P)\X o XT|sqZo(r, 0)) = P[-0]§3qzh 0]Vv 



Ssqz [2^ 



coscp 
— sini 



smt 

COS( 







-2r 



vac'^sqz L 
COS(j) 

sind) 



— Sine 
cos 6 



(78) 



where we introduced squeezing parameter r = pt and used a short notation for the symmetrized 
outer product of vector X with itself: 



XoX~^ 



XoX XoY 

YoX YoY 



The squeezing parameter r is the quantity reflecting the strength of the squeezing. This way 
of characterizing the squeezing strength, though convenient enough for calculations, is not very 
ostensive. Conventionally, squeezing strength is measured in decibels (dB) that are related to the 
squeezing parameter r through the following simple formula: 



TdB = lOlogioe^'' = 20rlogioe <^=^ r = rdB/(201ogio e) . 
For example, 10 dB squeezing corresponds to r ~ 1.15. 



(79) 
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a) Initial state 
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Figure 16: Schematic plot of a vacuum state transform.ation under the action of the squeezing 
operator ^[r, 0]. Eqs. (74) demonstrate the equivalence of the general squeezing operator 5'[r, (/)] 
to a sequence of phase plane counterclockwise rotation by an angle (f> (transition from a) to b)), 
phase plane squeezing and stretching by a factor e^ (transition from b) to c)) and rotation back 
by the same angle (p (transition from c) to d)). Point P' tracks how transformations change the 
initial state marked with point P. 



The covariance matrix (781 refers to a unique error ellipse on a phase plane with semi-major 
axis e''/v2 and semi-minor axis e" 
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'"/"v/S rotated by angle clockwise as featured in Figure 
It would be a wise guess to make, that a squeezed vacuum Wigner function can be obtained 
from that of a vacuum state, using these simple geometric considerations. Indeed, for a squeezed 
vacuum state it reads: 



W^|sqz>(^, y) 



1 



27r^detVs, 



: exp 



1 



{X, YyW;'{X, Y} 



(80) 



where the error ellipse refers to the level where the Wigner function value falls to 1/y/e of the max- 
imum. The corresponding plot and phase plane picture of the squeezed vacuum Wigner function 



are featured in Figure 17 



Another important state that arises in GW detectors is the displaced squeezed state |sqz^(r, (j>)) 
that is obtained from the squeezed vacuum state in the same manner as the coherent state yields 
from the vacuum state, i.e., by the application of the displacement operator (equivalent to the 
action of a classical force) : 



|sqz„(r,(/))) = Z)[a]|sqzo(r,0)) = D[a]S[r, 



(81) 



The light leaving a GW interferometer from the signal port finds itself in such a state, if a classical 
GW-like force changes the difference of the arm lengths, thus displacing a ponderomotively squeezed 
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Figure 17: Left panel: Wigner function of a squeezed vacuum state with squeeze parameter r = 0.5 
(5 dB) and rotation angle <j> — ~7i'/4. Right panel: Error ellipse corresponding to that Wigner 
function. 



vacuum state in phase quadrature Y by an amount proportional to the magnitude of the signal 
force. Such a displacement has no other consequence than simply to shift the mean values of X 
and Y by some constant values dependent on shift complex amplitude a: 

(X)sqz = V2Re[a] , (Y),^, = V2lm[a] . 

Let us now generalize the results of a two-mode consideration to a continuous spectrum case. 
Apparently, quadrature operators X{t) and Y{t) are similar to adt) and ds(t) for the traveling 
wave case. Utilizing this similarity, let us define a squeezing operator for the continuum of modes 
as: 



27r 



a+a^e 



-2j0(n) 



-t ;.t p2i0(O) 



(82) 



where r(D,) and (f>(S^) are frequency-dependent squeezing factor and angle, respectively. Acting 
with this operator on a vacuum state of the travelling wave yields a squeezed vacuum state of 



a continuum of modes in the very same manner as in Eq. (76). The result one could get in the 



measurement of the electric field amplitude of light in a squeezed state as a function of time is 
presented in Figure [Ts} Quadrature amplitudes for each frequency ft transform in accordance with 
Eqs. (73). Thus, we are free to use these formulas for calculation of the power spectral density 



matrix for a traveling wave squeezed vacuum state. Indeed, substituting ac^si^) — ^ '^cT(^) ^^ 



Eq. ( 64 ) and using Eq. ( 74 ) one immediately gets 



Ssq.lf^) = Ph(/.(r!)]S3qz[K^),O]§vac(^)§sqz[r(r!),O]P[0(f^)] - §,qz(r(f2),0(f])) 



(83) 



Note that entries of §sqz(^) might be frequency dependent if squeezing parameter r(D,) and squeez- 
ing angle 4'i^) ^-re frequency dependent as is the case in all physical situations. This indicates 
that quantum noise in a squeezed state of light is not Markovian and this can easily be shown by 
calculating the the covariance matrix, which is simply a Fourier transform of S(r2) according to 
the Wiener-Khinchin theorem: 



' sqz 



it-t') 



,2^ 



Ssqz(^)e~'"(*-*') 



(84) 
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h) amplitude squeezing (10 clB, 0=— , ^^=5) c) phase squeezing (10 dB, 0=0, /1,,=5) 
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Figure 18: Light field in a squeezed state |sqz^(r, 0)). Upper row features time dependence of the 
electric field strength E{t) in three different squeezed states (10 dB squeezing assumed for all): a) 
squeezed vacuum state with squeezing angle 4> = 7r/4; b) displaced squeezed state with classical 
amplitude Ac = h (mean field strength oscillations (£^(i))sqz are given by red dashed line) and 
amplitude squeezing (</) = 7r/2); c) displaced squeezed state with classical amplitude Ac = h and 
phase squeezing (0 = 0). Lower row features error ellipses (red dashed lines) for the corresponding 
plots in the upper row. 



Of course, the exact shape of Vsqz(t — i') could be obtained only if we specify r{Vt) and 0(^2). Note 
that the noise described by Vsqz(i — t') is stationary since all the entries of the covariance matrix 
(correlation functions) depend on the difference of times t — t' . 

The spectral density matrix allows for pictorial representation of a multimode squeezed state 
where an error ellipse is assigned to each sideband frequency fl. This effectively adds one more 
dimension to a phase plane picture already used by us for the characterization of a two-mode 
squeezed states. Figure [19] exemplifies the state of a ponderomotively squeezed light that would 



leave the speedmeter type of the interferometer (see Section 6.2 1 



3.3 How to calculate spectral densities of quantum noise in linear optical 
measurement? 

In this section, we give a brief introduction to calculation of the power spectral densities of quantum 
noise one usually encounters in linear optical measurement. In optomechanical sensors, as we have 
discussed earlier, the outgoing light carries the information about the measured quantity (e.g., the 
displacement due to GW tidal forces) in its phase and (sometimes) amplitude quadratures. The 
general transformation from the input light characterized by a vector of quadrature amplitudes 
a{n) = {(ici^), as(D,)} to the readout quantity of a meter is linear and can be written in spectral 
form as: 

Y{n) = y'f{n)d{n) + Gin) = y*cin)acin) + y;{n)as{n) + Gin), (85) 
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Figure 19: Example of a squeezed state of the continuum of modes: output state of a speedmeter 
interferometer. Left panel shows the plots of squeezing parameter rdB(^) and squeezing angle 4'{^) 
versus normalized sideband frequency 0/7 (here 7 is the interferometer half-bandwidth). Right 
panel features a family of error ellipses for different sideband frequencies Q that illustrates the 
squeezed state defined by rdB(^) and (/"(f^) drawn in the left panel. 

where G{Q) is the spectrum of the measured quantity, y^si^) are some complex- valued functions 
of fi that characterize how the light is transformed by the device. Quantum noise is represented 
by the terms of the above expression not dependent on the measured quantity G, i.e., 



Ny{n)^yHn)d(n) = (y*{n) y*,{n) 




(86) 



The measure of quantum noise is the power spectral density Syi^) that is defined by the 
following expression: 



27rS'y(rj)(5(^ - n') = {^|J\NY{n) o 7V|(f7')|V') = {Ny{^) o N^in')) 



(87) 



Here ji/') is the quantum state of the light wave. 

In our review, we will encounter two types of quantum states that we have described above, 
i.e., vacuum |vac) and squeezed vacuum |sqzg(r, 0)) states. Let us show how to calculate the power 
(double-sided) spectral density of a generic quantity F(il) in a vacuum state. To do so, one should 



substitute Eq. (86) into Eq. (87) and obtain that: 

snn) = y\n){a{n)oa^{n)),,,y{n)=y\n)s,,,in)y{n) 

y\n)y{n) ^ |yc(^)P + |y.(^)P 
2 2 

where we used the definition of the power spectral density matrix of light in a vacuum state ( 65 ^ 
Similarly, one can calculate the spectral density of quantum noise if the light is in a squeezed state 



sqzQ(r, 0)), utilizing the definition of the squeezed state density matrix given in Eq. (83) 

1 



sp'in) 



yHn){d{n) o a^n)),^,y{n)^y\n)Ssq.myin) = -y^n)v[-(t>ps^,[2r,o]F[^]y{fi) 



' ^^ ^' (cosh 2r + sinh 2r cos 2(/)) 
-Re [ycin)y;{n)] sinh2r sin2(/) . 



13^.(^)1' 



(cosh 2r — sinh 2r cos 20) 



(89) 



^Hereafter wc will omit, for the sake of brevity, the factor 27r5(f2 — $1') in equations that define the power (double- 
sided) spectral densities of relevant quantum observables, as well as assume Q = Q' , though keeping in mind that a 
mathematically rigorous definition should be written in the form of Eq. l|87f . 
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It might also be necessary to calculate also cross-correlation spectral density Syzi^) of Y{^) 
with some other quantity Z{il.) with quantum noise defined as: 

Using the definition of cross-spectral density Syzi^) similar to ([87|: 



2nSYzi^)Sin - n') = {tp\NY{n) o Nl{n')\ij) = (Nyin) o Nl{n')) , (90) 

one can get the following expressions for spectral densities in both cases of the vacuum state: 

y^{n)z{n) _ y:in)z,{n) + y:{n)Zs{n) 

" 2 " 2 ' ^^ ' 

and the squeezed state: 

s'^'zi^) = y^n){a{n) o a\n)),^,zin) = y\n)Ssq,{n)zin) 
= ^y\n)v[-4,p,^.,[2r,o]v[(i>]z{n) 

= yc(^)^c(») ^^^^^ 2r + sinh 2r cos 2</.) + y^*(^)^'^(^) (^osh 2r - sinh 2r cos 20) 
-&nM^ymM^,,^^,rsin2^, ' (92) 

Note that since the observables Y{t) and Z{t) that one calculates spectral densities for are 
Hermitian, it is compulsory, as is well known, for any operator to represent a physical quantity, 
then the following relation holds for their spectral coefficients 3^c,s(^) and Zc,s{^)'- 

y*c,si^) = ycA-^) , and Zl^{n) = Z,A-n) . (93) 

This leads to an interesting observation that the coefficients 3^c.s(^) and Zc^s{^) should be real- 
valued functions of variable s = iil. 

Now we can make further generalizations and consider multiple light and vacuum fields com- 
prising the quantity of interest: 

N 

Y{n) ^ iVy(f])-^3;J(f])d,(f7), (94) 

1=1 

where ai(f2) stand for quadrature amplitude vectors of N independent electromagnetic fields, 
and yl{il) are the corresponding complex- valued coefficient functions indicating how these fields 
are transmitted to the output. In reality, the readout observable of a GW detector is always a 
combination of the input light field and vacuum fields that mix into the output optical train as a 
result of optical loss of various origin. This statement can be exemplified by a single lossy mirror 



I/0-relations given by Eq. (34) of Section 2.2.4 



Thus, to calculate the spectral density for such an observable one needs to know the initial 
state of all light fields under consideration. Since we assume ai{fl) independent from each other, 
the initial state will simply be a direct product of the initial states for each of the fields: 

N 

I*) = (g) \i'^) , 
i=l 
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and the formula for the corresponding power (double-sided) spectral density reads: 

N N N 

Svm - ^3;l(r!)(i^,|d,(r!) o al{n)\^,)y,{n) = J2ylm^^my^m = E^^.(") ' (^5) 

with §i(ri) standing for the i-th input field spectral density matrix. Hence, the total spectral 
density is just a sum of spectral densities of each of the fields. The cross-spectral density for two 
observablcs Y{Q) and Z{i}) can be built by analogy and we leave this task to the reader. 
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4 Linear Quantum Measurement 

In Section[3J we discussed the quantum nature of light and fluctuations of the light field observables 
like phase and amplitude that stem thereof and yield what is usually called the quantum noise of 
optical measurement. In GW detection applications, where a sensitivity of the phase measurement 



is essential, as discussed in Section 2.1.3 the natural question arises: is there a limit to the 
measurement precision imposed by quantum mechanics? A seemingly simple answer would be 
that such a limit is set by the quantum fluctuations of the outgoing light phase quadrature, which 
are, in turn, governed by the quantum state the outgoing light finds itself in. The difficult part is 
that on its way through the interferometer, the light wave inflicts an additional back-action noise 
that adds up to the phase fluctuations of the incident wave and contaminates the output of the 
interferometer. The origin of this back action is in amplitude fluctuations of the incident light, 
giving rise to a random radiation pressure force that acts on the interferometer mirrors along with 
the signal GW force, thus effectively mimicking it. And it is the fundamental principle of quantum 
mechanics, the Heisenberg uncertainty principle, that sets a limit on the product of the phase 
and amplitude uncertainties (since these are complementary observables), thus leading up to the 
lower bound of the achievable precision of phase measurement. This limit appears to be a general 
feature for a very broad class of measurement known as linear measurement and is referred to as 
the SQL EI1[57|. 

In this section, we try to give a brief introduction to quantum measurement theory, starting 
from rather basic examples with discrete measurement and then passing to a general theory of 
continuous linear measurement. We introduce the concept of the SQL and derive it for special 
cases of probe bodies. We also discuss briefly possible ways to overcome this limit by contriving 
smarter ways of weak force measurement then direct coordinate monitoring. 

4.1 Quantum measurement of a classical force 

4.1.1 Discrete position measurement 

Let us consider a very simple measurement scheme, which, nevertheless, embodies all key features 
of a general position measurement. In the scheme shown in Figure [20} a sequence of very short light 
pulses are used to monitor the displacement of a probe body M . The position x of M is probed 
periodically with time interval d. In order to make our model more realistic, we suppose that 
each pulse reflects from the test mass F > 1 times, thus increasing the optomechanical coupling 
and thereby the information of the measured quantity contained in each reflected pulse. We also 
assume mass M large enough to neglect the displacement inflicted by the pulses radiation pressure 
in the course of the measurement process. 

Then each j'-th pulse, when reflected, carries a phase shift proportional to the value of the 
test-mass position x{tj) at the moment of reflection: 

4>J'' = 4>j - 2Fkpx{tj) , (96) 

where kp = Wp/c, cOp is the light frequency, j = ..., —1,0,1, .. . is the pulse number and (/)j is the 
initial (random) phase of the j-th pulse. We assume that the mean value of all these phases is 
equal to zero, {(f>j) — 0, and their root mean square (RMS) uncertainty ((0^) — {4>)'^Y^'^ is equal 
to A(/). 

The reflected pulses are detected by a phase-sensitive device (the phase detector). The imple- 
mentation of an optical phase detector is considered in detail in Section |2.3.1| Here we suppose 
only that the phase cj/^^ measurement error introduced by the detector is much smaller than the 
initial uncertainty of the phases A</). In this case, the initial uncertainty will be the only source of 
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Figure 20: Toy example of a linear optical position measurement. 



the position measurement error: 



Ax„ 



A(j) 



For convenience, we renormalize Eq. (|96|) as the equivalent test-mass displacement: 

X{tj) + Xfi{tj) , 



^r 



where 



Zr /Cp 



X{l{tj) 



Zr rh'Q 



are the independent random values with the RMS uncertainties given by Eq. (97) 



(97) 

(98) 
(99) 



Upon reflection, each light pulse kicks the test mass, transferring to it a back-action momentum 
equal to 

^'^"' (100) 



^after ^before ^b 



vr - v^ ' 



where p''°^°'''^ and fff^^'' are the test-mass momentum values just before and just after the light 
pulse reflection, and Wj is the energy of the j-th pulse. The major part of this perturbation is 
contributed by classical radiation pressure: 



^) 



b.a.' 



2r 



w. 



(101) 



with W the mean energy of the pulses. Therefore, one could neglect its effect, for it could be either 
subtracted from the measurement result or compensated by an actuator. The random part, which 
cannot be compensated, is proportional to the deviation of the pulse energy: 



(W,-W) 



(102) 



and its RMS uncertainly is equal to 



Apb.a. = 



2rAW 



(103) 
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with AW the RMS uncertainty of the pulse energy. 

The energy and the phase of each pulse are canonically conjugate observables and thus obey 
the following uncertainty relation: 

AWA0 > ^ . (104) 



Therefore, it follows from Eqs. (97 and 103) that the position measurement error Axmcas and the 



momentum perturbation Apb.a. due to back action also satisfy the uncertainty relation: 

AXmcasApb.a. > X • (1^5) 

This example represents a simple particular case of a linear measurement. This class of mea- 



surement schemes can be fully described by two linear equations of the form (98) and ( |100[ ), 
provided that both the measurement uncertainty and the object back-action perturbation (xnltj) 
and p^-^-{tj) in this case) are statistically independent of the test object initial quantum state and 
satisfy the same uncertainty relation as the measured observable and its canonically conjugate 
counterpart (the object position and momentum in this case). 

4.1.2 From discrete to continuous measurement 

Suppose the test mass to be heavy enough for a single pulse to either perturb its momentum notice- 
ably, or measure its position with the required precision (which is a perfectly realistic assumption 
for the kilogram-scale test masses of GW interferometers). In this case, many pulses should be 
used to accumulate the measurement precisions; at the same time, the test-mass momentum per- 
turbation will be accumulated as well. Choose now such a time interval T, which, on the one hand, 
is long enough to comprise a large number of individual pulses: 

iV = ^ » 1 , (106) 

and, on the other hand, is sufBciently short for the test-mass position x not to change considerably 
during this time due to the test-mass self-evolution. Then one can use all the N measurement 
results to refine the precision of the test-mass position x estimate, thus getting \/]V times smaller 
uncertainty 

AxT = ' !^'''' = Ax n^cas y Tf; ■ (107) 



At the same time, the accumulated random kicks the object received from each of the pulses 



random kicks, see Eq. (102), result in random change of the object's momentum similar to that of 



Brownian motion, and thus increasing in the same diffusive manner: 

Apt = Apb.a.\/^ = Apb.a.y ^ . (108) 

If we now assume the interval between the measurements to be infinitesimally small (i? — >■ 0) , 
keeping at the same time each single measurement strength infinitesimally weak: 

Aa;mcas ^ OO <^ Apb.a. -^ , 

then we get a continuous measurement of the test-mass position x{t) as a result. We need more 
adequate parameters to characterize its 'strength' than Axmcas and Apb.a.- For continuous mea- 
surement we introduce the following parameters instead: 



,2, ^0 ,, ,. (Apb.a.)^ 4r2& 



S, = lim {Axn^e^sY^ = T^IW > Sf = lim ^ ,' = ^^ , (109) 
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with 



S^ = lim {Acbfd , Si - lim ^^^ 



(110) 



This allows us to rewrite Eqs. ( 107[ ) and ( 108 1 in a form that does not contain the time interval d: 

AXT = ^, ApT^./S?f. (Ill) 



To clarify the physical meaning of the quantities Sx and 5'^ let us rewrite Eq. ( 98 1 in the continuous 
limit: 

x{t) ^ x{t) + Xii{t) , (112) 



where 



x^{t) 



m 

Z/- Kr, 



(113) 



stands for measurement noise, proportional to the phase 0(i) of the light beam (in the continuous 
limit the sequence of individual pulses transforms into a continuous beam). Then there is no 
difficulty in seeing that Sx is a power (double-sided) spectral density of this noise, and Srj, is a 
power double-sided spectral density of 4>{t). 

If we turn to Eq. (100), which describes the meter back action, and rewrite it in a continuous 



limit we will get the following differential equation for the object momentum: 

dp{t) 



dt 



Ffiit) 



(114) 



where ffl(i) is a continuous Markovian random force, defined as a limiting case of the following 
discrete Markov process: 



Ffiit j) = lim 



. p^-^-itj) 2r ,. Wj-W 2F 
— lim — 



C i9->-0 1? 



[i(t,)-io] 



(115) 



with T{t) the optical power, lo its mean value, and '. . . ' here meaning all forces (if any), acting 
on the object but having nothing to do with the meter (light, in our case). Double-sided power 
spectral density of Fb.a. is equal to Sp, and double-sided power spectral density of I is Si. 

We have just built a simple model of a continuous linear measurement, which nevertheless 
comprises the main features of a more general theory, i.e., it contains equations for the calculation 



of measurement noise (112) and also for back action (114). The precision of this measurement and 



the object back action in this case are described by the spectral densities Sx and Sp of the two 
meter noise sources, which are assumed to not be correlated in our simple model, and thus satisfy 



the following relation (cf. Eqs. (109l) 



St.Sf ~ 



Set, Si 



> 



(116) 



This relation (as well as its more general version to be discussed later) for continuous linear 



measurements plays the same role as the uncertainty relation ( 105 ) for discrete measurements 



establishing a universal connection between the accuracy of the monitoring and the perturbation 
of the monitored object. 
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Simple case: light in a coherent state. Recall now that scheme of representing the quantized 
light wave as a sequence of short statistically-independent pulses with duration e = -d we referred 



to in Section 3.2 It is the very concept we used here, and thus we can use it to calculate the 



spectral den sities of the measurement and back-action noise sources for our simple device featured 
in Figure 20 assuming the light to be in a coherent state with classical amplitude Ac — -\/2Xo/(^p) 



(we chose Ag = thus making the mean phase of light (cj)) =0). To do so we need to express phase 
(j) and energy W in the pulse in terms of the quadrature amplitudes cic.sit). This can be done if 
we refer to Eq. ( |61[ ) and make use of the following definition of the mean electromagnetic energy 
of the light wave contained in the volume v^ = Ac-d (here, A is the effective cross-sectional area of 
the light beam): 

dTE^{T) = W + SyV, (117) 

where W = v^CiAl/{8n) 



v^- 



v^ 



47r ^ ' 47r^ 



Iq^ is the mean pulse energy, and 



SW 



AdZl 

47r 



2Ar 



■djl 



-■0/2 



dracir) = ^2hL0pIa'dXij{t) = ^/2hujpWX.f,{t) 



(118) 



is a fluctuating part of the pulse energjPl We used here the definition of the mean pulse quadrature 
amplitude operators introduced in Eqs. (67). In the same manner, one can define a phase for each 



pulse using Eqs. (14) and with the assumption of small phase fluctuations (A(/) <C 1) one can get: 

1 /-^/^ r 



)?/2 



Thus, since in a coherent state A^J — AF^ 
equal to 



:Y,,^ 



huJri 



2W 



(119) 



A0= - 
^ 2 



1 hu,. 



W 



1/2 the phase and energy uncertainties are 
AW = y/hiUpW , (120) 



and hence 



Ax„ 




2F , 

Apb.a. = — \^huJp^^l . 

c 



(121) 



Substituting these expressions into Eqs. (110 1091, we get the following expressions for the power 



(double-sided) spectral densities of the measurement and back-action noise sources: 



Sdi 



^ 

4Xo 



hujpTQ , 



and 



0_7. 



l&ujplQF"^ 



Sp = 



Ahuipi^r"^ 



(122) 



(123) 



We should emphasize that this simple measurement model and the corresponding uncertainty 



relation (116) are by no means general. We have made several rather strong assumptions in the 



course of derivation, i.e., we assumed: 

1. energy and phase fluctuations in each of the light pulses uncorrelated: {yV{tj)(j){tj)) — 0; 

2. all pulses to have the same energy and phase uncertainties AW and A(/), respectively; 



*Here, we omitted the terms of <5VV proportional to a^ ^(i) since their contribution to the integral is of the second 
order of smallness in dc,s/-Ac compared to the one for the first order term. 
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3. the pulses statistically independent from each other, particularly taking {W{ti)W{tj)) = 
{4>{U)4>{tj)) = {W{U)^{tj)) = with U 7^ tj. 

These assumptions can be mapped to the following features of the fluctuations xn (i) and i^b.a. (t) 
in the continuous case: 

1. these noise sources are mutually not correlated; 

2. they are stationary (invariant to the time shift) and, therefore, can be described by spectral 
densities S^ and Sp; 

3. they are Markovian (white) with constant (frequency-independent) spectral densities. 

The features 1 and 2, in turn, lead to characteristic fundamentally-looking sensitivity limita- 
tions, the SQL. We will call linear quantum meters, which obey these limitations (that is, with 
mutually non-correlated and stationary noises xa and -Fb.a.), Simple Quantum Meters (SQM). 

4.2 General linear measurement 

In this Section, we generalize the concept of linear quantum measurement discussed above and give 
a comprehensive overview of the formalism introduced in '27] and further elaborated in '38) I48j . 
This formalism can be applied to any system that performs a transformation from some unknown 
classical observable (e.g., GW tidal force in GW interferometers) into another classical observable 
of a measurement device that can be measured with (ideally) arbitrarily high precision (e.g., in 
GW detectors, the readout photocurrent serves such an observable) and its value depends on the 
value of unknown observable linearly. For definiteness, let us keep closer to GW detectors and 
assume the continuous measurement of a classical force. 



Meter 



T 



Probe 



G 



-*x 



Figure 21: General scheme of the continuous linear measurement with G standing for measured 
classical force, X the measurement noise, JF the back-action noise, O the meter readout observable, 
X the actual probe's displacement. 



The abstract scheme of such a device is drawn in Figure 21 It consists of a probe V that is 
exposed to the action of a classical force G(t), and the meter. The action of this force on the 
probe causes its displacement x that is monitored by the meter (e.g., light, circulating in the 
interferometer). The output observable of the meter O is monitored by some arbitrary classical 
device that makes a measurement record o{t). The quantum nature of the probe-meter interaction 
is reflected by the back-action force F that randomly kicks the probe on the part of the meter 
(e.g., radiation pressure fluctuations). At the same time, the meter itself is the source of additional 
quantum noise Ofl(t) in the readout signal. Quantum mechanically, this system can be described 
by the following Hamiltonian: 

^=^!,"obe+C°L- + l>W, (124) 

where 'Hp^.^i-.j, and Hooter &re the Hamiltonians describing the free evolution of the probe and the 
meter, respectively, i.e., when there is no coupling between these systems, and V{t) = —x{G{t) + F) 
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is the interaction Hamiltonian. The evolution of this system can be found, in general, by solving 
Heisenberg equations for all of the system observables. However, it is convenient to rewrite it first 
in the interaction picture factoring out the free evolution of the probe and the meter (see Appendix 
3.7 of [48j for detailed derivation): 



H(°)(i - to) U -x^"\t){G{t) + F(")(t)) , (125) 



ni{t) = exp \ ^n^°\t - to) )■ t/(t)exp 

where 'H'-^' — 'Hpi.obo + 'Hrnctcr; ^iid x^-^\t) and F'^^\t) are the Heisenberg operators of the probe's 
displacement and the meter back-action force, respectively, in the case of no coupling between 
these systems, i.e., the solution to the following system of independent Heisenberg equations: 



dx(o)(f) 
Jt 



n 



(0) ^(0) 
probe' 



it) 



dF(o)(t) 
dt 



n 



(0) 



,F<^'\t) 



and to is the arbitrary initial moment of time that can be set to — oo without loss of generality. 

The following statement can be proven (see |98J, Section VI of 27J, and Theorems 3 and 4 
in Appendix 3.7 of [48J for proof): For a linear system with Hamiltonian (124 1, for any linear 



observable A of the probe and for any linear observable B of the meter, their full Heisenberg 
evolutions are given by: 



A{t) = m{t)+ dt'xA.{t,t')[F{t') + G{t')]. 



to 



m 



m 



(t) + / dt' 

Jtn 



XBF{t,t')x{t') , 



(126) 



where A^^'[t) and B^'~''{t) stand for the free Heisenberg evolutions in the case of no coupling, and 
the functions XAxit^t') and XBF(t,t') are called (time-domain) susceptibilities and defined as: 



XAx{t,t') 



XBF{t,t') 



i(o)(t),f(o)(t') 

t < t' 

SW(i),F(")(t')' 
t < t' 



t ^ t' 



t ^ t' 



(127) 



The second clauses in these equations maintain the causality principle. 

For time independent Hamiltonian W^' and operator F (in the Schrodinger picture), the sus- 
ceptibilities are invariant to time shifts, i.e., x(f, t') = x(t + T,t' -\- r), therefore they depend only 
on the difference of times: x(t, t') 
domain as: 



x{t — t'). In this case, one can rewrite Eqs. ( |126| in frequency 

(128) 



A{n) = m{n) -f xA.{mFm + G{n)] , B{^) - misi) + XBFmm) , 

where the Fourier transforms of all of the observables are defined in accordance with Eq. ^ . 

Let us now use these theorems to find the full set of equations of motion for the system of 
linear observables x, F and O that fully characterize our linear measurement process in the scheme 
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featured in Figure [2T| 



0(i) = 6(")(i)+ / dt' xoF{t ~ t')x(t') , 
Fit) - F(°)(t)+ / dt'xFF{t-t')xit'), 

Jto 

x{t) = i(°)(i) + / dt' xUt - t') \G{t') + Fit') 



(129) 



where time-domain susceptibilities are defined as 



XoFit-t') 

XFFit-t') 



6(°)(i),F(")(i') 
F(")(t),F(o)(t') 



(130) 



The meaning of the above equations is worth discussing. The first of Eqs. ( 129 ) describes how 



the readout observable 0(t) of the meter, say the particular quadrature of the outgoing light field 
measured by the homodyne detector (cf. Eq. (39)), depends on the actual displacement xit) of 
the probe, and the corresponding susceptibility XOFit — t') is the transfer function for the meter 
from X to O. The term O^'^^t) stands for the free evolution of the readout observable, provided 
that there was no coupling between the probe and the meter. In the case of the GW detector, 
this is just a pure quantum noise of the outgoing light that would have come out were all of the 
interferometer test masses fixed. It was shown explicitly in P¥j and we will demonstrate below 
that this noise is fully equivalent to that of the input light except for the insignificant phase shift 
acquired by the light in the course of propagation through the interferometer. 

The following important remark should be made concerning the meter's output observable 
0(i). As we have mentioned already, the output observable in the linear measurement process 
should be precisely measurable at any instance of time. This implies a simultaneous measurability 
condition [351 SHI 11481 \T7[ H51 155] on the observable 0(i) requiring that it should commute with 
itself at any moment of time: 



0(i), Oit') 



= , Vt, i' . 



(131) 



Initially, this condition was introduced as the definition of the quantum non-demolition (QND) 
observables by Braginsky et al. [32, 33 . In our case it means that the measurement of 0(ti) at 
some moment of time ti shall not disturb the measurement result at any other moments of time 

and therefore the sample data < 0(ti), 0(ii ),..., 0(t„) ^ can be stored directly as bits of classical 
data in a classical storage medium, and any noise from subsequent processing of the signal can be 
made arbitrarily small. It means that all noise sources of quantum origin are already included in 
the quantum fluctuations of 0(i) [351 |3S]. And the fact that due to (131) this susceptibihty turns 



out to be zero reflects the fact that 0(i) should be a classical observable. 

The second equation in ( |129 ) describes how the back-action force exerted by the meter on the 
probe system evolves in time and how it depends on the probe's displacement. The flrst term, 
F^^\t), meaning is rather obvious. In GW interferometer, it is the radiation pressure force that 
the light exerts on the mirrors while reflecting off them. It depends only on the mean value and 
quantum fluctuations of the amplitude of the incident light and does not depend on the mirror 
motion. The second term here stands for a dynamical back-action of the meter and since, by 
construction, it is the part of the back-action force that depends, in a linear way, from the probe's 
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displacement, the meaning of the susceptibiHty XFF{t — t') becomes apparent: it is the generaUzed 
rigidity that the meter introduces, effectively modifying the dynamics of the probe. We will see 
later how this effective rigidity can be used to improve the sensitivity of the GW interferometers 
without introducing additional noise and thus enhancing the SNR of the GW detection process. 



The third equation of ( 129 ) concerns the evolution of the probe's displacement in time. Three 



distinct parts comprise this evolution. Let us start with the second and the third ones: 

Xs{t)= j dt'xxx{t-t')G{t'), and x^.^X*) ^ f dt' Xx.{t ' t')F{t') . (132) 

Jto Jto 

Here Xs{t) is the probe's response on the signal force G{t) and is, actually, the part we are mostly 
interested in. This expression also unravels the role of susceptibility Xxx{i — i')'- it is just the 
Green's function of the equation of motion governing the probe's bare dynamics (also known as 
impulse response function) that can be shown to be a solution of the following initial value problem: 

T)xxx{t - t') = 5(t - t') , Xxx{t-t')\t^t'=Q, ... 
5"-2xxx(i-i')| _Q 5"-^X.x(t-t')| _ 1 



where D = X]fe=o '^*=^ ^^ ^^^ linear differential operator that is governed by the dynamics of the 
probe, e.g., it is equal to Df.^. — M^ for a free mass M and to Dqsc = ^^^ + ^'-^^m for a 
harmonic oscillator with eigenfrequency ^m ■ Apparently, operator D is an inverse of the integral 
operator Xxx whose kernel is Xxx{t — t'): 

Xs(t)^-D-^G{t)^XxxG{t)^ f dt'xxx{t~t')G{t'). 



The second value, Xb.a. [i), is the displacement of the probe due to the back-action force exerted 
by the meter on the probe. Since it enters the probe's response in the very same way the signal 
does, it is the most problematic part of the quantum noise that, as we demonstrate later, imposes 
the SQL [211 Eg. 

And finally, x'"' (t) simply features a free evolution of the probe in accordance with its equations 
of motion and thus depends on the initial values of the probe's displacement x'^^^Iq), momentum 
p^°^(io), and, possibly, on higher order time derivatives of x^^\t) taken at to, as per the structure 
of the operator D governing the probe's dynamics. It is this part of the actual displacement 
that bears quantum uncertainties imposed by the initial quantum state of the probe. One could 
argue that these uncertainties might become a source of additional quantum noise obstructing the 
detection of GWs, augmenting the noise of the meter. This is not the case as was shown explicitly 
in [23], since our primary interest is in the detection of a classical force rather than the probe's 
displacement. Therefore, performing over the measured data record o{t) the linear transformation 
corresponding to first applying the operator Xof '^^ ^^^ readout quantity that results in expressing 
o(i) in terms of the probe's displacement: 

i(i) = Xo'^pOit) = Xfi{t) + a;("'(t) + Xs{t) + xy,,,Xt) , 

with Xfl(t) — XoF^ i^) standing for the meter's own quantum noise (measurement uncertainty), 
and then applying a probe dynamics operator D that yields a force signal equivalent to the readout 
quantity o{t): 

Fit) = D£(i) = Bxait) + Fb.a.(i) + G(i) . 

The term 'Dx'^^\t) vanishes since x'^°^(t) is the solution of a free-evolution equation of motion. 
Thus, we see that the result of measurement contains two noise sources, Xfl(t) and Fb.a.(i), which 
comprise the sum noise masking the signal force G{t). 
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Since we can remove initial quantum uncertainties associated with the state of the probe, it 



would be beneficial to turn to the Fourier domain and rewrite Eqs. ( 129 ) in the spectral form 

d{n) = d(°'>{n) + xoFmm), 
F{n) = FW(r!) + XFF(^)i(^), 



where spectral susceptibilities are defined as: 

XABm^ / dTXAB{r)e'''\ 
Jo 



(133) 



(134) 



with {A, B) => (O, F, x), and we omit the term x'"' (il) in the last equation for the reasons discussed 
above. The solution of these equations is straightforward to get and reads: 



6{n) ^ 

F{n) = 
x{n) = 



d^"\n) 



Xxx{^)XOFi^) 



1 - Xxx{^)XFFi^) 



1 



1 - Xxxi^)XFF{^) 
Xxxi^) 

^~Xxx{^)xff{^) 



G{n) + P'-'^^n) 
F(°'>{n) + XFFmxxxmG{n) 
G{n) + F'^"\n) 



(135) 
(136) 
(137) 



It is common to normalize the output quantity of the meter 0{n,) to unit signal. In GW 
interferometers, two such normalizations are popular. The first one tends to consider the tidal 



force G as a signal and thus set to 1 the coefficient in front of G{Q) in Eq. (1351. The other one 



takes GW spectral amplitude h{Q) as a signal and sets the corresponding coefficient in 0{CI) to 
unity. Basically, these normalizations are equivalent by virtue of Eq. (|12p as: 



Mn^xhin) = -MQ^^^^ = Gin) 



h{n) = -2G{n)/{MLV? 



(138) 



In both cases, the renormalized output quantities can be considered as a sum of the noise and 
signal constituents: 

C'^=7V'-^ + G or 6'' = Af" + h{n) . (139) 

And it is the noise term in both cases that we are seeking to calculate to determine the sensitivity 
of the GW detector. Let us rewrite 0{n) in force normalization: 



o^(r!) 



_ 1 -XFF{n)Xxx[n) 
XOF{n)Xxx{n) 

x{n) 



o{n) 



O(o)(r!) 



XOFin)Xxxin) 



+ {F^'>in) 



^^d(")(17))+G(f^) 



j'in) + G{n) 



Xxx{n) 

where we introduce two new linear observables X and J^ of the meter defined as: 



XOFin) 

that have the following meaning: 



^{n) 



XOFi^) 



'm, 



(140) 



(141) 



• A" is the effective output fluctuation of the meter not dependent on the probe. Henceforth, 
we will refer to it as the effective measurement noise (shot noise, in the GW interferometer 
common terminology); 
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• J^ is the effective response of the output at time t to the meter's back-action force at earher 
times t < t' . In the foUowing we wiU refer to J- as the effective hack-action noise (radiation- 
pressure noise, in the GW interferometer common terminology). 

These two new observables that embody the two types of noise inherent in any Unear measurement 
satisfy the following commutation relations: 






Pm,Pm 



0, 



-2mh5{Q. - ^') , 



X{t),X{t') ^ P{t),P{t') =0,(142) 
X{t),P\t')\^-ih6{t~t'), (143) 



that can be interpreted in the way that X{t) and J^{t) can be seen at each instance of time 
as the canonical momentum and the coordinate of different effective linear measuring devices 
(meter -f probe) , thus defining an infinite set of subsequent measurements similar to the successive 
independent monitors of von Neumann's model [157]. In this case, however, the monitors described 
by X{t) and J^{t) are not, generally speaking, independent. In GW detectors, these monitors appear 
correlated due to the internal dynamics of the detector, i.e., the noise processes they describe are 
non-Markovian. 

In particular, this can be seen when one calculates the power (double-sided) spectral density 
of the sum noise N^ (t) : 



s^{n) 



dt {M^{t) o ^^(i'))e'"(*-*') = p^^^ + SM^) + 2Re 






where spectral densities: 



Sxxm = 



dt{X{t)oX{t'))e'^^''''\ 
dt{f{t)oT{t'))e'^'^'-*~''\ 
dt{X{t)oT{t'))e'^^*-''\ 



(144) 

(145) 
(146) 
(147) 



are not necessarily constant and, thus, describe non-Markovian random processes. It can also be 



shown that since X{i) and P{t) satisfy commutation relations (142), their spectral densities shall 
satisfy the Schrodinger-Rohertson uncertainty relation: 



Sxx{n)S^A^)~\SxA^)?^ 



(148) 



that is the generalization of a Heisenberg uncertainty relation in the case of correlated observables. 
The general structure of quantum noise in the linear measurement process, comprising two 



types of noise sources whose spectral densities are bound by the uncertainty relation (148), gives 



a clue to several rather important corollaries. One of the most important is the emergence of the 
SQL, which we consider in detail below. 



4.3 Standard Quantum Limit 



RecaU the SQM in Section [4X2 



The SQM has non-correlated effective measurement and back-action noises that results in 
Sxt{^) = 0. Apparently, under these conditions X and J- turn into ig and Ff^ of Eqs. (112) and 
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(114), respectively. Hence, we will use S^i^) instead of Sxx{^) and Sp{il) instead of Sjrjr[^) 



Then the uncertainty relation ( 148 1 transforms into: 



s^in)SFin) > ^ . (149) 

The SQL is the name for an ultimate lower bound of a sum noise spectral density the SQM 
can, in principle, have at any given frequency fl. To derive this limit we assume noise sources xn 
and i^b.a. to have minimal values allowed by quantum mechanics, i.e. 

s^{n)Sp{n) = ^ . (150) 

Then, using this condition, one can minimize SQM's sum noise: 
to yield: 

'«^'"' = EM <'^'' 

, that is achieved when contributions of measurement noise and back-action noise to the sum noise 
are equal to each other, i.e., when 

s^A^) = ^\x..m\, ^^ Spin) = ^|^^^^^)| ■ (152) 

It is instructive to cite the forms of the SQL in other normalizations, i.e., for /i-normalization 
and for ^-normalization. The former is obtained from (1511 via multiplication by 4/(M^L^i7^): 



^ 4S'fQL(^) 4fi 

^sql(") = M^L^n^ = APL^n^xUm ■ ^^^^^ 



The latter is obtained fromEq. ( 151[ ) using the obvious connection between force and displacement 

x{n) ^ xxxi^)F{n): 

siQLm - ixxx(^)P5|ql(") = h\xum ■ (154) 

These limits look fundamental. There are no parameters of the meter (only ft as a reminder 
of the uncertainty relation ( |116| )), and only the probe's dynamics is in there. Nevertheless, this is 
not the case and, actually, this limit can be beaten by more sophisticated, but still linear, position 
meters. 

At the same time, the SQL represents an important landmark beyond which the ordinary brute- 
force methods of sensitivity improving cease working, and methods that allow one to blot out the 
back-action noise J^{t) from the meter output signal have to be used instead. Due to this reason, 
the SQL, and especially the SQL for the simplest test object - free mass - is usually considered as 
a borderline between the classical and the quantum domains. 

4.3.1 Free mass SQL 

In the rest of this section, we consider in more detail the SQLs for a free mass and for a harmonic 



oscillator. We also assume the minimal quantum noise requirement (150) to hold. 

The free mass is not only the simplest model for the probe's dynamics, but also the most 
important class of test objects for GW detection. Test masses of GW detectors must be isolated 
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as much as possible from the noisy environm.ent. To this end, the design of GW interferometers 
imphes suspension of the test masses on thin fibers. The real suspensions are rather sophisticated 
and comprise several stages slung one over another, with mechanical eigenfrequencies fm in ^ IHz 
range. The sufficient degree of isolation is provided at frequencies much higher than fm, where the 
dynamics of test masses can be approximated with good precision by that of a free mass. 

Let us introduce the following convenient measure of measurement strength (precision) in the 

first place: 

/ o \ 1/4 



Using the uncertainty relation (1501, the noise spectral densities Sx and Sp can be expressed 
through rtq as follows: 

Therefore, the larger flq is, the smaller Sx is (the higher is the measurement precision), and the 
larger Sp is (the stronger the meter back action is). 

In the case of interferometers, 17^ is proportional to the circulating optical power. For example, 
for the toy optical meter considered above. 






^. = \/HV^> (157) 



see Eqs. (123). Using this notation, and taking into account that for a free mass Af, 



xL"(^) = -^, (158) 

the sum quantum noise power (double-sided) spectral density can be written as follows: 



The SQL optimization (152 1 takes the following simple form in this case: 

nq^n, (160) 

giving: 



This consideration is illustrated by Figure 22 (left), where power (double-sided) spectral den- 
sity (159) is plotted for three different values of Vtq. It is easy to see that these plots never dive 



under the SQL line (161), which embodies a common envelope for them. Due to this reason, the 
sensitivities area above this line is typically considered as the 'classical domain', and below it - as 
the 'quantum domain'. 

4.3.2 Harmonic oscillator SQL 



The simplest way to overcome the limit ( 161 ), which does not require any quantum tricks with the 
meter, is to use a harmonic oscillator as a test object, instead of the free mass. It is easy to see 
from Eq.( |151[ ) that the more responsive the test object is at some given frequency fi (that is, the 
bigger Xxx{^)) is, the smaller its force SQL at this frequency is. In the harmonic oscillator case. 
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Figure 22: Sum quantum noise power (double-sided) spectral densities of the Simple Quantum 
Meter for different values of measurement strength fig (red) < Sl^ (green) < J7g(blue). Thin black 
line: SQL. Left: free mass. Right: harmonic oscillator 



with J7o standing for the oscillator mechanical eigenfrequency, and the sum quantum noise power 
(double-sided) spectral density equal to 



s',jn)^M'int-n'rs, + SF 



hMfil 



2^2 






1 



(163) 



Due to a strong response of the harmonic oscillator near resonance, the first (measurement noise) 



term in Eq. (159) goes to zero in the vicinity of flo. Therefore, reducing the value of flq, that is, 
using weaker measurement, it is possible to increase the sensitivity in a narrow band around VIq- 
At the same time, the smaller f2g is, the more narrow the bandwidth is where this sensitivity is 



achieved, as can be seen from the plots drawn in Figure 22 (right). 

Consider, in particular, the following minimax optimization of the narrow-band sensitivity. Let 
J/ = fi — Oo be the detuning from the resonance frequency. Suppose also 



W\ <^o. 



(164) 



In this case, the sum noise power (double-sided) spectral density ( 163 1 can be approximated as 
follows: 



'S'oscC^O 



V K 



1 



(165) 



Then require the maximum of S'^^ in a given frequency range Afi be as small as possible. 

It is evident that this frequency range has to be centered around the resonance frequency iloj 
with the maximums at its edges, v = ±Arj/2. The sum noise power (double-sided) spectral density 
at these points is equal to 



hM^l /nll\n'' 



Cc(f^o±Afl/2) 
The minimum of this expression is provided by 



^^ 



1 



(166) 



(167) 
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Substitution of this value back into Eq. ( 165 ) gives the fohowing optimized power (double-sided) 
spectral density: 

Si(". + ^)^'^('^ + An1, (168, 



with 



s^^^i^o ± Ar2/2) == hunoAn . (i69) 

Therefore, the harmonic oscillator can provide a narrow-band sensitivity gain, compared to the 



free mass SQL (1611, which reads 



Sosc 



OF 
"-"SQLf.m. 



(^o) 



1 /4iy2 






and can be further written accounting for the above optimization as: 



^o''sc(^0+^) 



•^SQLf.m.V^o) |i/|<An/2 

Of course, the oscillator SQL, equal to 



< 



S^,,{no ± A»/2) ^ A» 

•5'sQLf.m.(^o) ^0 



^sQLosc = fiM\nl - n'\ « 2hMnoW\ 



(170) 



(171) 



(172) 



cannot be beaten is this way, and the question of whether the sensitivity ( 171 1 is the 'true' beating 



of the SQL or not, is the question to answer (and the subject of many discussions). 



4.3.3 Sensitivity in different normalizations. Free mass and harmonic oscillator 

Above, we have discussed, in brief, different normalizations of the sum noise spectral density and 



derived the general expressions for the SQL in these normalizations (cf. Eqs. (153) and (154)). 
Let us consider how these expressions look for the free mass and harmonic oscillator and how the 
sensitivity curves transform when changing to different normalizations. 



/i-normalization: The noise spectral density in /i-normalization can be obtained using Eq. ( 12 1 
Where the SQM is concerned, the sum noise in /i-normalization reads 



,{n) = u''{a) 



%(») 



Fain) 



ML9? 
In the case of a free mass with Xxx{^) = —l/{Mfl^) the above expression transforms as: 

and that results in the following power (double-sided) spectral density formula: 



Str.M = 



L2 



Sr. 



Sp 




ML^ni 



and results in the following formula for free mass SQL in /i-normalization: 



■^SQLf.m.l^) — 



Ah 



(173) 



(174) 
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Figure 23: Sum quantum noise power (double-sided) spectral densities of the Simple Quantum 
Meter in the /i- normalization for different values of measurement strength: r2g(red) < fJq (green) < 
J7q(blue). Thin black line: SQL. Left: free mass. Right: harmonic oscillator. 



The plots of these spectral densities at different values of fig are given in the left panel of 
Figure [23| 

As for the harmonic oscillator, similar formulas can be obtained taking into account that 
X°f^i^) = ^/{M{nl - n^)). Thus, one has: 



, osc 
sum 



m 



2(r2§ - i72)£fl(r2) 



2^fl 



that results in the following power (double-sided) spectral density formula: 



'-'oscl") — 



i2 



1 



"0 



t^x 



Sp 



ivpn^ 



2h 



ML^ni 






2^2 






and results in the following formula for free mass SQL in /i-normalization: 



'SQLosc 



m 






(175) 



(176) 



The corresponding plots are drawn in the right panel of Figure 23 Despite a quite different 
look, in essence, these spectral densities are the same force spectral densities as those drawn in 
Figure 22 yet tilted rightwards by virtue of factor 1/J7^. In particular, they are characterized by 



the same minimum at the resonance frequency, created by the strong response of the harmonic 



oscillator on a near-resonance force, as the corresponding force-normalized spectral densities ( 163 



172). 



^-normalization: Another normalization that is worth considering is the actual probe displace- 
ment, or a;-normalization. In this normalization, the sum noise spectrum is obtained by multiplying 



noise term J\f^{n) in Eq. (139) by the probe's susceptibility 



isum(^) = ifl(f^) + Xxxi^)Ffi{Q) . 



(177) 
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It looks rather natural at a first glance; however, as we have shown below, it is less heuristic 
than the force normalization and could even be misleading. Nevertheless, for completeness, we 
consider this normalization here. 

Spectral density of Xsum(r2) and the corresponding SQL are equal to 






(178) 
(179) 



In the free mass case, the formulas are the same as in ^.-normalization except for the multipli- 
cation by 4/L^: 



.m 



with SQL equal to: 



\s + ^' 1 


h 


fi + 51 


2Nmi 


SSQLL^M 


h 

n rn2 ' 





In the harmonic oscillator case, these equations have the following form: 



5:,,(r!) = 5,+ 



AP{nl-n^f 



2Kmi 









s, 



SQL osc 



m 



Tn\ni~n^\ 



(180) 
(181) 

(182) 
(183) 



The corresponding plot of the harmonic oscillator power (double-sided) spectral density in 
a:-normalization is given in Figure [24| 



3 
Co 




Figure 24: Sum quantum noise power (double-sided) spectral densities of Simple Quantum Meter 
and harmonic oscillator in displacement normalization for different values of measurement strength: 
J7g(red) < rig(green) < riq(blue). Thin black line: SQL. 
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Note that the curves display a sharp upsurge of noise around the resonance frequencies. How- 
ever, the resonance growth of the displacement due to signal force G have a long start over this 
seeming noise outburst, as we have shown already, leads to the substantial sensitivity gain for a 
near-resonance force. This sensitivity increase is clearly visible in the force and equivalent dis- 
placement normalization, see Figures [22] and [23) but completely masked in Figure [24) 

4.4 Beating the SQL by means of noise cancellation 

The SQL is not a fundamental limitation as we have mentioned already, and the clue to how 
to overcome it can be devised from the expression for the general linear measurement sum noise 



spectral density (144). One can see that a properly constructed cross-correlation between the 
measurement noise X = O'^'^^ /xof{^) and back-action noise F'^"\ i.e., the right choice of xff{^) 
that should be at any frequency equal to: 

Xff{^) = SxFm/Sxxm (184) 

can compensate the back-action term and leave only the measurement noise-related contribution 
to the final sum quantum noise: 

This spectral density could be arbitrarily small, providing the unbound measurement strength. 
However, there is a significant obstacle on the way towards back-action free measurement: the 
optimal correlation should be frequency dependent in the right and rather peculiar way. Another 
drawback of such back-action evasion via noise cancellation resides in the dissipation that is always 
present in real measurement setups and, according to Fluctuation-Dissipation Theorem [?51 [S3] is 
a source of additional noise that undermines any quantum correlations that might be built in the 
ideal system. 

The simplest way is to make the relation ( 184[) hold at some fixed frequency, which can always be 



done either (i) by preparing the meter in some special initial quantum state that has measurement 
and back-action fluctuations correlated (Unruh |149l I148J proposed to prepare input light in a 
squeezed state to achieve such correlations), or (ii) by monitoring a linear combination of the 
probe's displacement and momentum |1621 11591 11581 11601 11611 [551 [55] that can be accomplished, 
e.g., via homodyne detection, as we demonstrate below. 

We consider the basic principles of the schemes, utilizing the noise cancellation via building 
cross-correlations between the measurement and back-action noise. We start from the very toy 
example discussed in Section [4.1.1| 

The advanced version of that example is shown in Figure |25| The only difference between this 



scheme and the initial one (see Figure 20 ) is that here the detector measures not the phase of light 



pulses, but linear combination of the phase and energy, parametrized by the homodyne angle ^lo 



(cf. Eq. ([39])): 

6{t,) - -C" sin </)Lo + "'" cos (/-Lo , (185) 



(we subtracted the regular term proportional to the mean energy W from the output signal 0{tj)) 



Similar to Ec^. (98), renormalize this output signal as the equivalent test object displacement: 



2rfcpSin(/iLo 
where the noise term has in this case the following form: 



iitj) = — - — ^- — = x{tj) + xn{tj) , (186) 



i^itj) = ^ {^^4>{tj) + ^%y^ cot^LO ) . (187) 
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Figure 25: Toy example of a linear optical position measurement. 



RMS uncertainty of this value (the measurement error) is equal to 

1 



AXn 



Zr rCr] 



(A0)- (^>^)^..^. 



4W2 



cot^ 



'LO ■ 



(188) 



At first glance, it seems like we just obtained the increased measurement error, for the same value 



of the test object perturbation, which is still described by Eq. (103 1. However, the additional term 



in Eq. ( 185 ) can be viewed not only as the additional noise, but as the source of information about 



the test object perturbation. It can be used to subtract, at least in part, the terms induced by 
this perturbation from the subsequent measurement results. Quantitatively, this information is 
characterized by the cross-correlation of the measurement error and the back action: 



A{xp) = {x4t,)op^---{tj)) 



jAwy 



cot ( 



^LO ■ 



(189) 



It is easy to see that the uncertainties (103 188 1891 satisfy the following Schrodinger-Robertson 
uncertainty relation: 



^2,A ^2 TA. ^12 _ (A0)2(AW)2 . h^ 



(AXmoas)nApb.a.)^ - [A{xp)Y ^ 



> 



(190) 



Now we can perform the transition to the continuous measurement limit as we did in Sec- 
tion gTl 



op 



lim(Aa;,„oas)^^ = TF?T? 
,. (Apb.a.)' ^r'Si 



Sd, 



SxF = lim A{xp) = - — - 



COt( 



'LO 



^COt 0LO 



(191) 



66 



which transform inequaUty ( 190 1 to the Schodinger-Robertson uncertainty relation for continuous 
measurements: 



i^x^F ^xF 



Sif,Si h 



(192) 



In the particular case of the light pulses in a coherent quantum state (1201, the measurement 



error (188), the momentum perturbation (1031, and the cross-correlation term (189) are equal to: 

Apb.a. = y/hoJpW , A(xp) = -COt</)LO (193) 



4r y ujpW sm^ (/)LO 

(the momentum perturbation Apb.a. evidently remains the same as in the uncorrelated case, and 
we provided its value here only for convenience) , which gives the exact equality in the Schrodinger- 
Robertson uncertainty relation: 



(Aa;n,eas)^(Apb.a.)^ " [A(xp)]2 = 



(194) 



Correspondingly, substituting the coherent quantum state power (double-sided) spectral densi- 



ties (122) into Eqs. (191), we obtain 



i'x 



WlUpIqF'^ sin^ ^lo 



with 



Sf 



^x^F ^xF 



, SxF = 2 '^°* "^LO 



(195) 



(196) 



[compare with Eqs. ( 123 )] 



The cross-correlation between the measurement and back-action fluctuations is equivalent to 
the virtual rigidity xPf(^) = ^-^vht = const, as one can conclude looking at Eqs. (141). Indeed, 

Xfi{t) = X{t) , and F(") = T{t) + K^MXs{t) . 

The sum noise does not change under the above transformation and can be written as: 

Fsum(t) = Df fi(i) + F^^\t) = Deffifl(t) + P{t) , 

where the new effective dynamics that correspond to the new noise are governed by the following 
differential operator 

Deff = D + ifvirt . (197) 

The above explains why we refer to -ffvirt as 'virtual rigidity'. 

To see how virtual rigidity created by cross-correlation of noise sources can help beat the SQL 
consider a free mass as a probe body in the above considered toy example. The modified dynamics: 



Deff =M-^-hii:virt, 



SxF — K^irtSx 7^ , 



(198) 



correspond to a harmonic oscillator with eigenfrequency Hq = K^irt/M, and as we have demon- 



strated in Eq. ( 171 ) provide a narrow-band sensitivity gain versus a free mass SQL near the 



resonance frequency f2o • 

However, there is a drawback of virtual rigidity compared to the real one: it requires higher 
measurement strength, and therefore higher power, to reach the same gain in sensitivity as provided 
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by a harmonic oscillator. This becomes evident if one weighs the back-action spectral density Sp, 



which is a good measure of measurement strength according to Eqs. (1561, for the virtual rigidity 



against the real one. For the latter, to overcome the free mass SQL by a factor 



e 



Art 



(199) 



(see Eq. (171)) at a given frequency flo, the back-action noise spectral density has to be reduced 
by this factor: 

^2^optf.m. ^ (200) 



hMnl 



see Eqs. (156 160 167). Here, 5^^' ■™' ~ hMD,o/2 is the back-action noise spectral density, which 



allows one to reach the free mass SQL (161) at frequency JIq- Such a sensitivity gain is achieved 



at the expense of proportionally reduced bandwidth: 

An = frio ■ 



For the virtual rigidity, the optimal value of Sp results from Eq. (167) 



Sp^ 



h^/A + si 

Sx 



hM 



ni 



"0 



Qopt f .ni. 

— Op 



e 



1 



(201) 



(202) 



Hence, the better the sensitivity (the smaller ^^), the larger Sp must be and, therefore, measure- 
ment strength. 

Another evident flaw of the virtual rigidity, which it shares with the real one, is the narrow-band 
character of the sensitivity gain it provides around Oq and that this band shrinks as the sensitivity 
gain rises (cf. Eq. (199)). In order to provide a broadband enhancement in sensitivity, either the 



real rigidity K = A/fip, or the virtual one K^irt = Sxp/Sp should depend on frequency in such 
a way as to be proportional (if only approximately) to 51^ in the frequency band of interest. Of 
all the proposed solutions providing frequency dependent virtual rigidity, the most well known are 
the quantum speedmeter [26j and the filter cavities |94| schemes. Section 4.5 we consider the basic 



principles of the former scheme. Then, in Section [6] we provide a detailed treatment of both of 
them. 



4.5 Quantum speedmeter 

4.5.1 The idea of the quantum speedmeter 

The toy scheme that demonstrates a bare idea of the quantum speedmeter is shown in Figure |26[ 



The main difference of this scheme from the position meters considered above (see Figures 20 25 ) 
is that each light pulse reflects from the test mass twice: first from the front and then from the 
rear face after passing the delay line with delay time r. An outgoing pulse acquires a phase shift 
proportional to the difference of the test-object positions at time moments separated by r, which 
is proportional to the test- mass average velocity v(tj) in this time interval (tj indicates the time 
moment after the second reflection): 



where 



0''°"(tj) = 4>{tj) + 2FkpTv{tj) 



(203) 



(204) 



We omit here mathematical details of the transition to the continuous measurement limit as 



they are essentially the same as in the position measurement case, see Section 4.1.2 and start 
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Figure 26: Toy example of the quantum speedmeter scheme. 



directly with the continuous time equations. The output signal of the homodyne detector in the 
speedmeter case is described by the following equations: 



^"«^^lsinrA.^ + ^i^l:i^ 
2Io 



-0"fl(t)sin0LO + 



0(i) = 

^™fl(i) = ^{t) + 2rkp[x{t)-xit-T)]. 
In spectral representation these equations yield: 



cos 0LO , 



x{n) 



o{n) 

2rfcpSin0LO 



x{n) + ifl(ri). 



where 



£fi(^) 



1 



2rfcp(l~e*f2^) 






(205) 
(206) 

(207) 
(208) 



is the equivalent displacement measurement noise. 

The back-action force with account for the two subsequent light reflections off the faces of the 
probe, can be written as: 

Fb.a. (t) = — [i{t + r) - i(t)] (209) 



and in spectral form as: 



2r. 



i^b.a.(^)-— I(r!)(e-"'^-l). 
c 



(210) 



Then one can make a reasonable assumption that the time between the two reflections r is 
much smaller than the signal force variation characteristic time (~ fl~^) that spills over into the 
following condition: 

nT<l, (211) 



and allows one to expand the exponents in Eqs. (208 209) into a Taylor series: 

%(f]) 



ifl(O) 



-in 



fb.a.(t)=-«fipb.a.(rJ), 



(212) 
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where 



vain) = 



Pb.a.(^) = 



1 

2rfc„r 



H^) - -^-^ cot (/.LO 



fb.a.(») ^ 2FTi{n) 
—iQ c 



Spectral densities of theses noises are equal to 

by 
?]2 ' 



s.m 



Spi^) = n^Sp{n) , s^f(^) = -Sypin) , 



(213) 



(214) 



(215) 



where 



Dy 



Or, — 



-'vp 






Si 
2ujpX 



Scb 



Si 2 

^COt 0LO 



cot 0LO 



Note also that 



s^{n)SFm- 



sipm 



bybp 



s: 



vp 



S,f,Si h 

,2 - ^ 



CJ, 



(216) 



(217) 



p 



The apparent difference of the spectral densities presented in Eq. (2151 from the ones describing 



the 'ordinary' position meter (see Eqs. ( 191 )) is that they now have rather special frequency depen- 
dence. It is this frequency dependence that together with the cross-correlation of the measurement 
and back-action fluctuations, SxF ^ 0, allows the reduction of the sum noise spectral density to 



arbitrarily small values. One can easily see it after the substitution of Eq. (2151 into Eq. (144) 
with a free mass XxxiS^ = — l/(Afri^) in mind: 



2Af 5„ 



bj, 



(218) 



If there was no correlation between the back-action and measurement fluctuations, i.e., Syp — 0, 



then by virtue of the uncertainty relation, the sum sensitivity appeared limited by the SQL ( 161 1 



S' 



n' 



'iSy 



Sy ^ hMfl^ 



(219) 



One might wonder, what is the reason to implement such a complicated measurement strategy 
to find ourselves at the same point as in the case of a simple coordinate measurement? However, 
recall that in the position measurement case, a constant cross-correlation SxF oc cot 0lo allows one 
to get only a narrow-band sub-SQL sensitivity akin to that of a harmonic oscillator. This effect we 
called virtual rigidity, and showed that for position measurement this rigidity Xvirt = S^f/Sx is 
constant. In the speedmeter case, the situation is totally different; it is clearly seen if one calculates 
virtual rigidity for a speedmeter: 

SxF 






_ r)2 Syp 

bxK^i) by 



(220) 



It turns out to be frequency dependent exactly in the way that is necessary to compensate the free 
mass dynamical response to the back-action fluctuations. Indeed, in order to minimize the sum 



noise spectral density (2181 conditioned on uncertainty relation (217), one needs to set 



OxF '~*vp 



Sf 



s, 



= —^ = const. 
M 



(221) 
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which allows one to overcome the SQL simply by choosing the right fixed homodyne angle: 



cot0,o- '^';y° . (222) 



Then the sum noise is equal to 



and, in principle, can be made arbitrarily small, if a sufficient value of Sp is provided; that is, ifthe 
optomechanical coupling is sufficiently strong. 

Simple case: light in a coherent state. Let us consider how the spectral density of a speed- 
meter will appear if the light field is in a coherent state. The spectral densities of phase and power 



fluctuations are given in Eqs. (122), hence the sum noise power (double-sided) spectral density for 



the speedmeter takes the following form 

IGiOpTciF^T 



HM 


fny 


oF 
■^SQLf.m. 


2t2 


y^J 


2f]2r2 



^"(»)- .rv;L2 -E^(^) =^^^i^> (224) 



where flq for our scheme is defined in Eq. ( 157 1. This formula indicates the ability of a speedmeter 
to have a sub-SQL sensitivity in all frequencies provided high enough optical power and no optical 
loss. 

4.5.2 QND measurement of a free mass velocity 

The initial motivation to consider speed measurement rested on the assumption that a velocity 
■0 of a free mass is directly proportional to its momentum p = Mi). And the momentum in turn 
is, as an integral of motion, a QND-observable, i.e., it satisfies the simultaneous measurability 



condition (131 ) 



[p{t),p{t')]^0 yt,t'. 

But this connection between p and v holds only if one considers an isolated free mass not coupled 
to a meter. As the measurement starts, the velocity value gets perturbed by the meter and it 
is not proportional to the momentum anymore. Let us illustrate this statement by our simple 
velocity measurement scheme. The distinctive feature of this example is that the meter probes the 
test position x twice, with opposite signs of the coupling factor. Therefore, the Lagrangian of this 
scheme can be written as: 

i\/{,2 

£ = ^- + /3(t)£AA + r^eter , (225) 

where 

is the test-mass velocity, N stands for the meter's observable, which provides coupling to the test 
mass, £metor IS the self-Lagrangian of the meter, and j3{t) is the coupling factor, which has the 



form of two short pulses with the opposite signs, separated by the time r, see Figure [27j We 
suppose for simplicity that the evolution of the meter observable M can be neglected during the 
measurement (this is a reasonable assumption, for in real schemes of the speedmeter and in the 
gedanken experiment considered above, this observable is proportional to the number of optical 
quanta, which does not change during the measurement). This assumption allows one to omit the 
term £meter from consideration. 
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Figure 27: Real (/3(i)) and effective (a(t)) coupling constants in the speedmeter scheme. 



This Lagrangian does not satisfy the most well-known sufficient (but not necessary!) condition 
of the QND measurement, namely the commutator of the interaction Hamiltonian Hint = —(3{t)xJ\f 
with the operator of measured observable x does not vanish j33| . However, it can be shown that 
a more general condition is satisfied: 

[U{G,t),x] =0 

where C/(0, r) is the evolution operator of probe-meter dynamics from the initial moment istart = 
when the measurement starts till the final moment iond = t when it ends. Basically, the latter 
condition guarantees that the value of x before the measurement will be equal to that after the 
measurement, but does not say what it should be in between (see Section 4.4 of [57] for details). 

Moreover, using the following nice trick R the Lagrangian (2251 can be converted to the form, 
satisfying the simple condition of [55] : 



L' = C-'^-^N, 



where 



a{t) 



dt 



P{t') dt' , 



(227) 



(228) 



see Figure [27] It is known that two Lagrangians are equivalent if they differ only by a full time 
derivative of and arbitrary function of the generalized coordinates. Lagrangian equations of motion 
for the coordinates of the system are invariant to such a transformation. 

The new Lagrangian has the required form with the interaction term proportional to the test- 
mass velocity: 

C=^-a{t)m. (229) 

Note that the antisymmetric shape of the function j3{t) guarantees that the coupling factor a{t) 
becomes equal to zero when the measurement ends. The canonical momentum of the mass M is 

equal to 

dC 

p ^ ^- ^ Mv - a(tW , (230) 

ov 



and the Hamiltonian of the system reads 

H = pv — C 

^Personal communication with Yanbei Chen. 



[p + a{t)J^^' 
2M 



(231) 
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The complete set of observables describing our system includes in addition apart to x, p, and 
J\f, the observable $ canonically conjugated to A/": 

[^,N]^ih (232) 

(if A/" is proportional to the number of quanta in the light pulse then $ is proportional to its phase), 
which represents the output signal O of the meter. The Heisenberg equations of motion for these 
observables are the following: 

dx{t) .,^, p + a(t)J\f dp{t) 

= v[t) — , = U , 

dt M ' dt ' ^233) 

dm a{t)[p + a{t)U] dMit) 

— ; — — 7T — a(t)v{t) , — ; — = . 

These equations show clearly that (i) a canonical momentum p is preserved by this measurement 
scheme while (ii) the test-mass velocity v, as well as its kinematic momentum Mv, are perturbed 
during the measurement (that is, while a ^ 0), yet restore their initial values after the measure- 
ment, and (iii) the output signal of the meter $ carries the information about this perturbed value 
of the velocity. 

Specify a{t) to be of a simple rectangular shape: 

[0, ift<Oori>T. 
This assumption does not affect our main results, but simplifies the calculation. In this case, the 



solution of Eqs. (233) reads: 



with 



x{r) 


= X(0) + VT 


$(r) 


= ^{0)+VT 


7! 


P+M 



(234a) 
(234b) 



(235) 



the test- mass velocity during the measurement [compare with Eqs. (203 204)]. 

Therefore, by detecting the variable $(t), the perturbed value of velocity v is measured with 
an imprecision 

AUmoas = , (236) 

r 
where A<i> is the initial uncertainty of $. The test-mass position perturbation after this measure- 
ment is proportional to the initial uncertainty of A/": 

Axb.a. = Al^b.a.T = ■ (237) 

It follows from the last two equations that 

A^AA/" h ,„„„, 

A-WincasAXb.a. = Ty > ^ ■ (238) 

The sum error of the initial velocity estimate Wjnit — p/m yielding from this measurement is thus 
equal to: 

A$\^ /AA/'\^ FY 



^-- = W( — ) +{-Jf} >\jjrr- (239) 
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As we see, it is limited by a value of the velocity measurement SQL: 



A«sQL = V^. (240) 

To overcome this SQL one has to use cross-correlation between the measurement error and 
back-action. Then it becomes possible to measure iJinit with arbitrarily high precision. Such a 
cross-correlation can be achieved by measuring the following combination of the meter observables 

a>(r) - Afr/m ^ $(0) + — (241) 

instead of ^(t), which gives a sum measurement uncertainty for the initial velocity p/m, propor- 
tional to the initial uncertainty of $ only: 

AWsum - , (242) 

T 

and hence not limited by the SQL. 
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5 Quantum Noise in Conventional GW Interferometers 

FYIn Section [4j we have talked about the quantum measurement, the general structure of quan- 
tum noise implied by the quantum mechanics and the restrictions on the achievable sensitivity it 
imposes. In this section, we turn to the application of these general and lofty principles to real life, 
i.e., we are going to calculate quantum noise for several types of the schemes of GW interferometers 
and consider the advantages and drawbacks they possess. 

To grasp the main features of quantum noise in an advanced GW interferometer it would be 
elucidating to consider first two elementary examples: (i) a single movable mirror coupled to a 
free optical field, reflecting from it, and (ii) a Fabry-Perot cavity comprising two movable mirrors 
and pumped from both sides. These two systems embody all the main features and phenomena 
that also mold the advanced and more complicated interferometers' quantum noise. Should one 
encounter these phenomena in real-life GW detectors, knowledge of how they manifest themselves 
in these simple situations would be of much help in successfully discerning them. 

5.1 Movable mirror 

The scheme of the mirror is drawn in Figure |28| It is illuminated from both sides by the two 
independent laser sources with frequency Wp, and mean power values Xi and Xi. In terms of 
the general linear measurement theory of Section |4.2| we have two meters represented by these 
two incident light waves. The two arbitrary quadratures of the reflected waves are deemed as 
measured quantities 0\ and 02- Measurement can be performed, e.g., by means of two independent 
homodyne detectors. Let us analyze quantum noise in such a model keeping to the scheme given 
by Eqs. 



(1331 




■■ X 

Figure 28: Scheme of light reflection off the single movable mirror of mass M pulled by an external 
force G. 



5.1.1 Optical transfer matrix of the movable mirror 



The first one of Eqs. ( 133 ) in our simple scheme is represented by the input-output coupling 



equations (30 36) of light on a movable mirror derived in Section 2.2.4 We choose the transfer 



matrix of the mirror to be real: 



Ureal — 






(243) 



according to Eq. (29 1. Then we can write down the coupling equations, substituting electric field 



strength amplitudes E\^i and ei_2(^) by their dimensionless counterparts as introduced by Eq. (61 ) 
of Section 13.11 



Bi(r!) 

B2(17) 



= iVJlreal 



Ai(f]) 
A2(f)) 



and 



bi(fi) 
b2(0) 



= Mfcal 



di(rj) 



Ri 

R2 



x{n) , (244) 
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where 



c 






and 



R^ = 2x/i?^ 



A2. 

-A2c 



(245) 



Without any loss of generahty one can choose the phase of the hght incident from the left to 
be such that Ai^ = and Ai,, = yJIXxjihw^). Then factoring in the constant phase difference 
between the left and the right beams equal to $0i one would obtain for the left light {A2C, ^2s} = 
^2Z2/(^p){cos $0, sin $o}- 

The two output measured quantities will then be given by the two homodyne photocurrents: 



6x(Sl) = H^[4)i]bi=bic{n)cos(Pi+bis{^)sm(l)i, 
62 (f^) = -H''^[(/'2]b2 =^2c(^)cos02+f'2s(f^)sin(/)2, 



where vector i? [(/)] was first introduced in Section 2.3.2 after Eq. p7| ) as: 



COS( 

sine 



Then, after substitution of ( 243[ ) into (244), and then into (2461, one gets 

d^°\n) = H^[(bi]{-VRdi + Vfd2), 
d^^\n) = H^[(b2]iVTdi + ^d2), 

and 



XOif(^) 

X02f{^) 



H-[^,]R,^-h/^^ 



sin0i 



H' 



li?2 



-^— — ^^ sm(</.2 - $0) 
he 



(246) 

(247) 



(248) 



(249) 
(250) 



(251) 
(252) 



5.1.2 Probe's dynamics: radiation pressure force and ponderomotive rigidity 



Now we can write down the equation of motion for the mirror assuming it is pulled by a GW tidal 
force G: 

mx{t) = Fr.p. (t) + G{t) =^ -M9?x{Vt) = F^.p. (^) + G{Vt) , (253) 

that gives us the probe's dynamics equation (the third one) in (|133|: 



x{n) = xx.m[h...m + G(o)] = -^[F,,^x^) + G{n)] 



(254) 



where the free mirror mechanical susceptibility Xxx{^) = — 1/(MX1^). The term i^b.a. (^) stands for 
the radiation pressure force imposed by the light that can be calculated as 



F,..p.(f^)=Fo+Fb.a.(f^) = 



la, m + lb, (r!) - la, (n) - ih, (n) 



nKf) 



(A^Ai + B]"Bi - AJA2 - BjBa) + 2(A]"ai(t) + Bjbiit) ~ Aja2(<) - Bjb2itl^5) 



where kp = ojp/c 
hk. 



Fo = ^(A^Ai + B^Bi - AjAa - BTB2) = — (Ii -I2) - ^^^ V^^cos^o , (256) 
2 c c 
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is the regular part of the radiation pressure forcer^ It is constant and thus can be compensated 
by applying a fixed restoring force of the same magnitude but with opposite direction, which can 
be done either by employing an actuator, or by turning the mirror into a low-frequency pendulum 
with ujm ^ ^GW by suspending it on thin fibers, as is the case for the GW interferometers, that 
provides a necessary gravity restoring force in a natural way. However, it does not change the 
quantum noise and thus can be omitted from further consideration. The latter term represents a 
quantum correction to the former one 

.Fb.a.(f^) ^ hkp{Ajdi{n)+Bjbi{n)-Ala2{n)-B]b2{n)) = Fjdi{n)+Fja2in)~Kx{n) (257) 

where Fj^ j = F^ai{Q) + Fja2{i^) is the random part of the radiation pressure that depends on 
the input light quantum fluctuations described by quantum quadrature amplitudes vectors ai(f2) 
and 0.2(51) with coefficients given by vectors: 



2^2hujj,R 



y/BI[- ^^7% cos ^Q 
-^/T22sin$o 



and F-) 



2^2hujpR 



/TX|"+ VMaCos^o 
\/KL2 sin $0 



and the term —Kx{Q?j represents the dynamical back action with 

8wpVi?TZil2sin$o 



(258) 



being a ponderomotive rigidity that arises in the potential created by the two counter propagating 
light waves. Eq. (257) gives us the second of Eqs. (133). Here xff{^) = —K. 



5.1.3 Spectral densities 



We can reduce both our readout quantities to the units of the signal force G according to Eq. ( 140 1 



Of 



xfA^) 



'F{n)+G, 



6^ 



X2(P)_ 



+ j-(rj) + G 



and define the two effective measurement noise sources as 



and an effective force noise as 



o- 



(0) 



XOif(^) ' 



and X2{n) 



O 



(0) 



X02f{^) 



m) = Fil 



absorbing optical rigidity K into the effective mechanical susceptibility: 

Xxx{^) _ 1 



X. 



tm 



i + xxr.x(fi) K-Am^ 



(259) 

(260) 
(261) 

(262) 



^''Note the second term proportional to ^XxX^ cosiJ>Oi which owes its existence to the interference of the two 
traveling waves running in opposite directions. An interesting consequence of this is that the radiation pressure 
does not vanish even if the two waves have equal powers, i.e., X\ = I2, that is, in order to compensate for the 
radiation pressure force of one field on the semi-transparent mirror, the other one should not only have the right 
intensity but also the right phase with respect to the former one: 



COS<J>i 



■ 



1 /iJ Xi - X2 

2 V T ^/T^ 
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One can then easily calculate their power (double-sided) spectral densities according to Eq. (144 1: 









Qf Qf * 

'-'12 — '-'21 — 
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2Re 
2Re 



Sx.^m 



xfA^) 

Sx.^m 



where, if both lights are in coherent states that implies 

{d,{n) o a](^')) = ^irS^^cSijSin - n') 

, one can get: 
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S. 
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xf^n) 



(263) 



16ujpl2Rsm'^{(t)2-<^o) ' 



SxM^) 



■ cot ( 



Sx,^{n)^^COti(j)2 



sx,x,m^o, 

$0) ■ (264) 



Comparison of these expressions with the Eqs. (|195| shows that we have obtained the results 

mpinj 
the Sc 

4>LO- 



similar to that of the toy example in Section |4.4[ If we switched one of the pumping carriers off, 
say the right one, the resulting spectral densities for 0[{fl) would be exactly the same as in the 

RIi and 4 



simple case of Eqs. (195), except for the substitution of F Iq 



5.1.4 Full transfer matrix approach to the calculation of quantum noise spectral 
densities 

It was easy to calculate the above spectral densities by parts, distinguishing the effective measure- 
ment and back-action noise sources and making separate calculations for them. Had we considered 
a bit more complicated situation with the incident fields in the squeezed states with arbitrary 



squeezing angles, the calculation of all six of the above individual spectral densities (264 1 and 



subsequent substitution to the sum spectral densities expressions (2631 would be more difhcult 



Thus, it would be beneficial to have a tool to do all these operations at once numerically. 

It is achievable if we build a full transfer matrix of our system. To do so, let us first consider 
the readout observables separately. We start with Oi and rewrite it as follows: 



61 = H^[(^i] (-\/i?ai + xf fli-Fi^ai + VTa2 + xf HiF^da) + xf Jf^i^ilHiG 



H' 



xZH^[<tyi]RiG., 



(265) 



where we omitted the frequency dependence of the constituents for the sake of brevity and intro- 
duced a 2 X 4 full transfer matrix Mi for the first readout observable defined as 



111 = 



-Vm + xf.RiF^ Vti + xf.RiFj 



(266) 



with outer product of two arbitrary vectors a = {ai,a2}~^ and ^ — {/3i,/32}^ written in short 
notation as: 

ai/3i ai/32 

a2/3i a2/32 



T _ 



a^' = 



(267) 



78 



In a similar manner, the full transfer matrix for the second readout can be defined as: 



112 = 



n + xf.R2F^ Vm + xf,R2Fj 



(268) 



Having accomplished this, one is prepared to calculate all the spectral densities (263) at once using 
the following matrix formulas: 



sfm = 



1 H^m 



l] ■ H[cbj] + H^[cj,,] ■ M*§i„MT . H[q^^] 



2|x; 



.cff|2 

XX I 



H^[<j),]R,R]H[^^] 



ihj) = {l,2) 



(269) 



where M* = (M^^)^ and §in is the 4 x 4-matrix of spectral densities for the two input fields: 



§in = 







Ssqz [''2 I 



(270) 



with Ssqz[fi,Oi\ defined by Eq. (83 1 



Thus, we obtain the formula that can be (and, actually, is) used for the calculation of quantum 
noise spectral densities of any, however complicated, interferometer given the full transfer matrix 
of this interferometer. 



5.1.5 Losses in a readout train 

Thus far we have assumed that there is no dissipation in the transition from the outgoing light 
to the readout photocurrent of the homodyne detector. This is, unfortunately, not the case since 
any real photodetector has the finite quantum efficiency rjd < I that indicates how many photons 
absorbed by the detector give birth to photoelectrons, i.e., it is the measure of the probability 
for the photon to be transformed into the photoelectron. As with any other dissipation, this loss 
of photons gives rise to an additional noise according to the FDT that we should factor in. We 



have shown in Section 2.2.4 that this kind of loss can be taken into account by means of a virtual 
asymmetric beamsplitter with transmission coefficients y/rfd and \/l — rjd for the signal light and 
for the additional noise, respectively. This beamsplitter is set into the readout optical train as 
shown in Figure IS] and the i-th readout quantity needs to be modified in the following way: 

o'rm = y%o.(fi) + v/r^%"*(^) , (271) 

where 

ni{n) ^ H''' [(f>i\hi{n) = n.,^c(^)coS(/)i + ni_5(ri)sin</>i 

is the additional noise that is assumed to be in a vacuum state. Since the overall factor in front of 
the readout quantity does not matter for the final noise spectral density, one can redefine Oj°''''(r2) 
in the following way: 



o'r{n) = o^,in) + edn^,in). 



where 



e-d 



1 

Vd 



1 



(272) 



The influence of this loss on the final sum spectral densities ( 269 1 is straightforward to calculate 



if one assumes the additional noise sources in different readout trains to be uncorrelated. If it is 



so, then Eq. ( 269 1 modifies as follows 



gF^U^ss^ 



1 



if[0, 



H' 



lI*§i„MT . H[cj),] 



2|XS' 



ff 12 



H^[<l),]R,R]H[q^,] 



H^\<b,]R,R]H[ 



(*,J) = (1,2). 



(273) 
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Figure 29: Fabry-Perot cavity 



Now, when we have considered aU the stages of the quantum noise spectral densities calculation 
on a simple example of a single movable mirror, we are ready to consider more complicated systems. 
Our next target is a Fabry-Perot cavity. 



5.2 Fabry— Perot cavity 

The schematic view of a Fabry-Perot cavity with two movable mirrors is drawn in Figure |29| 
This simple scheme is important for at least two reasons: (i) it is the most common element for 
more sophisticated interferometer configurations, which are considered below; and (ii) the analysis 
of real-life high-sensitivity interferometers devoted, in particular, to detection of GWs, can be 
reduced to a single Fabry-Perot cavity by virtue of the 'scaling law' theorem jSHj, see Section 5.3 



A Fabry-Perot cavity consists of two movable mirrors that are separated by a distance L 
Xi + X2, where L = cr is the distance at rest with r standing for a single pass light travel time, 
and xi and X2 are the small deviations of the mirrors' position from the equilibrium. Each of 
the mirrors is described by the transfer matrix Mi. 2 with real coefficients of reflection \J R\;i and 
transmission jTx^i according to Eq. ( 243 1 . As indicated on the scheme, the outer faces of the 



mirrors are assumed to have negative reflectivities. While the intermediate equations depend on 
this choice, the flnal results are invariant to it. The cavity is pumped from both sides by two laser 
sources with the same optical frequency Wp. 

The coupling equations for the ingoing and outgoing fields at each of the mirrors are absolutely 
the same as in Section |2.2.5| The only new thing is the free propagation of light between the 
mirrors that adds two more field continuity equations to the full set, describing the transformation 
of light in the Fabry-Perot cavity. It is illuminating to write down input-output relations first in 
the time domain: 



h, it) = -^/RiEa, [t + 2xi/c) + ^/IY4l (t) 
El, (i) = ./RlE,, {t - 2x1 1 c) + VTiEa, (i) , 
E,,{t)^Ef^{t-L/c), 



^2EeM. 

(274) 

Further we use notation t = L/c for the light travel time between the mirrors. 

The frequency domain version of the above equations and their solutions are derived in Ap- 



4. W = -\/7?2^a. (i + 2.T2/C) 

Ef^{t) = ^R~2Ea,{t- 2x2/ c) + 
KM^Ef,{t-L/c). 



pendix A.l We write these I/0-relations given in Eqs. (545) in terms of complex amplitudes 



instead of 2 photon amplitudes, for the expressions look much more compact in that representa- 
tion. However, one can simplify them even more using the single-mode approximation. 



Single-mode approximation. We note that: 

(i) in GW detection, rather high-finesse cavities are used, which implies low transmittance 
coefficients for the mirrors 

^1,2 « 1 ; (275) 
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(ii) the cavities are relatively short, so their Free Spectral Range (FSR) /fsr ~ (2t) ^ is much 
larger than the characteristic frequencies of the mirrors' motion: 

|f7|r<l; (276) 

and (iii) the detuning of the pump frequency from one of the cavity eigenfrequencies: 

(5 = ojp ( 71 is an integer) (277) 

r 

is also small in comparison with the FSR: 

|(5|t<1. (278) 

In this case, only this mode of the cavity can be taken into account, and the cavity can be treated 
as a single-mode lumped resonator. Note also that, while our intermediate equations below depend 
on whether n is even or odd, the final results do not. Therefore, we assume for simplicity that n 
is even. 



Expanding the numerators and denominators of Eqs. (540 545) into Taylor series in r and 
keeping only the first non-vanishing terms, we obtain that 

Bi,2 = 7^l^2(0)Al,2+r(0)A2,l, 



and 



Er,2 = F,,2^E^^^^L_i_^, (279) 



bi,2(c^) = 7^l^2(^)al^2(c^)+r(^!)a2,l(c^)+ ^\,^. , (280) 



where 



■^ ( \ i ( \ ~( \ \/7rai(w) + V72a2(a;)+A:(f7) 

ei,2(w) = fi,2(a;)=e(w) = ^^ m)J^ ' ^ ^ 



71,2 = ^ , (282) 

7 = 71 + 72 (283) 



is the cavity half-bandwidth. 



i(Sl) = 7 - i(5 + r?) , (284) 

are the cavity left and right refiectivities and its transmittance, 

XiSl) ^ ^"^^ , (286) 

and 

X = xi + X2 (287) 

is the sum variation of the cavity length. 
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5.2.1 Optical transfer matrix for a Fabry Perot cavity 



The above optical I/0-relations are obtained in terms of the complex amplitudes. In order to 
transform them to two-photon quadrature notations, one needs to employ the following linear 
transformations: 

1. change frequency a; — >■ Wp ± 17 and rewrite the relations between the input a(uj) and output 
operators /3(w) in the form: 



(3[(jjp + Q) = f{ujp + D?)a{ujp -\-VL) — /+Q+ and 
p\ujp -n) = riujp - n)a\ujp ^n) = flal ; (288) 

where /(f2) is an arbitrary complex-valued function of sideband frequency fl; 



£ 



2. use the definition (57) to get the following relations for two-photon quadrature operators: 

(289) 



/3c(^) 



(/+ + /-) «(/+-/-) 






Applying transformations (289) to Eqs. (280), we rewrite the I/0-relations for a Fabry-Perot 
cavity in the two-photon quadratures notations: 






,2(r!)ai,2(f7) + T(r!)d24(r!) + 277i;^L(r!)x(n) , 



1 



L(^) V7i;^ai.2(f^) + y7^a2,i(^)+X(r!) 



where 



X(f7) = 






I\"qX t i i ) 



Ec = V2ReE, E3 = V2ImE, 



(290) 
(291) 

(292) 
(293) 



L(r2) 



1 7 - i^i -5 
^M [ 5 7 - zf} 
V(Vt) = i{9.)t{-9.) = (7 - inf + 5^ , 
Mi,2(r!) = 27i,2L(r!) - I , T{^) = 2V^I^L(fl) . 



Therefore, the I/0-relations in standard form read: 
'biifl) 

b2m_ 

with optical transfer matrix defined as: 



c 



ai(f7) 

02(0) 






x{n) 



Cm = 



Mi(0) T(0) 
T(0) M2(0) 



and the response to the cavity elongation x{il) defined as: 



R^^{n) = 2/fcp^/— L(0) 



-E, 
E, 



and -R2 (^) = 2fcj; 



L(0) 



-E3 
E, 



(294) 

(295) 
(296) 



(297) 



(298) 



(299) 
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Note that due to the fact that (Mpp)"^ = (Mpp) '^ the reflectivity and the transmission matrices 
Mi_2 and T satisfy the following unitarity relations: 

MiRJ + TT^ = MaK^ + TT^ = I , MiT^ + TM| = . (300) 

5.2.2 Mirror dynamics, radiation pressure forces and ponderomotive rigidity 

The mechanical equations of motion of the Fabry-Perot cavity mirrors, in spectral representation, 
are the following: 

ix^U^)hm^F,{n) + G,{9) 1 = 1,2, (301) 

where Xxx,i are the mechanical susceptibilities of the mirrors, Gi stand for any external classical 
forces that could act on the mirrors (for example, a signal force to be detected), 

are the radiation pressure forces acting on the mirrors, and Xeii ^fi, T-a.i, T-\,i are the powers of 
the waves e^, f^, etc. The signs for all forces are chosen in such a way that the positive forces are 
oriented outwards from the cavity, increasing the corresponding mirror displacement x\_i. 

In the spectral representation, using the quadrature amplitudes notation, the radiation pressure 
forces read: 

F\.1^) = -;^ (^1^2^1,2 + F^_2Fl,2 " Aj 2-^1,2 ^ Bl,2Bl,2) 



+hkp 



Ej^e.^^m + F[2/i,2(0) - A|.2di,2(r!) - Bj^b.^^m] . (303) 



The first group, as we have already seen, describes the regular constant force; therefore, we omit 
it henceforth. 

In the single-mode approximation, the radiation pressure forces acting on both mirrors are 
equal to each other: 

Fi,2(rJ) = Fb.U^) = 2hkpE^e{n) , (304) 

and the optical field in the cavity is sensitive only to the elongation mechanical mode described by 
the coordinate x. Therefore, combining Eqs. (301[), we obtain for this mode: 



X-^ mx{n) = /'b.a. m + G{Q) , (305) 

where 

XxA^) = [Xxxa{^) + XxxA^)] (306) 

is the reduced mechanical susceptibility and 

^,^. ^ x...imGi{n) + x..,2mG2{n) 

Xxxi^) 

is the effective external force. 

In the simplest and at the same time the most important particular case of free mirrors: 

XxccA^)^ W i = l,2, (308) 

the reduced mechanical susceptibility and the effective external force are equal to 

X..(f^) = -^, (309) 
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and 



G{n) = n 



Gi(n) , G2{n) 



mi 



7712 



where 

1 1 

\ 

^777i 7772 

is the effective mass of the elongation mechanical mode 



fJ- 



(310) 
(311) 



It follows from Eqs. (|291|) and (|304|) that the radiation pressure force can be written as a sum 

(312) 



of the random and dynamical back-action terms, similarly to the single mirror case: 
with the random component equal to 



^b"i.(^) = ^^^L(f7) i^iaiin) + ^^Mm 



and the ponderomotive rigidity that reads 



K{n) 



MJS 



V{Vt) ' 
We introduced here the normalized optical power 



4^fc^ 
Mt 



McL 



(313) 



(314) 



(315) 



with 



Ic = hwp\E\^ 



(316) 

standing for the mean optical power circulating inside the cavity, and M is some (in general, 
arbitrary) mass. Typically, it is convenient to make it equal to the reduced mass fi. 



Substitution of the force (312) into Eq. (305) gives the following final equation of motion 

ix^^in) + Kman) = F^^lin) + Gin) . 



(317) 



Thus, the effective mechanical susceptibility X^ {^) for ^ Fabry-Perot cavity reads: 

xf.-''''m = {x^^m + K{n))-' 



K{n) - ^f72 



(318) 



5.3 Fabry— Perot— Michelson interferometer 

In real-life high-precision experiments with mechanical test objects, interferometer schemes that 
are much more sophisticated than the simple Fabry-Perot cavity are used. In particular, the best 
sensitivity in mechanical displacement measurements is achieved by the laser GW detectors. The 
typical scheme of such a detector is shown in Figure |30[ It is this scheme which is planned to be 
used for the next generation Advanced LIGO [Mil [Ml [l39] , Advanced VIRGO [11^, and LCGT [4] 
GW detectors, and its simplified versions are (or were) in use in the contemporary first generation 
detectors: Initial LIGO [10, 5^, VIRGO [7], GEO 600 [Hi [3], and TAMA [13, 6,. 

This scheme works similar to the ordinary Michelson interferometer considered briefly in Sec- 



tion 2.1.2 The beamsplitter BS distributes the pump power from the laser evenly between the 
arms. The beams, reflected off the Fabry-Perot cavities are recombined on the beamsplitter in 
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Figure 30: Power- and signal-recycled Fabry-Perot-Michclson interferometer. 
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such a way that, in the ideal case of perfect symmetry of the arms, ah the hght goes back to the 
laser, i.e., keeping the signal ('south') port dark. Any imbalance of the interferometer arms, caused 
by signal forces acting on the end test masses (ETMs) makes part of the pumping light leak into 
the dark port where it is monitored by a photodetector. 

The Fabry-Perot cavities in the arms, formed by the input test masses (ITMs) and the end 
test masses, provide the increase of the optomechanical coupling, thus making photons bounce 
many times in the cavity and therefore carry away a proportionally-amplified mirror displacement 
signal in their phase (cf. with the F factor in the toy systems considered in Section H. The two 
auxiliary recycling mirrors: the PRM and the signal recycling (SRM) allow one to increase the 
power, circulating inside the Fabry-Perot cavities, for a given laser power, and provide the means 
for fine-tuning of the quantum noise spectral density |105l I156J , respectively. 

It was shown in |39] that quantum noise of this dual (power and signal) recycled interferometer 
is equivalent to that of a single Fabry-Perot cavity with some effective parameters (the analysis in 
that paper was based on earlier works |114l 11301 , where the classical regime had been considered) . 
Here we reproduce this scaling law theorem, extending it in two aspects: (i) we factor in optical 
losses in the arm cavities by virtue of modeling it by the finite transmissivity of the ETMs, and (ii) 
we do not assume the arm cavities tuned in resonance (the detuned arm cavities could be used, in 
particular, to create optical rigidity in non-signal-recycled configurations). 



5.3.1 Optical I/O-relations 

We start with Eqs. (( 279[ )) and (( 280 )) for the arm cavities. The notation for the field amplitudes is 
shown in Figure 30 The fields referring to the interferometer arms are marked with the subscripts 
N ('northern') and E ('eastern') following the convention of labeling the GW interferometer parts 
in accordance with the cardinal directions they are located at with respect to the drawing (up- 
direction coincides with north). In order to avoid subscripts, we rename some of the field amplitudes 
as follows: 

aiAf,B — > aAT^B , hiN.E ^^N,E 1 &2N,E ^ gN,E', (319) 



see also Figure 30 Note that the fields gN.E now describe the noise sources due to optical losses 



in the arm cavities. 



Rewrite those Eqs. (279), (2801, (281 1 that are relevant to our consideration, in these notations: 



^N,E 



— T^arni{0)'^N,E : 
1 



.(0) 



^A 



N,E ■ 



(320) 



where 



bAr,_E(w) = TZi^,.ni{fl)aN,E{^) +%im{^)gN.E{^) + 



2-^7larmA'jv,B i^) 



^arml,^ 'J 



eN,Ei^) = 



^/lliirmaN,Ei^) + ^/l2gN,Ei^) + XN,Ei^) 



4rm(^)v^ 



^ arm -^arin 

Tlarin ^ ; 12 ~. 

4r 4t 



(321) 



(322) 



Tkrm is the input mirrors power transmittance, Aarm is the arm cavities power losses per bounce, 

7arm = 7larm + 72 (323) 

is the arm cavities half-bandwidth, 

4rm(^) = 7arm -i(<5arm + ^), (324) 
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^arm IS the ariTL cavities detuning, 



and 



7^a™(f^) - 1^ - 1 , %U^) = "^^^ (325) 

X^,E = ^fcpE^^.^f_^.^(") . (326) 



Assume then that the beamsphtter is described by the matrix (32), with R = T = 1/2. Let 
Iw = CT\Y be the power recychng cavity length (the optical distance between the power recycling 
mirror and the input test masses) and Is — cts - power recychng cavity length (the optical distance 
between the signal recycling mirror and the input test masses). In this case, the light propagation 
between the recycling mirrors and the input test masses is described by the following equations for 
the classical field amplitudes: 



A 



N,E — 



/2 



C^,s = ^^i^e^^-^^, (327) 

where 

4>w,s = ^pTw,s , (328) 

are the phase shifts gained by the carrier light with frequency ujp passing through the power and 
signal recycling cavities, and the similar equations: 



^^.^(ii^) 



dw{uj)e"^^«' ±d5(w)e*"^s 
72 



c^M = i^^lM^M^e— (329) 

v2 

appling to the quantum fields' amplitudes. 

The last group of equations that completes our equations set is for the coupling of the light 
fields at the recycling mirrors: 

hw — —y/Rwi^w + vTwCw , <iw — ^jTw^w + -\ARwCw ■, ('x^f\\ 

where Rwy Tw and Rs, Tg are the reflectivities and transmissivities of the power and signal 
recycling mirrors, respectively. These equations, being linear and frequency independent, are valid 
both for the zeroth-order classical amplitudes and for the first-order quantum ones. 

5.3.2 Common and differential optical modes 

The striking symmetry of the above equations suggests that the convenient way to describe this 
system is to decompose all the optical fields in the interferometer arms into the superposition of 
the symmetric (common) and antisymmetric (differential) modes, which we shall denote by the 
subscripts + and — , respectively: 

ajvias ' bATibb^; , Sat ± e^ ^ gN±gE ,„„,^ 
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It follows from Eqs. (329) that the symmetric mode is coupled solely to the 'western' (bright) port, 
while the antisymmetric one couples exclusively to the 'southern' (dark) port of the interferometer. 
It is easy to see that the classical field amplitudes of the antisymmetric mode are equal to zero. 
For the common mode, combining Eqs. (320), (327), (3301, (331), it is easy to obtain the following 
set of equation: 



B+ = 

E+ = 



7^a„n(0)A+ , 
1 /7la 



-A, 



4rm(0) V T 

A+ = Dvye**'"' , 
Cw = B+e*"^"-, 
BvF = —yRwAw'i'yTwGw, 

J^w — V TvyAvi/ + \/ Rw^^w ■ (332) 

Its solution is equal to (only those amplitudes are provided that we shall need below): 

7^arm(0)e2''^»' - /R^ 



B 



w 



A 



w 



E. 



l-7^arm(0)V^M/e2'^«' 

1 / 7l arm V IW 



\-W ■ 



(333a) 



(333b) 



4rm(0)V T 1 - 7^arm(0)yi?H^e2»'^v 

In the differential mode, the first non-vanishing terms are the first-order quantum-field ampli- 



tudes. In this case, using Eqs. (321), (329), (330), (331), and taking into account that 



E 



N 



E, 



E+ 
V2 



(334) 



we obtain: 



b_(c.) 



e^(w) = 



7?.larm(^)a-(w) +7Irin(^)g-(w) + 



2^711^^- (^) 



,in) 



V7larma- (w) + ^g_ (w) + X_ (Q) 



4rm(f^)\/T 

a_(c^) = d5(c^)e''^^^ 

csiio) = b_(c.)e''^^^ 

bs(w) = -y^ks + \/t^'Cs{uj) , 

ds(a;) = y^as + \/Rscs{u!) , 



(335) 



where 



x^{n) 



ifcpE_|_x_(57) 



X- 



_2P 

xn - xe 



(336) 

(337) 
(338) 



The solution of this equation set is the following: 

[7^a„„(f^)e2■'"^-^ - VlhY^sm + VT^e'^^^ 
bs(w) = 

e_(a;) = 



7;rm(^)g-(t^) + 






1 - 7earm(^)/R^e2»--« 

/75e*""^V7lI^as(fi) + [1 + VR^e^'"^^][y^g-iu) + X^{Q)] 



4rm(^)V^[l - 7^a,.m(f^)/R^e2--s] 



(339) 
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Figure 31: Effective model of the dual- recycled Fabry-Perot-Michelson interferometer, consisting 
of the common (a) and the differential (b) modes, coupled only through the mirrors displacements. 



Eqs. (332 1 and (333), on the one hand, and (335) and (339), on the other, describe two almost 



independent optical configurations each consisting of the two coupled Fabry-Perot cavities as 



featured in Figure 31 'Almost independent' means that they do not couple in an ordinary linear 
way (and, thus, indeed represent two optical modes). However, any variation of the differential 
mechanical coordinate X- makes part of the pumping carrier energy stored in the common mode 
pour into the differential mode, which means a non-linear parametric coupling between these 
modes. 

5.3.3 Interferometer dynamics: mechanical equations of motion, radiation pressure 
forces and ponderomotive rigidity 

The mechanical elongation modes of the two Fabry-Perot cavities are described by the following 



equations of motion [see Eq. (305)]: 



(340) 



where fx — A//2 is the effective mass of these modes and Gtv.b are the external classical forces 



acting on the cavities end mirrors. The differential mechanical mode equation of motion (338) 
taking into account Eq. (334) reads: 

where 



FlP-(r2) ^ 2nkpE\e_{VL) 
is the differential radiation-pressure force and 

G— ^ Gn — Ge 
is the differential external force. 



(341) 
(342) 
(343) 



Equations (339) and (341 1 together form a complete set of equations describing the differential 



optomechanical mode of the interferometer featured in Figure 31 'b). Eq. (341) implies that the 
effective mass of the differential mechanical degree of freedom coincides with the single mirror 
mass: 

2/x = Af , (344) 
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which prescribes the mirrors of the effective cavity to be twice as heavy as the real mirrors, i.e. 



2M. For the same reason Eq. (334) imphes for the effective optical power a value twice as high as 



the power of light, circulating in the arm cavities: 



x^ — hUpSlii — 2xaini — 2hoJprj 



N.E ■ 



(345) 



5.3.4 Scaling law theorem 



Return to Eqs. (333) for the common mode. Introduce the following notations: 

1- 



ow 



liw 



7la 



7la 




7larmTw 

1 + 2^/Rw COS 2(j)iY + Rw ' 

2jiay:mV Rw sin 2cj)w 
1 + 2y/Rw COS 2(pw + Rw 

llW + 72 , 



(346) 

(347) 

(348) 
(349) 



In these notation, Eqs. (3331 have the following form: 

BvK 

E, 
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(350) 

(351) 
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It is easy to see that these equations have the same form as Eqs. (279) for the single Fabry-Perot 
cavity, with the evident replacement of 71 and 5 with the effective parameters 71 vi/ and 5w- The 
only difference is an additional phase shift aw- Thus, we have shown that the power recycling 
cavity formed by the PRM and the ITMs can be treated as a single mirror with some effective 



parameters defined implicitly by Eqs. (350), complemented by light propagation over length aw /kp. 
Note also that the phase shift aw can be absorbed into the field amplitudes simply by renaming 
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which yields: 



BvF — TZw^w 
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w ■ 



(354) 



The corresponding equivalent model of the common mode is shown in Figure |32[ a) . 

Taking into account that the main goal of power recycling is the increase of the power Xc — 
hjjjp\E^\'^ circulating in the arm cavities, for a given laser power Xq = fiwpj^vKp: the optimal tuning 
of the power recycling cavity corresponds to the critical coupling of the common mode with the 
laser: 

liw = 72 , 5w = Q- (355) 
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Figure 32: Common (top) and differential (bottom) modes of the dual-recycled Fabry-Perot- 
Michelson interferometer, reduced to the single cavities using the scaling law model. 



In this case, 
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Note that this regime can be achieved even with the detuned arm cavities, 6^0. 
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Consider now the differential mode quantum field amplitudes as given in Eqs. (339). Note a 



factor e"^'^^ that describes a frequency-dependent phase shift the sideband fields acquire on their 
pass through the signal recycling cavity. It is due to this frequency-dependent phase shift that 
the differential mode cannot be reduced, strictly speaking, to a single effective cavity mode, and a 
more complicated two-cavity model of Figure [ST] should be used instead. The reduction to a single 
mode is nevertheless possible in the special case of a short signal-recycling cavity, i.e., such that: 



\n\Ts « 1 



(357) 



The above condition is satisfied in a vast majority of the proposed schemes of advanced GW 
interferometers and in all current interferometers that make use of the recycling techniques [TJ [3] . 
In this case, the phase shift (ps can be approximated by the frequency-independent value: 



ujTs ~ ujpTs = <ps ■ 
It allows one to introduce the following effective parameters: 

1 - y/Rse^"l>^ 
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and to rewrite Eqs. (339) as follows 
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Figure 33: The differential mode of the dual-recycled Fabry-Perot-Michelson interferometer in 



simplified notation (364). 



where the reflectivity and transmittance of the equivalent Fabry-Perot cavity are still defined by 



the Eqs. (285), but with new effective parameters (3591, and 



as = arg 



1 



/Rse^^'I's 



(362) 



is the phase shift introduced by the signal recycled cavity. Along similar lines as in the common 
mode case, we make the following change of variables 



age 



as: 



bs. 



(363) 



Perot cavity (280). Thus, we have successfully built a single cavity model for the differential mode. 



Eqs. (360) and (3611 have exactly the same form as the corresponding equations for the Fabry- 



see Figure 32 'b). 

The mechanical equations of motion for the effective cavity are absolutely the same as for an 
ordinary Fabry-Perot cavity considered in Section |5.2| except for the new values of the effective 
mirrors' mass 2M and effective circulating power Ic = 2Iarm- Bearing this in mind, we can procede 
to the quantum noise spectral density calculation for this interferometer. 

5.3.5 Spectral densities for the Fabry Perot Michelson interferometer 

The scaling law we have derived above enables us to calculate spectral densities of quantum noise 
for a dual-recycled Fabry-Perot-Michelson featured in Figure |30] as if it were a bare Fabry-Perot 
cavity with movable mirrors pumped from one side, similar to that shown in Figure [32] 
We remove some of the subscripts in our notations, for the sake of notational brevity: 



as -J> a, 
a;_ — > x , 
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e_ — >■ e, 



^b.a. ^ p 



(364) 



b.a. ; 



(compare Figures 32 and 33 1 . We also choose the phase of the classical field E amplitude inside 

ImE = 0=^E = V2E 



the arm cavities to be zero: 



(365) 



that obviously does not limit the generality of our consideration, yet sets the reference point for 
all the classical and quantum fields' phases. 



With this in mind, we rewrite I/0-relations (360) and (361) in the two-photon quadratures 
notation: 
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where the matrices L, Mi, T are defined by Eqs. (294) and (296) 



D(f7) =X'(f])L(f7) 
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J — ifl 



oU^pX^pj^ 



Mt McL McL 

is the normahzed optical power, circulating in the interferometer arms. 



(368) 



(369) 



Suppose that the output beam is registered by the homodyne detector; see Section[2.3.1| Com- 



bining Eqs. (366) and (272), we obtain for the homodyne detector a readout expressed in units of 
signal force: 



o 



F.loss 
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Xxx (") 



-H^) + ^^, (370) 



where 



A"°""(f7) 



0(0), loss (f^) 

Xof(^) 



2Af J71 H^[(f>LoM^) 



(371) 
stands for the measurement noise, which is typically referred to as shot noise in optomechanical 
measurement with the interferometer opto-mechanical response function defined as 



Xof{^) 



2MJ71 H^[(f>wM^) 



h v{VL) 

and the back-action noise caused by the radiation pressure fiuctuations equal to 



(372) 



P{n) = h.^XQ) = VmMj 



M^VTiam + V^gi^) 
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The dynamics of the interferometer is described by the effective susceptibility x'^ {^) that 



is the same as the one given by Eq. (318) where 



K{n) 



MJS 



(374) 



is the frequency-dependent optical rigidity that has absolutely the same form as that of a single 



Fabry-Perot cavity given by Eq. (314) 



Suppose then that the input field of the interferometer is in the squeezed quantum state that 
is equivalent to the following transformation of the input fields: 



a^E.s^^[r,e]a' 



(375) 



where the squeezing matrix is defined by Eq. ( 74 ) , and the quadrature vector d™'^ corresponds to 
the vacuum state. 



Using the rules of spectral densities computation given in Eqs. (89) and (92), taking into 



account unitarity conditions (300), one can get the following expressions for the power (double- 
sided) spectral densities of the dual-recycled Fabry-Perot-Michelson interferometer measurement 
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and back-action noise sources as well as their cross-correlation spectral density: 

h \'D{n)\^ 



Sxxm 



STTm 



Sx^m 
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These spectral densities satisfy the Schrodinger-Robertson uncertainty relation: 

sxxmsMn)-\sx^m^>^ 



(376) 

(377) 

(378) 
(379) 



of the same form as in the general linear measurement case considered in Section 4.2 see Eq. ( 148[ ), 
with the exact equality in the ideal lossless case: 



72 = , % = 1 



(380) 



see Appendix |A.2| 



5.3.6 Full transfer matrix approach to calculation of the Fabry Perot Michelson 
interferometer quantum noise 

In order to compute the sum quantum noise spectral density one has to first calculate Sxxi^), 



Sj^j^{fl) and Sxj^i^) using Eqs. (376), (377), and (378) and then insert them into the general 
formula (|144p. 



However, there is another option that is more convenient from the computational point of 
view. One can compute the full quantum noise transfer matrix of the Fabry-Perot-Michelson 
interferometer in the same manner as for a single mirror in Section 5.1.4[ The procedure is rather 
straightforward. Write down the readout observable of the homodyne detector in units of signal 
force: 



0^,0,. 4^321. /,«,H-£M = 2i2) 
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where matrices Ci^2(^) can be computed using the fact that 
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which yields: 
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In the GW community, it is more common to normalize the signal of the interferometer in units 



of GW amplitude spectrum h{Q). This can easily be done using the simple rule given in Eq. (138) 
and taking into account that in our case G{fl) -^ G{il)/2: 



O'^in) = hcwm + Km = hcwm + T\UTT^TlTU-^^F^77^f^^ 



2 I h 1 
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where /igw(^) is the spectrum of the GW signal and the second term /i„(ri) stands for the sum 

1/9 

quantum noise expressed in terms of metrics variation spectrum units, i.e., in Hz ' . 
The power (double-sided) spectral density of the sum quantum noise then reads: 
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(385) 



In conclusion, we should say that the quantum noise of the Fabry-Perot-Michelson interferom- 
eter has been calculated in many papers, starting from the seminal work by Kimble et al. [94' where 
a resonance-tuned case with 6 = was analyzed, and then by Buonanno and Chen in [37. 39 , who 
considered a more general detuned case. Thus, treading their steps, we have shown that the quan- 
tum noise of the Fabry-Perot-Michelson interferometer (as well as the single cavity Fabry-Perot 
one) has the following distinctive features: 

• It comprises two effective noise sources as in any quantum linear measurement device. These 
are measurement noise X^°^^, more frequently called quantum shot noise in the GW commu- 
nity, and the back-action noise T, often referred to as quantum radiation-pressure noise. 

• These noise sources are correlated and this correlation depends not only on the homodyne 
angle (^lo or the correlations in the input light (e.g., squeezing angle in case of squeezed 
input), but also on the interferometer effective detuning S, which, according to the scaling 
law theorem, can be changed by varying signal-recycling cavity parameters. 

• The scaling law theorem also shows that changing the arm cavities' detuning is equivalent to 
the modification of the signal recycling cavity parameters in terms of effective detuning and 
bandwidth of the interferometer. 

• Another important corollary of the scaling law is that the effective bandwidths and detunings 
for the common and differential optical modes can be chosen independently, thus making it 
possible to tune the former in resonance with the pumping laser to keep as high a value of the 
circulating optical power in the arms as possible, and to detune the latter one to modify the 
test masses dynamical response by virtue of the introduction of optical rigidity that arises in 
the detuned cavity as we have shown. 

All of these features can be used to decrease the quantum noise of the interferometer and reach 
a sensitivity below the SQL in a decent range of frequencies as we show in Section [6] 
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6 Schemes of GW Interferometers with Sub-SQL Sensitivity 

6.1 Noise cancellation by means of cross-correlation 

6.1.1 Introduction 

In this section, we consider the interferometer configurations that use the idea of the cross- 
correlation of the shot and the radiation pressure noise sources discussed in Section |4.4| This 
cross-correlation allows the measurement and the back-action noise to partially cancel each other 
out and thus effectively reduce the sum quantum noise to below the SQL. 



As we noted above, Eq. (378) tells us that this cross-correlation can be created by tuning 
either the homodyne angle (pho, the squeezing angle 6, or the detuning S. In Section [4.4[ the 
simplest particular case of the frequency-independent correlation created by means of measurement 
of linear combination of the phase and amplitude quadratures, that is, by using the homodyne angle 
0LO 7^ ""/S, has been considered. We were able to obtain a narrow-band sensitivity gain at some 
given frequency that was similar to the one achievable by introducing a constant rigidity to the 
system, therefore such correlation was called effective rigidity. 

However, the broadband gain requires a frequency-dependent correlation, as it was first demon- 
strated for optical interferometric position meters |149j , and then for general position measurement 
case |HS] . Later, this idea was developed in different contexts by several authors [85] 11201 11591 [Ml 
[711 [751 11501 113]. In particular, in [94 , a practical method of creation of the frequency-dependent 
correlation was proposed, based on the use of additional filter cavities^ which were proposed to be 
placed either between the squeeze light source and the main interferometer, creating the frequency- 
dependent squeezing angle (called pre- filtering), or between the main interferometer and the ho- 
modyne detector, creating the effective frequency-dependent squeezing angle (post- filtering) . As 
we show below, in principle, both pre- and post-filtering can be used together, providing some 
additional sensitivity gain. 

It is necessary to note also an interesting method of noise cancellation proposed by Tsang and 
Caves recently |147j . The idea was to use matched squeezing; that is, to place an additional cavity 
inside the main interferometer and couple the light inside this additional cavity with the differential 
mode of the interferometer by means of an optical parametric amplifier (OPA) . The squeezed light 
created by the OPA should compensate for the ponderomotive squeezing created by back-action 
at all frequencies and thus decrease the quantum noise below the SQL at a very broad frequency 
band. However, the thorough analysis of the optical losses influence, that as we show later, are 
ruinous for the subtle quantum correlations this scheme is based on, was not performed. 

Coming back to the filter-cavities-based interferometer topologies, we limit ourselves here by 
the case of the resonance-tuned interferometer, 6 = 0. This assumption simplifies all the equations 
considerably, and allows one to clearly separate the sensitivity gain provided by the quantum noise 
cancellation due to cross-correlation from the one provided by the optical rigidity, which will be 
considered in Section [6731 

We also neglect optical losses inside the interferometer, assuming that 72 = 0. In broadband 
interferometer configurations considered here, with typical values of 7 > 10'^ s~^, the infiuence of 
these losses is negligible compared to those of the photodetector inefficiency and the losses in the 
filter cavities. Indeed, taking into account the fact that with modern high-reflectivity mirrors, the 
losses per bounce do not exceed Aarm ^ 10"'*, and the arms lengths of the large-scale GW detectors 
are equal to several kilometers, the values of 72 ^ 1 s~^, and, correspondingly, 72/7 ^ 10"'^, are 
feasible. At the same time, the value of photodetector quantum inefficiency e^ ~ 1 — ?7d ~ 0.05 
(factoring in the losses in the interferometer output optical elements as well) is considered quite 
optimistic. Note, however, that in narrow-band regimes considered in Section [6. 3[ the bandwidth 
7 can be much smaller and influence of 72 could be significant; therefore, we take these losses into 
account in Section [6^ 



96 



Using these assumptions, the quantum noises power (double-sided) spectral densities (3761 

ag explicit form: 

[cosh2r + sinh2rcos2(0- 0lo) + ^d] j 



(377) and (378) can be rewritten in the following explicit form: 
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is the convenient optomechanical coupling factor introduced in |94) . 

Eq. (381) and (3851 for the sum quantum noise and its power (double-sided) spectral density 



in this case can also be simplified significantly: 
2 I h 1 
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In Section 6.1.2 we consider the optimization of the spectral density (3911, assuming that the 



arbitrary frequency dependence of the homodyne and/or squeezing angles can be implemented. As 
we see below, this case corresponds to the ideal lossless filter cavities. In Section [6. 1.3[ we consider 
two realistic schemes, taking into account the losses in the filter cavities. 

6.1.2 Frequency- dependent homodyne and/or squeezing angles 

Classical optimization. As a reference point, consider first the simplest case of frequency 
mdependent homodyne and squeezing angles. We choose the specific values of these parameters 



following the classical optimization, which minimizes the shot noise ( 386 ) without taking into 



account the back action. Because, the shot noise dominates at high frequencies, therefore, this 
optimization gives a smooth broadband shape of the sum noise spectral density. 
It is evident that this minimum is provided by 
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61 = 0. 



In this case, the sum quantum noise power (double-sided) spectral density is equal to 



s''{n) 



2h 
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(394) 



It is easy to note the similarity of this spectral density with the ones of the toy position meter con- 



sidered above, see Eq. (173). The only significant differences introduced here are the optical losses 



and the decrease of the optomechanical coupling at high frequencies due to the finite bandwidth 7 
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of the interferometer. If ed 
evident correspondence 



and 7 ^ 57, then Eqs. (173 1 and (394 1 become identical, with the 
loF^ O — . (395) 



7T 



In particular, the spectral density (394 1 can never be smaller than the free mass SQL S^q^i ^ (il) 
(see (174)). Indeed, it can be minimized at any given frequency fl by setting 



and in this case. 



em 
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The spectral density ( 394 ) was first calculated in the pioneering work of ,43} , where the existence 



of two kinds of quantum noise in optical interferometric devices, namely the measurement (shot) 
noise and the back action (radiation pressure) noise, were identified for the first time, and it was 
shown that the injection of squeezed light with 9 = into the interferometer dark port is equivalent 
to the increase of the optical pumping power. However, it should be noted that in the presence of 
optical losses this equivalence holds unless squeezing is not too strong, e~^ > e^. 
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Figure 34: Examples of the sum quantum noise spectral densities of the classically-optimized 
(</'LO = ''''/2, = 0) resonance-tuned interferometer. 'Ordinary': J = JaLiGO) no squeezing. 
'Increased power': J = lOJaLiGO: no squeezing. 'Squeezed': J = JaLiGO, 10 dB squeezing. For all 
plots, J = 2n X 500 s^^ and ijd = 0.95. 

The noise spectral density curves for the resonance-tuned interferometer are drawn in Figure [34] 
The default parameters for this and all subsequent similar plots are chosen to be close to those 
planned for the Advanced LIGO interferometer: the value of J = JaLioo = (27r x 100)^ s~'^ corre- 
sponds to the circulating power of Xarm = 840 kW, L = A km, and M = 40 kg; the interferometer 
bandwidth 7 = 27r x 500 s^^ is close to the one providing the best sensitivity for Advanced LIGO 
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in the presence of technical noise [57]; 10 dB squeezing (e^'' — 10), which corresponds to the best 
squeezing available at the moment (2011) in the low-frequency band |104l 1151 1 1152) ): i]ci = 0.95 
can be considered a reasonably optimistic estimate for the real interferometer quantum efhciency. 
Noteworthy is the proximity of the plots for the interferometer with 10 dB input squeezing and 
the one with 10-fold increased optical power. The noticeable gap at higher frequencies is due to 
optical loss. 



Frequency dependent squeezing angle. Now suppose that the homodyne angle can be 
frequency dependent, and calculate the corresponding minimum of the sum noise spectral den- 



sity (391 ). The first term in square brackets in this equation can be rewritten as: 
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It is evident that the minimum of ( 398 1 is provided by 
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Thus, we obtaine a well-known result 94J that, using an optimal squeezing angle, the quantum 
noise spectral density can be reduced by the squeezing factor e^^J" jn comparison with the vacuum 
input case. Note, however, that the noise contribution due to optical losses remains unchanged. 
Concerning the homodyne angle ^lO) we use again the classical optimization, setting 



0LO=2. 



(402) 



In this case, the sum noise power (double-sided) spectral density and the optimal squeezing angle 
are equal to 
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and 



tan6l = /C(r2) 
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The sum quantum noise power (double-sided) spectral density ( 403 ) is plotted in Figure 35 for 
the ideal lossless case and for rjd = 0.95 (dotted line). In both cases, the optical power and the 
squeezing factor are equal to J = •/aLiGO and e^'' — 10, respectively. 



Frequency dependent homodyne angle. Suppose now that the squeezing angle corresponds 
to the classical optimization: 

0^0 (405) 
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Figure 35: Examples of the sum quantum noise power (double-sided) spectral densities of the 
resonance-tuned interferometers with frequency-dependent squeezing and/or homodyne angles 
Left: no optical losses, right: with optical losses, rjd = 0.95. 'Ordinary': no squeezing, (/)lo = 7''/2 
'Squeezed': 10 dB squeezing, 9 = 0, (/)lo = tt/S (these two plots are provided for comparison) 
Dots [pre-filtering, Eq. ( 403[)]: 1 dB squeezing, (/)lo = 7r/2, frequency-dependent squeezing angle 
Dashes [post-filtering, Eq. (408)]: 10 dB squeezing, 9 = 0, frequency-dependent homodyne angle 
Dash-dots [pre- and post-filtering, Eq. (410l]: 10 dB squeezing, frequency-dependent squeeze and 
homodyne angles. For all plots, J = JaLiGO a-nd j = 2tt x 500 s^^. 



and minimize the resulting sum noise spectral density: 
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and is equal to 



2h 



ML'^Q? 



cosh 2r -\- sinh 2r cos 2(/)lo + e^ 



2e^''cot(/)LO+/C(0)e 



2r 



/C(J7) sin (/)Lo 
!)Lo- The minimum is provided by the following dependence 



cot ( 



'LO 



1 



S^{^) = 



2h 



IC{n) 



1 



ic{n) 



(406) 



(407) 



(408) 



The sum quantum noise spectral density ( 408 1 is plotted in Figure p5l for the ideal lossless case 
and for rjd = 0.95 (dashed lines). 

Compare this spectral density with the one for the frequency-dependent squeezing angle (pre- 



filtering) case, see Eq. (403). The shot noise components in both cases are exactly equal to each 
other. Concerning the residual back-action noise, in the pre-filtering case it is limited by the 
available squeezing, while in the post-filtering case - by the optical losses. In the latter case, 
were there no optical losses, the back-action noise could be removed completely, as shown in 
Figure [35](left). For the parameters of the noise curves presented in Figure 35 (right), the post- 
filtering still has some advantage of about 40% in the back-action noise amplitude ^/S. 



Note that the required frequency dependences ( 404 ) and ( 407 ) in both cases are similar to each 



other (and become exactly equal to each other in the lossless case Cd = 0). Therefore, similar 
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setups can be used in both cases in order to create the necessary frequency dependences with 
about the same implementation cost. From this simple consideration, it is possible to conclude 
that pre-filtering is preferable if good squeezing is available, and the optical losses are relatively 
large, and vice versa. In particular, post-filtering can be used even without squeezing, r = 0. 



Frequency dependent homodyne and squeezing angles. And, finally, consider the most 
sophisticated configuration: double-filtering with both the homodyne angle 0lo &nd the squeezing 
angle 9 being frequency dependent. 

Concerning the squeezing angle, we can reuse Eqs. (400) and (401). The minimum of the 



spectral density (401) in 0lo corresponds to 



cot 0LO = 



/C(f7) 



(409) 



and is equal to 



S''{n) 



2h 



ML^n^ 



-2r 



ic{n) 



ic{n) 



(410) 



It also follows from Eqs. (400) and (409) that the optimal squeezing angle in this case is given by 



tan6' 



/C(f}), 



(411) 



It is easy to see that in the ideal lossless case the double-filtering configuration reduces to a 
post-filtering one. Really, if e^ = 0, the spectral density (410) becomes exactly equal to that for 



the post- filtering case (408), and the frequency dependent squeezing angle (411) degenerates into 
a constant value (405). However, if ed > 0, then the additional pre-filtering allows one to decrease 
more the residual back-action term. For example, if e^*" = 10 and rjd — 0.95 then the gain in the 
back-action noise amplitude \/S is equal to about 25%. 



We have plotted the sum quantum noise spectral density (410) in Figure 34 right (dash-dots) 



This plot demonstrates the best sensitivity gain of about 3 in signal amplitude, which can be 
provided employing squeezing and filter cavities at the contemporary technological level. 



Due to the presence of the residual back-action term in the spectral density ( 410 ) , there exists an 



optimal value of the coupling factor IC{n,) (that is, the optical power) which provides the minimum 
to the sum quantum noise spectral density at any given frequency Q: 



The minimum is equal to 



icin) = 



qh 



1 



Ah 



e^e 



Cde 



(412) 



(413) 



This limitation is severe. The reasonably optimistic value of quantum efficiency r]d = 0.95 that we 
use for our estimates corresponds to e^ ~ 0.23. It means that without squeezing (r = 0) one is 
only able to beat the SQL in amplitude by 



Smin 



mm 
•^SQL 



V^«0.5. 



(414) 



The gain can be improved using squeezing and, if r -> oo then, in principle, arbitrarily high 
sensitivity can be reached. But S, depends on r only as e"''/^, and for the 10 dB squeezing, only a 
modest value of 

Cmin « 0.27 (415) 



101 



can be obtained. 

In our particular case, the fact that the additional noise associated with the photodetector 
quantum inefficiency e^ > does not correlate with the quantum fluctuations of the light in the 
interferometer gives rise to this limit. This effect is universal for any kind of optical loss in the 
system, impairing the cross-correlation of the measurement and back-action noises and thus limiting 
the performance of the quantum measurement schemes, which rely on this cross-correlation. 



Noteworthy is that Eq. (410 1 does not take into account optical losses in the filter cavities. As we 
shall see below, the sensitivity degradation thereby depends on the ratio of the light absorption per 
bounce to the filter cavities length, Af/Lf. Therefore, this method calls for long filter cavities. In 
particular, in the original paper [94], filter cavities with the same length as the main interferometer 
arm cavities (4 km), placed side by side with them in the same vacuum tubes, were proposed. For 
such long and expensive filter cavities, the infiuence of their losses indeed can be small. However, 



as we show below, in Section 6.1.3 for the more practical short (up to tens of meters) filter cavities. 



optical losses thereof could be the main limiting factor in terms of sensitivity. 

Virtual rigidity for prototype interferometers. The optimization performed above can be 
viewed also in a different way, namely, as the minimization of the sum quantum noise spectral 
density of an ordinary interferometer with frequency-independent homodyne and squeezing angles, 
yet at some given frequency fig. In Section |44j this kind of optimization was considered for a 
simple lossless system. It was shown capable of the narrow-band gain in sensitivity, similar to the 
one provided by the harmonic oscillator (thus the term 'virtual rigidity'). 

This narrow-band gain could be more interesting not for the full-scale GW detectors (where 
broadband optimization of the sensitivity is required in most cases) but for smaller devices like 
the 10-m Hannover prototype interferometer [2|, designed for the development of the measurement 
methods with sub-SQL sensitivity. Due to shorter arm length, the bandwidth 7 in those devices is 
typically much larger than the mechanical frequencies 51. If one takes, e.g., the power transmissivity 
value of T > 10~^ for the ITMs and length of arms equal to L ^ 10 m, then 7 > 10^ s~^, which is 
above the typical working frequencies band of such devices. In the literature, this particular case 
is usually referred to as a bad cavity approximation. 

In this case, the coupling factor )C{Q) can be approximated as: 

ICin) « ^ , (416) 



where 



nl^-. (417) 

7 



Note that in this approximation, the noise spectral densities ( 386 1 , ( 387 1 and ( 388 ) turn out to be 
frequency independent. 



In Figure 36 the SQL beating factor 



is plotted for the sum quantum noise spectral density S'''{n) with the following values of homodyne 
and squeezing angles 

cot0LO = T^^, tane= f' JC{no), (419) 



and factoring in the "bad cavity" condition (|416). The four panes correspond to the following 



four combinations: (upper left) no losses {rjd — 1) and no squeezing [r — 0); (lower left) no losses 
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Figure 36: Plots of the locally-optimized SQL beating factor S,{fl) (418) of the interferometer 
with cross-correlated noises for the "bad cavity" case Jlo ^ 7, for several different values of the 
optimization frequency JIq within the range 0.1 x Qq < flo < \/T0 x fi^. Thick solid lines: the 



common envelopes of these plots; see Eq. (420). Left column: rjd ~ 1; right column: 7]d = 0.95. 
Top row: no squeezing, r — 0; bottom row: 10 dB squeezing, e^^ — 10. 
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iVd = 1) and 10 dB squeezing (e^'' = 10); (upper right) with losses (77^ = 0.95) and no squeezing 
(r — 0); (lower right) with losses (77^ — 0.95) and 10 dB squeezing(e^'" = 10). In each pane, the 
family of plots is shown that corresponds to different values of the ratio flo/flq, ranging from 0.1 
to yfo. 



The minima of these plots form the common envelope, given by Eqs. (410) and (416) 

1 



Cmin(^o) — 



q 



^ ^'"'ni + l + e^e^r^l 



(420) 



which is also plotted in Figure [36] It is easy to see that in the ideal case of e^ = 0, there is no 
limitation on the SQL beating factor, provided a sufficiently small ratio of Qo/flq: 



^min(f^) 



n^n 



2 ni 



However, if Cd > 0, then function (4201 has a minimum in flq at 



^n — f^o 



1 



fide' 



(421) 



(422) 



compare with Eq. (412|], equal to (4131 



6.1.3 Filter cavities in GW interferometers 

Input /output relations for the filter cavity. In essence, a filter cavity is an ordinary Fabry- 
Perot cavity with one partly transparent input/output mirror. The technical problem of how to 
spatially separate the input and output beam can be solved in different ways. In the original 
paper [94 1 the triangular cavities were considered. However, in this case, an additional mirror in 
each cavity is required, which adds to the optical loss per bounce. Another option is an ordinary 
linear cavity with additional optical circulator, which can be implemented, for example, by means 
of the polarization beamsplitter and Faraday rotator (note that while the typical polarization optics 
elements have much higher losses than the modern high-quality mirrors, the mirrors losses appear 
in the final expressions inflated by the filter cavity finesse). 

In both cases, the filter cavity can be described by the input/output relation, which can be 
easily obtained from Eqs. (290) and (291) by setting Ec,s = (there is no classical pumping in 
the filter cavity and, therefore, there is no displacement sensitivity) and by some changes in the 
notations: 

6{n) = Rf{nyi{n) + Tf{n)ci{n), (423) 

where i and 6 are the two-photon quadrature amplitude vectors of the input and output beams, 
q stands for noise fields entering the cavity due to optical losses (which are assumed to be in a 
vacuum state). 



Rf{n) 
Vf{n) 

7/1 

7/2 



1 



Vf{n) 



^h 



^h 



2V7/17/2 



7/ — ifl 

(7/ - in)' + 6} , 

AL 



-6} 
27/1(5/ 

"ff — in 



n'^ + 2in-ff2 



1% 



-fh 



-27/1(5/ 
5j + 02 



2in^f2 



/ 

cA/ 
ALf 



(424) 

(425) 

(426) 
(427) 

(428) 



104 



where Tf is the power transmittance of the input/output mirror, Af is the factor of power loss per 
bounce, L/ is the filter cavity length, 

7/ = 7/1 + 7/2 (429) 

is its half-bandwidth, and 5f is its detuning. 

In order to demonstrate how the filter cavity works, consider the particular case of the lossless 
cavity. In this case, 

6(f7) =M/(f7)i(f]), (430) 

and the reflection matrix describes field amplitude rotations with the frequency-dependent rotation 
angle: 

M/(f)) =P[6l/(f7)]e*'^^(") , (431) 



where 



9/(f}) = arctan , ^^^ { ^^ , (432) 



„^^/(f2) _ 



|P/(^)I 



(433) 



The phase factor e'*^^'^^^ is irrelevant, for it does not appear in the final equations for the spectral 
densities. 

Let us now analyze the influence of the filter cavities on the interferometer sensitivity in post- 
and pre-filtering variational schemes. Start with the latter one. Suppose that the light, entering 
the interferometer from the signal port is in the squeezed state with fixed squeezing angle 9 and 
squeezing factor r and thus can be described by the following two-photon quadrature vector 

a = §,q,[r,0]a™% (434) 

where the quadratures vector d™'^ describes the vacuum state. After reflecting off the filter cavity, 
this light will be described with the following expression 

6(0) =P[0/(f^)]§,q,[r,0]d™=e^^/(") =Ssq4r,0-f 0/(f^)]a™^' , (435) 

[see Eq. (75)], where a'"'"^ = P[0/(rj)]d™'^e*'^^(^^^ also describes the light field in a vacuum state. 



Thus, the pre-filtering indeed rotates the squeezing angle by a frequency-dependent angle 0/(f2). 
In a similar manner, we can consider the post-filtering schemes. Consider a homodyne detection 



scheme with losses, described by Eq. (272). Suppose that prior to detection, the light described by 
the quadrature vector b, reflects from the filter cavity. In this case, the photocurrent (in Fourier 
representation) is proportional to 

i.in) (X H''[<j,i^o\[V[ef{n)\b{n)e'P'^^^+ean{n)\ 

= H^[^^o-Ofm][b{n) + edfi'{n)yf'i'^''\ (436) 

where fi' — V\—9f{Vt)\he~^^f'^^^ again describes some new vacuum field. This formula demon- 
strates that post-filtering is equivalent to the introduction of a frequency-dependent shift of the 
homodyne angle by —0/(17). 

It is easy to see that the necessary frequency dependencies of the homodyne and squeezing 
angles (|404) or (407) (with the second-order polynomials in Sl^ in the r.h.s. denominators) cannot 



be implemented by the rotation angle (432[) (with its first order in i7^ polynomial in the r.h.s. 



denominator). As was shown in the paper 94 , two filter cavities are required in both these cases. 
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In the double pre- and post-filtering case, the total number of the filter cavities increases to four. 
Later it was also shown that, in principle, arbitrary frequency dependence of the homodyne and/or 
squeezing angle can be implemented, providing a sufficient number of filter cavities [40 . 

However, in most cases, a more simple setup consisting of a single filter cavity might suffice. 
Really, the goal of the filter cavities is to compensate the back-action noise, which contributes 
significantly in the sum quantum noise only at low frequencies 17 < 51^ = ^1J j^. However, when 



7 



> J^l\ 



(437) 



which is actually the case for the planned second and third generation GW detectors, the factor 



/C(f2) can be well approximated by Eq. (4161 in the low-frequency region. In such a case the single 
filter cavity can provide the necessary frequency dependence. Moreover, the second filter cavity 
could actually degrade sensitivity due to the additional optical losses it superinduces to the system. 
Following this reasoning, we consider below two schemes, each based on a single filter cavity 
that realize pre-filtering and post-filtering, respectively. 



Single-filter cavity-based schemes. The schemes under consideration are shown in Figure 37 



In the pre-filtering scheme drawn in the left panel of Figure 37 a squeezed light source emits 
frequency-independent squeezed vacuum towards the filter cavity, where it gets reflected, gaining a 
frequency-dependent phase shift 6*^(51), and then enters the dark port of the main interferometer. 
The light going out of the dark port is detected by the homodyne detector with fixed homodyne 
angle (/)lo in the usual way. 



Following the prescriptions of Section 6.1.2 we suppose the homodyne angle defined by Eq. (402 1 



The optimal squeezing angle should then be equal to zero at higher frequencies, see (404 1. Taking 



into account that the phase shift introduced by the filter cavity goes to zero at high frequencies, 
we obtain that the squeezing angle d of the input squeezed vacuum must be zero. Combining 



Eqs. (390) and (423) taking these assumptions into account, we obtain the following equation for 



the sum quantum noise of the pre-filtering scheme: 



/isum(^) = -y J^TiTTT ^^N2]{(7+zr!)K(f])[M/(f])§,q,[r,0]a™^(f])-fT/(r!)q(r!)]+(7-zri)edn(ri) 



which yields the following expression for a power (double-sided) spectral density 
2h 



(438) 



S''{VL) = 



AfL202^^{^^['^/2Wfi)[%(^)§sqz[2r,0]M}(fi)+T;(f])T}(f7)]Kt(f7)/f[V2]+e^}. 

(439) 



In the ideal lossless filter cavity case, taking into account Eq. (4311, this spectral density can 
be simplified as follows: 



s'^in) = 



2h 



ML^n^ic{n) 



[H^[Tr/2]Kin)S{2r,9f{n))K\n)H[TT/2]+el] 



(440) 



[compare with Eq. (391)]. In this case, the necessary frequency dependence of the squeezing 
angle (404) can be implemented by the following filter cavity parameters: 



where 




(441) 
(442) 
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Figure 37: Schemes of interferometer with the single filter cavity. Left: In the pre-filtering scheme, 
squeezed vacuum from the squeezor is injected into the signal port of the interferometer after the 
reflection from the filter cavity; right: in the post-filtering scheme, a squeezed vacuum first passes 
through the interferometer and, coming out, gets refiected from the filter cavity. In both cases the 
readout is performed by an ordinary homodyne detector with frequency independent homodyne 
angle (f>i^o- 



Along similar lines, the post-filtering scheme drawn in the right panel of Figure 37 can be 
considered. Here, the squeezed-vacuum produced by the squeezor first passes through the interfer- 
ometer and then, coming out, gets reflected from the filter cavity, gaining a frequency-dependent 
phase shift, which is equivalent to introducing a frequency dependence into the homodyne angle, 
and then goes to the fixed angle homodyne detector. Taking into account that this equivalent 
homodyne angle at high frequencies has to be tt/2, and that the phase shift introduced by the 
filter cavity goes to zero at high frequencies, we obtain that the real homodyne angle must also be 
7r/2. Assuming that the squeezing angle is defined by Eq. (405) and again using Eqs. (390) and 



(423), we obtain that the sum quantum noise and its power (double-sided) spectral density are 



equal to 






L]] 2MJj H^[TT/2]Rf{n){°) 
iO)M/(r2)K(f7)§sqz[r, 0]d™'=(f7) + (7 - ifl) [T/(r2)q(fi) + edn{n)] } (443) 
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and 



s'\n) = 



2% 



ML2f72/C(f7)|jjT[^/2]M/(j^)(0)|2 
X {i?"r[7r/2] [M/(rj)K(f})Ssqz[2r, 0]Kt(n)M|(17) + T/(f7)Tj.(rj)] i?[7r/2] + e^}(444) 
In the ideal lossless filter cavity case, factoring in Eq. (|431[), this spectral density takes a form 



similar to (391), but with the frequency-dependent homodyne angle: 
2h 



s''{n) = 



ML^n^ic{n)sm''(j)i^o{n) 



H^[cl)i^om]K{n)H,^,[2r,0]KHn)H''[<j>^om] + e^] , (445) 



where 



(t)wm^7r/2~9fin). 



(446) 



The necessary frequency dependence ( 407 ) of this effective homodyne angle can be implemented 



by the following parameters of the filter cavity: 

7/ =Sf = 



7/0 



yr 



eje-^r 



(447) 



Note that for reasonable values of loss and squeezing factors, these parameters differ only by a few 
percents from the ones for the pre-filtering. 



It is easy to show that substitution of the conditions (441) and (447) into Eqs. (440) and 



(445), respectively, taking Condition (437) into account, results in spectral densities for the ideal 



frequency dependent squeezing and homodyne angle, see Eqs. ( 403 ) and ( 408 ) 



In the general case of lossy filter cavities, the conditions (404) and (407) cannot be satisfied 



exactly by a single filter cavity at all frequencies. Therefore, the optimal filter cavity parameters 
should be determined using some integral sensitivity criterion, which will be considered at the end 
of this section. 

However, it would be a reasonable assumption that the above consideration holds with good 
precision, if losses in the filter cavity are low compared to other optical losses in the system: 



7/2 

7/ 



7/0 



(448) 



This inequality can be rewritten as the following condition for the filter cavity specific losses: 



47/0 



^/ < -* //o 2 
Lf c 



(449) 



In particular, for our standard parameters used for numerical estimates (J = JaLiGO: 7 = 27r x 
500 s^^, rjd = 0.95), we obtain 7/0 ~ 280 s^^, and there should be 



Af 



< 2 X 10~^ 



(450) 



(the r.h.s. corresponds, in particular, to a 50 m filter cavity with the losses per bounce Af = 10~^). 
Another more crude limitation can be obtained from the condition that 7^2 should be small 
compared to the filter cavity bandwidth 7/: 



Af_ ^ 47/0 
Lf ^ c 



(451) 
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Apparently, were it not the case, the filter cavity would just cease to work properly. For the same 
numerical values of J and 7 as above, we obtain: 



^^ < 4 X 10"*^ m~ 



(452) 



(for example, a very short 2.5 m filter cavity with Af = 10 '^ or 25 m cavity with Af = 10 ^). 



Numerical optimization of filter cavities. In the experiments devoted to detection of small 
forces, and, in particular, in the GW detection experiments, the main integral sensitivity measure is 
the probability to detect some calibrated signal. This probability, in turn, depends on the matched 
filtered SNR defined as 

|/i^(r2)|2 dn 
, s''{n) 2^ 



(453) 



with hs (r^) the spectrum of this calibrated signal. 

In the low and medium frequency range, where back-action noise dominates, and wherein our 
interest is focused, the most probable source of signal is the gravitational radiation of the inspiraling 
binary systems of compact objects such as neutron stars and/or black holes |1331I126] . In this case, 
the SNR is equal to (see [55] ) 



fco 



(454) 



S''{n) 27r 

where fco is a factor that does not depend on the interferometer parameters, and /max is the cut-off 
frequency that depends on the binary system components' masses. In particular, for neutron stars 
with masses equal to 1.4 solar mass, /max ~ 1-5 kHz. 

Since our goal here is not the maximal value of the SNR itself, but rather the relative sensitivity 
gain offered by the filter cavity, and the corresponding optimal parameters 7/1 and Sf, providing 
this gain, we choose to normalize the SNR by the value corresponding to the ordinary interferometer 
(without the filter cavities): 



pI 



fco 



2x/„_ ^-7/3 ^^ 



Sl^ifl) 27r 



with power (double-sided) spectral density 



s^{n) = 



2h 



ML^n^ 



/C(^) 



IC{Q) 



(455) 



(456) 



[see Eq. ( |394| )]. 

We optimized numerically the ratio p^ /p^, with filter cavity half-bandwidth 7/1 and detuning 
6f as the optimization parameters, for the values of the specific loss factor Af/Lf ranging from 
10~^ (e.g., very long 10 km filter cavity with Af — 10~^) to 10~^ (e-g-i 10 m filter cavity with 
Af — 10"^). Concerning the main interferometer parameters, we used the same values as in all 



our previous examples, namely, J = JaLiGO: 7 = 27r x 500 s 



The results of the optimization are shown in Figure 38 
of the filter cavity parameters 7/1 and Sf 



, and rid = 0.95. 
In the left pane, the optimal values 



are plotted, and in the right one the corresponding 
optimized values of the SNR are. It follows from these plots that the optimal values of 7/1 and Sf 



are virtually the same as 7/0, while the specific loss factor Af/Lf satisfies the condition (448), and 
starts to deviate sensibly from 7/0 only when Af/Lf approaches the limit (451 ). Actually, for such 
high values of specific losses, the filter cavities only degrade the sensitivity, and the optimization 
algorithm effectively turns them off, switching to the ordinary frequency-independent squeezing 
regime (see the right-most part of the right pane). 
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Figure 38: Left: Numerically-optimized filter-cavity parameters for a single cavity based pre- and 
post-filtering schemes: half-bandwidth 7^1 (solid lines) and detuning b^ (dashed lines), normalized 
by 7/0 [see Eq. (442)], as functions of the filter cavity specific losses Af/Lf. Right: the corre- 



sponding optimal SNRs, normalized by the SNR for the ordinary interferometer [see Eq. (455)]. 
Dashed lines: the normalized SNRs for the ideal frequency-dependent squeeze and homodyne angle 
cases, see Eqs. (403) and (408). 'Ordinary squeezing': frequency-independent 10 dB squeezing with 
9 = 0. In all cases, J = JaLiGO, 7 = 27r x 500 s~^, and r/^ = 0.95. 



It also follows from these plots that post-filtering provides slightly better sensitivity, if the 
optical losses in the filter cavity are low, while the pre-filtering has some advantage in the high-losses 
scenario. This difference can be explained in the following way [91]. The post-filtration effectively 
rotates the homodyne angle from 0lo = 7''/2 (phase quadrature) at high frequencies to 0lo — ^ 
(amplitude quadrature) at low frequencies, in order to measure the back-action noise, which 
dominates the low frequencies. As a result, the optomechanical transfer function reduces at low 
frequencies, emphasizing all noises introduced after the interferometer [see the factor sin^ (/)Lo(f^) 
in the denominator of Eq. (445)]. In the pre-filtering case there is no such effect, for the value of 
(I'LO = 7r/2, corresponding to the maximum of the optomechanical transfer function, holds for all 
frequencies (the squeezing angle got rotated instead). 

The optimized sum quantum noise power (double-sided) spectral densities are plotted in Fig- 
ure |39] for several typical values of the specific loss factor, and for the same values of the rest of 
the parameters, as in Figure 38 For comparison, the spectral densities for the ideal frequency- 



dependent squeezing angle Eqs. (403) and homodyne angle (408) are also shown. These plots 
clearly demonstrate that providing sufficiently-low optical losses (say, Aj/Lf < 10~®), the single 
filter cavity based schemes can provide virtually the same result as the abstract ones with the ideal 
frequency dependence for squeezing or homodyne angles. 



6.2 Quantum speedmeter 

6.2.1 Quantum speedmeter topologies 

A quantum speedmeter epitomizes the approach to the broadband SQL beating, in some sense, 
opposite to the one based on the quantum noises cross-correlation tailoring with filter cavities. 
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Figure 39: Examples of the sum quantum noise power (double-sided) spectral densities of the 
resonance-tuned interferometers with the single filter cavity based pre- and post-filtering. Left: 
pre-filtering, see Figure [37](left) ; dashes - 10 dB squeezing, (j)i,o = 7r/2, ideal frequency-dependent 

- 10 dB squeezing, (/)lo = '"'/S 
IQ-^ 10-'^ 10-6-5 and 10" " 



squeezing angle (404 1; thin solid 



numerically-optimized lossy 
Right: post-filtering, see Fig- 



0, ideal frequency-dependent homodyne angle (4071; 



pre-filtering cavity with Af/Lf 
urel37](right); dashes: 10 dB squeezing, 
thin solid - 10 dB squeezing, 9 = 0, numerically optimized lossy post-filtering cavity with 
Af/Lf = 10^^, 10^^ and 10~^. In both panes (for the comparison): 'Ordinary' - no squeez- 
ing, 0LO — tt/Z; 'Squeezed': 10 dB squeezing, 9 = 0, <j)LO = 7r/2 



considered above. Here, instead of fitting the quantum noise spectral dependence to the Fabry- 



Perot-Michelson interferometer optomechanical coupling factor (389), the interferometer topology 



is modified in such a way as to mold the new optomechanical coupling factor /CsM(f^) so that 
it turns out frequency-independent in the low- and medium-frequency range, thus making the 
frequency-dependent cross-correlation not necessary. 

The general approach to speed measurement is to use pairs of position measurements separated 
by a time delay r < l/fi, where f2 is the characteristic signal frequency (cf. the simplified consid- 
eration in Section [Xs] ). Ideally, the successive measurements should be coherent, i.e., they should 
be performed by the same photons. In effect, the velocity v of the test mass is measured in this 
way, which gives the necessary frequency dependence of the /Csm(^)- 

In Section [4. 5[ we have considered the simplest toy scheme that implements this principle and 
which was first proposed by Braginsky and Khalili in |26| . Also in this paper, a modified version 
of this scheme, called the sloshing- cavity speedmeter, was proposed. This version uses two coupled 
resonators (e.g., microwave ones), as shown in Figure 40 (left), one of which (2), the sloshing cavity, 



is pumped on resonance through the input waveguide, so that another one (1) becomes excited at 
its eigenfrequency Wg. The eigenfrequency of resonator 1 is modulated by the position x of the test 
mass and puts a voltage signal proportional to position x into resonator 2, and a voltage signal 
proportional to velocity dx/dt into resonator 1. The velocity signal flows from resonator 1 into an 
output waveguide, from which it is monitored. One can understand the production of this velocity 
signal as follows. The coupling between the resonators causes voltage signals to slosh periodically 
from one resonator to the other at frequency Vt. After each cycle of sloshing, the sign of the signal 
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is reversed, so the net signal in resonator 1 is proportional to the difference of the position at times 
t and t + 27r/J7, thus implementing the same principle of the double position measurement. 
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Figure 40: Left: schematic diagram of the microwave speedmeter on coupled cavities as given 
in [2 6) . Right: optical version of coupled-cavities speedmeter proposed in fT 



Later, the optical version of the sloshing-cavity speedmeter scheme suitable for large-scale laser 
GW detectors was developed [2S1 11281 I129J . The most elaborated variant proposed in |129) is 



shown in Figure 40 (right). Here, the differential mode of a Michelson interferometer serves as the 
resonator 1 of the initial scheme of j21] , and an additional kilometer-scale Fabry-Perot cavity - as 
the resonator 2, thus making a practical interferometer configuration. 

In parallel, it was realized by Chen and Khalili [47, 89J that the zero area Sagnac interfer- 
ometer |1421 [TBI 1146] actually implements the initial double-measurement variant of the quantum 
speedmeter, shown in Figure |26[ Further analysis with account for optical losses was performed 
in |55| and with detuned signal-recycling in |115| . Suggested configurations are pictured in Fig- 
ure 41 The core idea is that light from the laser gets split by the beamsplitter (BS) and directed 



to Fabry-Perot cavities in the arms, exactly as in conventional Fabry-Perot-Michelson interfer- 
ometers. However, after it leaves the cavity, it does not go back to the beamsplitter, but rather 
enters the cavity in the other arm, and only afterwards returns to the beamsplitter, and finally 
to the photo detector at the dark port. The scheme of [JT] uses ring Fabry-Perot cavities in the 
arms to spatially separate ingoing and outgoing light beams to redirect the light leaving the first 
arm to the second one evading the output beamsplitter. The variant analyzed in [89t [55] uses po- 
larized optics for the same purposes: light beams after ordinary beamsplitter, having linear (e.g., 
vertical) polarization, pass through the polarized beamsplitter (PBS), then meet the A/4 plates 
that transform their linear polarization into a circular one, and then enter the Fabry-Perot cavity. 
After reflection from the Fabry-Perot cavity, light passes through a A/4-plate again, changing its 
polarization again to linear, but orthogonal to the initial one. As a result, the PBS reflects it and 
redirects to another arm of the interferometer where it passes through the same stages, restoring 
flnally the initial polarization and comes out of the interferometer. With the exception of the 
implementation method for this round-robin pass of the light through the interferometer, both 
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Figure 41: Two possible optical realizations of zero area Sagnac speedmeter. Left panel: The ring 
cavities can be used to spatially separate the ingoing fields from the outgoing ones, in order to 
redirect output light from one arm to another [47]. Right panel: The same goal can be achieved 
using an optical circulator consisting of the polarization beamsplitter (PBS) and two A/4-plates [89l 



schemes have the same performance, and the same appellation Sagnac speedmeter will be used for 
them below. 

Visiting both arms, counter propagating light beams acquire phase shifts proportional to a sum 
of end mirrors displacements of both cavities taken with time delay equal to average single cavity 
storage time Tami: 



bji ex XN{t) +XE{t + Tarm) , 



CX XE{t) +XNit + Ti^rm) ■ 



(457) 



After recombining at the beamsplitter and photo detection the output signal will be proportional 
to the phase difference of clockwise (R) and counter clockwise (L) propagating light beams: 

S(I)R - S(j>L OC [xN{t) - XN{t + Tarm)] " [xE[t) - XE{t + Tarm)] OC XN[t) - XE{t) + O^T^rm) (458) 

that, for frequencies ^ t~1^^, are proportional to the relative velocity of the interferometer end test 
masses. 

Both versions of the optical speedmeter, the sloshing cavity and the Sagnac ones, promise about 
the same sensitivity, and the choice between them depends mostly on the relative implementation 
cost of these schemes. Below we consider in more detail the Sagnac speedmeter, which does not 
require the additional long sloshing cavity. 

We will not present here the full analysis of the Sagnac topology similar to the one we have 
provided for the Fabry-Perot-Michelson one. The reader can find it in [Ul [5S] . We limit ourselves 
by the particular case of the resonance tuned interferometer (that is, no signal recycling and 
resonance tuned arm cavities). It seems that the detuned Sagnac interferometer can provide a 
quite interesting regime, in particular, the negative inertia one |115| . However, for now (2011) the 
exhaustive analysis of these regimes is yet to be done. We assume that the squeezed light can be 
injected into the interferometer dark port, but consider only the particular case of the classical 
optimization, 9 — 0, which gives the best broadband sensitivity for a given optical power. 
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6.2.2 Speedmeter sensitivity, no optical losses 



In order to reveal the main properties of the quantum speedmeter, start with the simplified case of 
the lossless interferometer and the ideal photodetector. In this case, the sum quantum noise power 
(double-sided) spectral density of the speedmeter can be written in a form similar to the one for 
the Fabry-Perot-Michelson interferometer [see, e.g., Eqs. (386), (387) and (388)]: 



s^^n) = ^^^ 



e-^'' + e^'- cot^ (/)Lo 



2e2''cot0LO+^SM(^)e^ 



/Csm(^) 

but with a different form of the optomechanical coupling factor, see |47j : 

4J7 



A^sm(^) 



(72 + n^)2 



(459) 



(460) 



The factor J here is still defined by Eq. (369), but the circulating power is now twice as high as 



that of the position meter, for the given input power, because after leaving the beamsplitter, here 
each of the "north" and "east" beams visit both arms sequentially. 

The key advantage of speedmeters over position meters is that at low frequencies, fi < 7, /Csm 
is approximately constant and reaches the maximum there: 



4J 
JCsm(^ ^ 7) ~ arccot /CsM (0) — arccot ^5- 



(461) 



As a consequence, a frequency-independent readout quadrature optimized for low frequencies can 

be used: 

4J 
0LO = arccot /CsM (0) = arccot -^ , (462) 

which gives the following power (double-sided) spectral density 

;-2r r!4(272+rj2)2 



s''{n) 



'^'sQLV^) 



^sm(^) 



ryO 



lCsM{^)e^ 



(463) 



Here, the radiation-pressure noise (the second term in brackets) is significantly suppressed in low 
frequencies (fi ;6 7), and S'smlf '^^^ beat the SQL in a broad frequency band. 

This spectral density is plotted in Figure 42 (left). For comparison, spectral densities for the 
lossless ordinary Fabry-Perot-Michelson interferometer without and with squeezing, as well as for 



the ideal post-filtering configuration [see Eq. (408)] are also given. One might conclude from these 
plots that the Fabry-Perot-Michelson interferometer with the additional filter cavities is clearly 
better than the speedmeter. However, below we demonstrate that optical losses change this picture 
significantly. 



6.2.3 Optical losses in speedmeters 

In speedmeters, optical losses in the arm cavities could noticeably affect the sum noise at low 
frequencies, even if 

^'"«62, (464) 



72 
7i 



because the radiation pressure noise component created by the arm cavity losses has a frequency 
dependence similar to the one for position meters (remember that /C//Csm — >■ 00 if fi — ^ 0; see 
Eqs. (389), (460)). In this paper, we will use the following expression for the lossy speedmeter sum 
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Figure 42: Examples of the sum quantum noise power (double-sided) spectral densities of the 
Sagnac speedmeter interferometer (thick solid line) in comparison with the Fabry-Perot-Michelson 
based topologies considerd above (dashed lines). Left: no optical losses, right: with optical losses, 
Tjd = 0.95, the losses part of the bandwidth 72 = 1.875 s~^ (which corresponds to the losses 
^arm = 10~* per bouuce in the 4 km length arms). "Ordinary": no squeezing, (^lo = 7r/2. 
"Squeezed": 10 dB squeezing, = 0, (/)lo = '''/2. "Post-filtering": 10 dB squeezing, = 0, 
ideal frequency-dependent homodyne angle [see Eq. ( |408[ )]. For the Fabry-Perot-Michelson-based 
topologies, J = JaLiGO and 7 = 27r x 500 s~^. In the speedmeter case, J = 2JaLiGO and the 
bandwidth is set to provide the same high-frequency noise as in the other plots (7 = 27r x 385 s^^ 
in the lossless case and 7 = 27r x 360 s^^ in the lossy one). 
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noise, which takes these losses into account (more detailed treatment of the lossy speedmeter can 
be found in papers |1291l3F| ): 



s''{n) = 



^sV(^)/ 



1 



\/CsM(f^) 



+ e'"' cot-" (/.LO + — 



sm 0LO 



-2e^^ cot <^LO + ICsume^'' + el,^IC{n) 



The low-frequency optimized detection angle, in presence of loss, is 



»LO = arccot 



/Csm(O) 



= arccot 



4J/73 



l + ele-^r 



which gives 



s'\n) = ^^^ 



/Csm(^) 



/CsM(fi) 1 + ' 



n^{2-f^ + n^)^ 



„2r 



M^) 



(465) 



(466) 



(467) 



Eq. (4081]. 



[compare with Eq. ( 463 1 and note the additional residual back-action term similar to one in 



This spectral density is plotted in Figure 42 (right), together with the lossy variants of the same 

0.95, the 
— TO"'' 



configurations as in Figure 42 (left), for the same moderately optimistic value of rjd - 
losses part of the bandwidth and for 72 = 1.875 s""'^ [which corresponds to the losses A, 



per bounce in the 4 km length arms, see Eq. (322)]. These plots demonstrate that the speedmeter 
in more robust with respect to optical losses than the filter cavities based configuration and is able 
to provide better sensitivity at very low frequencies. 

It should also be noted that we have not taken into account here optical losses in the filter cavity. 
Comparison of Figure |42] with Figure [39) where the noise spectral density for the more realistic 
lossy-filter-cavity cases are plotted, shows that the speedmeter has advantage over, at least, the 
short and medium length (tens or hundred of meters) filter cavities. In the choice between very 
long (and hence expensive) kilometer scale filter cavities and the speedmeter, the decision depends, 
probably, on the implementation costs of both configurations. 



6.3 Optical rigidity 
6.3.1 Introduction 



We have seen in Section |4.3| that the harmonic oscillator, due to its strong response on near- 
resonance force, is characterized by the reduced values of the effective quantum noise and, therefore. 



by the SQL around the resonance frequency, see Eqs. (165 172 1 and Figure 22 However, practical 



implementation of this gain is limited by the following two shortcomings: (i) the stronger the 
sensitivity gain, the more narrow the frequency band in which it is achieved; see Eq. (171); (ii) 
in many cases, and, in particular, in a GW detection scenario with its low signal frequencies and 
heavy test masses separated by the kilometers-scale distances, ordinary solid-state springs cannot 
be used due to unacceptably high levels of mechanical loss and the associated thermal noise. 

At the same time, in detuned Fabry-Perot cavities, as well as in the detuned configurations 
of the Fabry-Perot-Michelson interferometer, the radiation pressure force depends on the mirror 



displacement (see Eqs. (312)), which is equivalent to the additional rigidity, called the optical 
spring, inserted between the cavity mirrors. It docs not introduce any additional thermal noise, 
except for the radiation pressure noise -Fb.a., and, therefore, is free from the latter of the above 
mentioned shortcomings. Moreover, as we shall show below, spectral dependence of the optical 
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rigidity K{fl) alleviates, to some extent, the former shortcoming of the 'ordinary' rigidity and 
provides some limited sensitivity gain in a relatively broad band. 

The electromagnetic rigidity was first discovered experimentally in radio- frequency systems |31) . 
Then its existence was predicted for the optical Fabry-Perot cavities [30j . Much later it was shown 
that the excellent noise properties of the optical rigidity allows its use in quantum experiments 
with macroscopic mechanical objects [55J [^Sl [2U] . The frequency dependence of the optical rigidity 
was explored in papers [371 1871 I38j . It was shown that depending on the interferometer tuning, 
either two resonances can exist in the system, mechanical and optical ones, or a single broader 
second-order resonance will exist. 

In the last decade, the optical rigidity has been observed experimentally both in the table-top 
setup [52] and in the larger prototype interferometer |113j . 

6.3.2 The optical noise redefinition 

In detuned interferometer configurations, where the optical rigidity arises, the phase shifts between 
the input and output fields, as well as between the input fields and the field, circulating inside the 
interferometer, depend in sophisticated way on the frequency Q. Therefore, in order to draw full 
advantage from the squeezing, the squeezing angle of the input field should follow this frequency 
dependence, which is problematic from the implementation point of view. Due to this reason, 
considering the optical-rigidity-based regimes, we limit ourselves to the vacuum-input case only, 
setting Ssqz[r,6l] = I in Eq. ([375]). 



In this case, it is convenient to redefine the input noise operators as follows: 






lia- 

rig- 

72". 



V729' 

V72"' 



Ed", 



(468) 



where 



i-1 



and 



ri 



7i. 

7 



Vd 



(469) 



is the unified quantum efficiency, which accounts for optical losses both in the interferometer and 
in the homodyne detector. 

Note that if the operators d, g, and n describe mutually-uncorrelated vacuum noises, then the 



same is valid for the new dnewj Anew: and rincw 



noises ( 468 ) and renaming them, for brevity, 



— ?► a, 



Expressing Eqs. (371 and 373) in terms of new 
rincw -^ n , (470) 



we obtain: 






v{n) 



2MJ7 HT(0Lo)D(r!) 

T 

L(r2)d(m, 





HT(0Lo)[R(^)a(^) + e"(^)] 



y/2hMJ-f 



(471) 
(472) 



is the lossless cavity reflection factor; see Eq. (550) 



where 

Thus, we have effectively reduced our lossy interferometer to the equivalent lossless one, but 
with less effective homodyne detector, described by the unified quantum efficiency rj < rjd- Now we 
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can write down explicit expressions for the interferometer quantum noises ( 376 ) , ( 377 ) and ( 378 1 



which can be calculated using Eqs. ( 552 1 

SxA^) 

where 



4M J7?7 r2 sin^ ip + Vfi sin^ 0lo ' 

hMj-f{r^ + n^) 

h r cos if — iQ COS 0LO 



2 r sin f — iQ sin (pi^o 



5 

arctan — 

7 



(473) 
(474) 



6.3.3 Bad cavities approximation 



We start our treatment of the optical rigidity with the "bad cavity" approximation, discussed in 
Section |6.1.2| for the resonance-tuned interferometer case. This approximation, in addition to its 
importance for the smaller-scale prototype interferometers, provides a bridge between our idealized 



harmonic oscillator consideration of Section 4.3.2 and the frequency-dependent rigidity case specific 
to the large-scale GW detectors, which will be considered below, in Section [6. 3. 4| 

In the "bad cavity" approximation F ^ fl, the Eqs. (473) for the interferometer quantum 



noises, as well as the expression (374 1 for the optical rigidity can be significantly simplified: 

hr^ „ hMJj 
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XT 



cot Lp , 



(475) 



and 



K = 



MJS 



Substituting these equations into the equation for the sum quantum noise (cf. Eq. (144l) 

S'^in) = \K-Mn^\^Sxx + '2Re{[K-Mn^]SxT}+Sjrjr 
= \K,s-Mn^\^Sxx + Sf^, 

where 
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stand for the effective back-action noise and effective optical rigidity, respectively, and dividing by 



'SQL fill defined by Eq. (161), we obtain the SQL beating factor (418) 
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(479) 



{5 — ^rj sin2p) 



(480) 



gpff ^ 

M n' 

is the effective resonance frequency (which takes into account both real and virtual parts of the 
effective rigidity Kcff). Following the reasoning of Section 4.3.2 it is easy to see that this spectral 
density allows for narrow-band sensitivity gain equal to 
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In the ideal lossless case (77 = 1), 



An 



= e{^m±An/2), 



(482) 



(483) 



in accord with Eq.( 171 ). However, if 77 < 1, then the bandwidth, for a given ^ lessens gradually as 



the homodyne angle ip goes down. Therefore, the optimal case of the broadest bandwidth, for a 
given ^, corresponds to 1^9 = tt/2, and, therefore, to S^f = [see Eqs. (475)], that is, to the pure 



'real' rigidity case with non-correlated radiation-pressure and shot noises. This result naturally 
follows from the above conclusion concerning the amenability of the quantum noise sources cross- 
correlation to the influence of optical loss. 

Therefore, setting (p = 7r/2 in Eq. (479) and taking into account that 
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where n is the normalized optical power defined in Eq. (417), we obtain that 
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(485) 



Consider now the local minimization of this function at some given frequency Oq, similar to 



one discussed in Section 6.1.2 Now, the optimization parameter is /3, that is, the detuning 5 of 
the interferometer. It is easy to show that the optimal /3 is given by the following equation: 
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This fifth-order equation for tan/3 cannot be solved in radicals. However, in the most interesting 
case of Oo ^ Oq, the following asymptotic solution can easily be obtained: 



thus yielding 
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(487) 



The function ( 485 ) , with optimal values of /3 defined by the condition ( 486 ) , is plotted in 



Figure |43] for several values of the normalized detuning. We assumed in these plots that the unified 
quantum efficiency is equal to 77 = 0.95. In the ideal lossless case 77 = 1, the corresponding curves 
do not differ noticeably from the plotted ones. It means that in the real rigidity case, contrary to 
the virtual one, the sensitivity is not affected significantly by optical loss. This conclusion can also 
be derived directly from Eqs. (485) and (488 1. It stems from the fact that quantum noise sources 
cross-correlation, amenable to the optical loss, has not been used here. Instead, the sensitivity 
gain is obtained by means of signal amplification using the resonance character of the effective 
harmonic oscillator response, provided by the optical rigidity. 
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Figure 43: Plots of the SQL beating factor ( |485 ) of the detuned interferometer, for different values 
of the normalized detuning: < /? = arctan((5/7) < 7r/2, and for unified quantum efficiency 
rj = 0.95. Thick solid line: the common envelope of these plots. Dashed lines: the common 
envelopes (420) of the SQL beating factors for the virtual rigidity case, without squeezing, r — 0, 
and with 10 dB squeezing, e^'' = 10 (for comparison). 



The common envelope of these plots (that is, the optimal SQL-beating factor), defined implicitly 



by Eqs. (485) and (486), is also shown in Figure 43 Note that at low frequencies, il <C rjq, it can 
be approximated as follows: 



^2 



(actually, this approximation works well starting from fl < 0.3Qq). It follows from this equation 
that in order to obtain a sensitivity significantly better than the SQL level, the interferometer 
should be detuned far from the resonance, (5 S> 7. 

For comparison, we reproduce here the common envelopes of the plots of ^^(il) for the virtual 
rigidity case with t] = 0.95; see Figure 36 (the dashed lines). It follows from Eqs. (489) and (420) 



that in absence of the optical loss, the sensitivity of the real rigidity case is inferior to that of 
the virtual rigidity one. However, even a very modest optical loss value changes the situation 
drastically. The noise cancellation (virtual rigidity) method proves to be advantageous only for 
rather moderate values of the SQL beating factor of ^ > 0.5 in the absence of squeezing and £, > 0.3 
with 10 dB squeezing. The conclusion is forced upon you that in order to dive really deep under 
the SQL, the use of real rather than virtual rigidity is inevitable. 

Noteworthy, however, is the fact that optical rigidity has an inherent feature that can complicate 
its experimental implementation. It is dynamically unstable. Really, the expansion of the optical 



rigidity (374) into a Taylor series in Q gives 






n 



pi 



(490) 



The second term, proportional to ifi, describes an optical friction, and the positive sign of this 
term (if J > 0) means that this friction is negative. 
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The corresponding characteristic instabihty time is equal to 



^- = 27^ ■ (4^1) 



In principle, this instability can be damped by some feedback control system as analyzed in |37 [ l48 j . 
However, it can be done without significantly affecting the system dynamics, only if the instability 
is slow in the timescale of the mechanical oscillations frequency: 

Taking into account that in real-life experiments the normalized optical power J is limited for 
technological reasons, the only way to get a more stable configuration is to decrease ^, that is, 
to improve the sensitivity by means of increasing the detuning. Another way to vanquish the 
instability is to create a stable optical spring by employing the second pumping carrier light with 
opposite detuning as proposed in j52l I131J . The parameters of the second carrier should be chosen 



so that the total optical rigidity must have both positive real and imaginary parts in Eq. (490): 



MJi(5i AfJ2(52\ .^ /2Af Ji7i(5i 2MJ2J2S2 



K^U^) = K,{n) + K^m « ^-^ + -^ j + ^n ^ ^y + rf ) +■■■ (493) 

that can always be achieved by a proper choice of the parameters J1.2, 71.2 and 61,2 {S162 < 0). 

6.3.4 General case 

Frequency-dependent rigidity. In the large-scale laser GW detectors with kilometer-scale arm 
cavities, the interferometer bandwidth can easily be made comparable or smaller than the GW sig- 



nal frequency Q. In this case, frequency dependences of the quantum noise spectral densities (3761 



( 377 ) and ( 378 1 and of the optical rigidity ( 374 1 influence the shape of the sum quantum noise 



and, therefore, the detector sensitivity. 

Most quantum noise spectral density is affected by the effective mechanical dynamics of the 



probe bodies, established by the frequency-dependent optical rigidity (374). Consider the charac- 
teristic equation for this system: 

- n^[{-f - inf + s^] + js = . (494) 

In the asymptotic case of 7 = 0, the roots of this equation are equal to 





(495) 

(hereafter we omit the roots with negative- valued real parts). The corresponding maxima of the 
effective mechanical susceptibility: 

are, respectively, called the mechanical resonance (f2m) and the optical resonance (flo) of the 
interferometer |37| . In order to clarify their origin, consider an asymptotic case of the weak 
optomechanical coupling, J ^ S^. In this case, 



"™"^ = V™' ''°"'^- ^'^'^ 
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Figure 44: Roots of the characteristic equation (494) as functions of the optical power, for ^/S 



0.03. Solid lines: numerical solution. Dashed hues: approximate solution, see Eqs. (498) 



It is easy to see that flm originates from the ordinary resonance of the mechanical oscillator 
consisting of the test mass M and the optical spring K [compare with Eq. (476)]. At the same 
time, fto, in this approximation, does not depend on the optomechanical coupling and, therefore, 
has a pure optical origin - namely, sloshing of the optical energy between the carrier power and 
the differential optical mode of the interferometer, detuned from the carrier frequency by S. 

In the realistic general case of 7 7^ 0, the characteristic equation roots are complex. For small 
values of 7, keeping only linear in 7 terms, they can be approximated as follows: 
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(498) 



Note that the signs of the imaginary parts correspond to a positive dumping for the optical reso- 
nance, and to a negative one (that is, to instability) for the mechanical resonance (compare with 



Eq. (4901) 



In Figure 44 



the numerically-calculated roots of Eq. (494) are plotted as a function of the 
normalized optical power J/ 6^, together with the analytical approximate solution (498), for the 
particular case of j/S = 0.03. These plots demonstrate the peculiar feature of the parametric 
optomechanical interaction, namely, the decrease of the separation between the eigenfrequencies of 
the system as the optomechanical coupling strength goes up. This behavior is opposite to that of 
the ordinary coupled linear oscillators, where the separation between the eigenfrequencies increases 
as the coupling strength grows (the well-known avoided crossing feature). 
As a result, if the optomechanical coupling reaches the critical value: 



•^" 4 ' 



(499) 
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then, in the asymptotic case of 7 = 0, the eigenfrequencies become equal to each other: 

S 



n(^^ = o(") = Oo 



n/2 



(500) 



If 7 > 0, then some separation retains, but it gets smaller than 27, which means that the cor- 
responding resonance curves effectively merge, forming a single, broader resonance. This second- 
order pole regime, described for the first time in [87', promises some significant advantages for 
high-precision mechanical measurements, and we shall consider it in more detail below. 

If J < Jcrit , then two resonances yield two more or less separated minima of the sum quantum 
noise spectral density, whose location on the frequency axis mostly depends on the detuning S, and 
their depth (inversely proportional to their width) hinges on the bandwidth 7. The choice of the 
preferable configuration depends on the criterion of the optimization, and also on the level of the 
technical (non-quantum) noise in the interferometer. 
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Figure 45: Examples of the sum noise power (double-sided) spectral densities of the detuned 
interferometer. 'Broadband': double optimization of the Advanced LIGO interferometer for NS-NS 
inspiraling and burst sources in presence of the classical noises [97] (J = JaLiGO = (2vr x 100)'^ s~^, 
r = 3100 s~^, (3 = 0.80, 0LO = 7i'/2 — 0.44). 'High-frequency': low-power configuration suitable 
for detection of the GW signals from the millisecond pulsars, similar to one planned for GEO 
HP dSH] [J = O.lJaLiGO, r = 27r X 1000 s-\ /? = 7r/2 - 0.01, 0lo = 0]. 'Second-order pole': 
the regime close to the second-order pole one, which provides a maximum of the SNR for the GW 
burst sources given that technical noise is smaller than the SQL [S'tcch = O.IS'sql, J = •/aLiGO, 
r = 1050 s~\ (3 = n/2 - 0.040, cfiLo = 0.91]. In aU cases, rjd = 0.95 and the losses part of the 
bandwidth 72 = 1.875 s~^ (which corresponds to the losses ^arm = 10~^ per bounce in the 4 km 
long arms). 

Two opposite examples are drawn in Figure |45| The first one features the sensitivity of a 
broadband configuration, which provides the best SNR for the GW radiation from the inspiraling 
neutron-star-neutron-star binary and, at the same time, for broadband radiation from the GW 
burst sources, for the parameters planned for the Advanced LIGO interferometer (in particular. 
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the circulating optical power Iam\ — 840 kW, L = 4 km, and M = 40 kg, which translates to 
J = JaLiGO = (27r X 100)'^ s^^ , and the planned technical noise). The optimization performed 



in [57] gave the quantum noise spectral density, labeled as 'Broadband' in Figure 45 It is easy to 
notice two (yet not discernible) minima on this plot, which correspond to the mechanical and the 
optical resonances. 

Another example is the configuration suitable for detection of the narrow-band GW radiation 
from millisecond pulsars. Apparently, one of two resonances should coincide with the signal fre- 
quency in this case. It is well to bear in mind that in order to create an optical spring with 
mechanical resonance in a kHz region in contemporary and planned GW detectors, an enormous 
amount of optical power might be required. This is why the optical resonance, whose frequency 
depends mostly on the detuning 5, should be used for this purpose. This is, actually, the idea 
behind the GEO HF project |168j . The example of this regime is represented by the curve labeled 
as 'High-frequency' in Figure |45| Here, despite one order of magnitude less optical power used 
(J = O.UaLiGo); several times better sensitivity at frequency 1 kHz, than in the 'Broadband' 
regime, can be obtained. Note that the mechanical resonance in this case corresponds to 10 Hz 
only and therefore is virtually useless. 

The second-order pole regime. In order to clarify the main properties of the second-order 
pole regime, start with the asymptotic case of 7 — >■ 0. In this case, the optical rigidity and the 



mechanical susceptibility (|496|) read 



K{^) - #^> (501) 



1 1 j;2 _ q2 

-^"^ ^ ' K{n) - MQ^ M js - <52J72 + n^ ■ ^ ' 



If condition (4991 is satisfied, then in the close vicinity of the frequency flo (see Eq. (500[)): 



\n-no\<^^o, (503) 

the susceptibility can be approximated as follows: 

Cff,dbl/^N % /QQ^N 

X.x ^"^ ~ Af (f^2 _ ^2)2 ■ ^^"*^ 

Note that this susceptibility is proportional to the square of the susceptibility of the ordinary 
oscillator, 

and has the second-order pole at the frequency fta (thus the name of this regime). 

This 'double-resonance' feature creates a stronger response to the external forces with spectra 
concentrated near the frequency J7o, than in the ordinary harmonic oscillator case. Consider, for 
example, the resonance force Fq sin QqI. The response of the ordinary harmonic oscillator with 
eigenfrequency Qq on this force increases linearly with time: 

FqI 
a^osc {t) = WTFT^ ^i^ ^0^ , (506) 

ZIVl i Iq 

while that of the second-order pole object grows quadratically: 

xMt) = -^ (t'cosnot - ^^"j . (507) 
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It follows from Eq. (151) that due to this strong response, the second-order pole test object has a 



reduced value of the SQL around Hq by contrast to the harmonic oscillator. 

Consider the quantum noise of the system, consisting of this test object and the SQM (that 
is, the Heisenberg's-uncertainty-relation-limited quantum meter with frequency-independent and 



non-correlated measurement and back-action noises; see Section 4.1.11, which monitors its posr 



tion. Below we show that the real-life long-arm interferometer, under some assumptions, can be 
approximated by this model. 

The sum quantum noise power (double-sided) spectral density of this system is equal to 



If the frequency ft is close to Hq: 
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(508) 



(509) 



(510) 



(511) 



(compare with Eq. (165)), where the frequency J7g is defined by Eq. (155) 



The same minimax optimization as performed in Section |4.3.2| for the harmonic oscillator case 
gives that the optimal value of i7g is equal to 



Vtq = i\n, 



(512) 



and in this case, 

ein^+v) ^^^^^^^ <e{^o±^n/2) = (^\ . (513) 

|iy|<AO/2 \ ilo J 

Comparison of Eqs. ( 513 ) and ( |171| ) shows that for a given SQL-beating bandwidth AJ7, the second- 
order pole system can provide a much stronger sensitivity gain (i.e., much smaller value of ^^), 
than the harmonic oscillator, or, alternatively, much broader bandwidth Afi for a given value of 
^^. It is noteworthy that the factor (5131 can be made smaller than the normalized oscillator SQL 



2|t/|/ilo (see Eq. (171)), which means beating not only the free mass SQL, but also the harmonic 
oscillator one. 

where the factors ^^ for the 



This consideration is illustrated by the left panel of Figure 
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harmonic oscillator ( 170 1 and of the second-order pole system ( 511 ) are plotted for the same value 
of the normalized back-action noise spectral density {flq/ilo)'^ = 0.01, as well as the normalized 
oscillator SQL (Iml). 



Now return to the quantum noise of a real interferometer. With account of the noises redefini- 



tion (468), Eq. (385) for the sum quantum noise power (double-sided) spectral density takes the 



following form: 



s''{n) = 



ML'^Jj fi2 sin2 ^^^ ^ p2 sin2 ip 
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Figure 46: Left panel: the SQL beating factors ^^ for fig/fio = 0.1. Thick sohd: the second-order 
pole system (5111; dots: the two-pole system with optimal separation between the poles (5281, 



(529); dashes: the harmonic oscillator (170); thin solid - SQL of the harmonic oscillator (1711 



Right: the normalized SNR (526). Solid line: analytical optimization, Eq. (5301; pluses: numerical 
optimization of the spectral density ( 514[ ) in the lossless case (77 = 1); diamonds: the same for 
the interferometer with J = (27r x lOOp s~^, r]d — 0.95 and the losses part of the bandwidth 
72 = 1.875 s~^ (which corresponds to the losses ^arm = 10^^ per bounce in the 4 km long arms). 



Suppose that the interferometer parameters satisfy approximately the second-order pole conditions. 
Namely, introduce a new parameter A defined by the following equation: 



J{5 - 7sin2(/)Lo) 



nt 



nW . 



where the frequency fio is defined by Eq. (500), and assume that 

Keeping only the first non- vanishing terms in v^ , A^, and 7 in Eq. (514), we obtain that 



S'-'ino 



where 




A 



2\2 



Av' 



A^ 



f^o7 . 



V2 



sin20LO 



2o2 



-AY^i 



fi2 = v^^f]o(l + COs2(/)Lo). 



(515) 



(516) 




(518) 



It follows from Eq. (517) that the parameter A is equal to the separation between the two poles of 
the susceptibility x'^^x^^^- 



It is evident that the spectral density (517) represents a direct generalization of Eq. (510) in 



two aspects. First, it factors in optical losses in the interferometer. Second, it includes the case 
of A ^ 0. We show below that a small yet non-zero value of A allows one to further increase the 
sensitivity. 
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Optimization of the signal-to-noise ratio. The peculiar feature of the second-order pole 
regime is that, while being, in essence, narrow-band, it can provide an arbitrarily-high SNR for 
the broadband signals, limited only by the level of the additional noise of non-quantum (technical) 
origin. At the same time, in the ordinary harmonic oscillator case, the SNR is fundamentally 
limited. 

In both the harmonic oscillator and the second-order pole test object cases, the quantum noise 
spectral density has a deep and narrow minimum, which makes the major part of the SNR integral. 
If the bandwidth of the signal force exceeds the width of this minimum (which is typically the case 
in GW experiments, save to the narrow-band signals from pulsars), then the SNR integral (453) 
can be approximated as follows: 



P 



%{n^)[' 



dv 



_^Sh{Sl^ + v) 



(519) 



It is convenient to normalize both the signal force and the noise spectral density by the corre- 
sponding SQL values, which gives: 



P 



^-^^IM^o)|V, 



where 



1 



dv 



is the dimensionless integral sensitivity measure, which we shall use here, and 

An 



For a harmonic oscillator, using Eq. (170), we obtain 



a' = l. 



(520) 
(521) 

(522) 
(523) 



This result is natural, since the depth of the sum quantum-noise spectral-density minimum (which 



makes the dominating part of the integral (521)) in the harmonic oscillator case is inversely pro- 



portional to its width AJl, see Ec^. (171). As a result, the integral does not depend on how small 
the minimal value of f ^ is. 

The situation is different for the second-order pole-test object. Here, the minimal value of ^^ is 
proportional to (Afi)"^ (see Eq. (513l) and, therefore, it is possible to expect that the SNR will 
be proportional to 



a ex 



Indeed, after substitution of Eq. (511) into (521), we obtain: 



V2^q 



1 



(525) 



Therefore, decreasing the width of the dip in the sum quantum noise spectral density and increasing 
its depth, it is possible, in principle, to obtain an arbitrarily high value of the SNR. 

Of course, it is possible only if there are no other noise sources in the interferometer except 
for the quantum noise. Consider, though, a more realistic situation. Let there be an additional 
(technical) noise in the system with the spectral density S'tcch(^)- Suppose also that this spectral 
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density does not vary much within our frequency band of interest Afl. Then the factor a^ can be 
approximated as follows: 

^^ (526) 



a^ = 



1 



where 
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<5'toch(f^o) 
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(527) 



'5'sQLf.m. (^o) 

Concerning quantum noise, we consider the regime close (but not necessarily exactly equal) to 
that yielding the second-order pole, that is, we suppose < A ^ r^g. In order to simplify our 
calculations, we neglect the contribution from optical loss into the sum spectral density (we show 



below that it does not affect the final sensitivity much). Thus, as follows from Eq. (5171, one gets 
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In the Appendix A. 3 we calculate integral Eq. (526) and optimize it over A and f2g. The opti 



mization gives the best sensitivity, for a given value of ^^cch; i^ provided by 



A = n„ 



noC 



tech ■ 
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In this case, 



2^2^ 



(530) 



tech 



The pure second-order pole regime (A = 0), with the same optimal value of J7g, provides slightly 
worse sensitivity: 

'2- ^ . (531) 



VeTse 



tech 



The optimized function (528) is shown in Figure|46|(left) for the particular case oi A = flq = O.lflo- 
In Figure 46 (right), the optimal SNR (530) is plotted as a function of the normalized technical 
noise ^'^ 



tech" 



In order to verify our narrow-band model, we optimized numerically the general normalized 
SNR for the broadband burst-type signals: 
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where S^ is the sum quantum noise of the interferometer defined by Eq. ( 514 ) . The only assumption 



we have made here is that the technical noise power (double-sided) spectral density 



„ ^ 2h 2 
•^tech ML^ni ^^''"^ 



(533) 



does not depend on frequency, which is reasonable, since only the narrow frequency region around 



flo contributes noticeably to the integral (532). The result is shown in Figure 46 (right) for two 

particular cases: the ideal (no loss) case with 77 = 1, and the realistic case of the interferometer 

with J = (27r X 100)'^ s~^, rjd = 0.95 and 72 = 1.875 s~^ (which corresponds to the loss factor of 
A 



10 per bounce in the 4 km long arms; see Eq. (322)). The typical optimized quantum 



noise spectral density (for the particular case of ^.^ 



It is easy to see that the approximations (528) work very well, even if ^1 



tech 



0.1) is plotted in Figure 45 



s,tech 



1 and, therefore, 
the assumptions (|516[) cease to be valid. One can conclude, looking at these plots, that optical 
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losses do not significantly affect the sensitivity of the interferometer, working in the second-order 
pole regime. The reason behind it is apparent. In the optical rigidity based systems, the origin of 
the sensitivity gain is simply the resonance increase of the probe object dynamical response to the 
signal force, which is, evidently, immune to the optical loss. 

The only noticeable discrepancy between the analytical model and the numerical calculations 
for the lossless case, on the one hand, and the numerical calculations for the lossy case, on the 



other hand, appears only for very small values of ^^^^j^ ~ 0.01. It follows from Eqs. (518) and (529) 



that this case corresponds to the proportionally reduced bandwidth of the interferometer, 

7 ^ rio^Lh - 10 S-' (534) 

(for a typical value of Hq ~ lO'^ s~^). Therefore, the loss-induced part of the total bandwidth 72, 



which has no noticeable effect on the unified quantum efficiency 77 [see Eq. (469)] for the 'normal' 
broadband values of 7 ^ 10'^ s~^, degrades it in this narrow-band case. However, it has to be 
emphasized that the degradation of a^, for the reasonable values of Ctech' i^ o^^y about a few 
percent, and even for the quite unrealistic case of ^^^^.j^ = 0.01, does not exceed ~ 25%. 
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7 Conclusion and Future Directions 

In this review, our primary goal was to tell in a clear and understandable way what is meant by 
quantum measurement in GW detectors. It was conceived as a comprehensive introduction to the 
quantum noise calculation techniques that are employed currently for the development of advanced 
interferometric detectors. The target audience are the young researchers, students and postdocs, 
who have just started their way in this field and need a guide that provides a step-by-step tutorial 
into the techniques and covers all the current achievements in the field. At the same time, we tried 
to make this manuscript interesting to all our colleagues from the GW community and, perhaps, 
from other branches of physics, who might be interested in getting themselves familiar with this 
area, not necessarily close to their own research field. 

However, the reality is crude and such a lofty ambition is always a pot of gold at the end of the 
rainbow. Thus, we could not claim this review to be a complete and comprehensive description 
of the field of quantum measurement. We present here a pretty detailed analysis of the quan- 
tum noise features in the first and second generation of GW interferometers, contemplating the 
techniques considered robust and established. However, many hot topics, related to the planned 
third generation of GW interferometers [49l 11271 182] remained uncovered. Here are only some 
of them: (i) xylophone configurations f83l, (ii) multiple-carrier detectors [1321 1131] . (iii) negative 
optical inertia [^Bj, (iv) intracavity detection schemes 23, 22, 88, 90, S7j, etc. It is our determined 
intention to enjoy the great advantages of the format of living reviews and include those topics in 
future revisions of this review. 

We would like to conclude our review by pointing out how the new swiftly-developing areas of 
modern science and technology, not directly related to GW astronomy and detector science, turn 
out to be deeply rooted in the quantum measurement theory developed by the GW community. It 
is amazing how sinuous the ways of scientific progress are. The history of how GW detection and 
quantum-measurement theory developed and interwove might serve as an example thereof. Indeed, 
from the very first steps towards the experimental observation of GWs made by Weber in the early 
1960s [1651 1166] . it was realized that the extreme weakness of interaction between the ripples of 
space-time and matter appeals for unprecedentedly precise measurement. And almost at the same 
time, Braginsky realized that the expected amplitude of the GW-induced oscillations of the bar 
detector signal mode would be on the order of the zero point oscillations of this mode, as predicted 
by quantum mechanics; that is, in order to observe GWs, one has to treat a detector quantum- 
mechanically and as a consequence there will be a quantum back action, setting a limitation on 
the achievable sensitivity, the SQL pT| . 

This serendipity had a powerful impact on the quantum measurement theory development, for 
it set an objective to contrive some ways to overcome this limitation. For decades up to this point, 
it was a purely theoretical discipline having little in common with experimental science, and, fancy 
that, become a vital necessity for GW astronomy. And again, for several decades, GW detection 
has been perhaps the only field where the results of quantum measurement theory were applied, 
mainly in the struggle with quantum noise, considered as a hindrance towards the noble goal of the 
detection of GWs. And only recently, the same optomechanical interaction, begetting quantum 
noise and the SQL in the interferometric GW detectors, has aroused a keen interest among wide 
circles of researchers studying the quantum behavior of macroscopic objects and testing the very 
foundations of quantum mechanics in the macroscopic world [951 I15j . 

All the techniques and concepts developed in the GW community turn out to be highly sought 
by this new field [SD]. Such methods, initially developed for future GW detectors, as back-action 
evasion via properly constructed cross correlation between the measurement and back-action noise 
sources [1621 11591 11581 11601 11611 [5H1 155] , find a use in the optomechanical experiments with micro- 
and nanoscale mechanical oscillators jSH HHl UMl HOll EOl HZQl M\- It turns out that GW 
detectors themselves fit extremely well for testing the fundamental principles of quantum mechanics 
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just for the record low values of the noise, having non-quantum origin, that owes to the ingenuity, 
patience and dedication by an entire generation of experimental physicists [133] . The very fact 
that the mechanical differential mode of the km-scale LIGO detector has been cooled down to 
Teff = 1.4 fxK without any special arrangement, just by modifying the standard feedback kernel 
of the actuators to provide a virtual rigidity, shifting the 10-kg suspended mirrors oscillation 
frequency from ilm/^TT = 0.74 Hz to 150 Hz, where the GW detector is most sensitive [3], tells 
its own tale. Noteworthy also is the experiment on cooling a several-ton AURIGA bar detector 
mechanical oscillation mode to Tcff ~ 0.17 mK |154| . In principle, some dedicated efforts might 
yield even cooling to ground state of these definitely macroscopic oscillators [SHI I108J . 

One might foresee even more striking, really quantum phenomena, to be demonstrated exper- 
imentally by future GW detectors, whose sensitivity will be governed by quantum noise and not 
limited by the SQL. It is possible, e.g., to prepare the mechanical degree of freedom of the inter- 
ferometer in a close-to-pure squeezed quantum state |116j . entangle the differential and common 
motion of the kg-scale mirrors in the EPR-like fashion [611 I117J . or even prepare it in a highly 
non-classical Schrodinger-cat state (1371 I92j . 
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A Appendices 

A.l Input/Output relations derivation for a Fabry— Perot cavity 



Here we consider the derivation of the I/0-relations for a Fabry-Perot cavity given by Eqs. (274) 



in Section 5.2 We start by writing down the input fields in terms of reduced complex amplitudes 



^1,2 for the classical part of the field and annihilation operators ai_2 for quantum corrections, 



respectively, as prescribed by Eqs. (48) and (51 1 
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with c.c. standing for "complex conjugate" and h.c. — for "Hermitian conjugate". It would be 
convenient for us to make the preliminary calculations in terms of complex amplitudes before going 
for two-photon quadrature amplitudes. 

Start with the equations that connect the classical field amplitudes of the fields shown in 



Figure 29 At each mirror, the corresponding fields are related according to Eqs. (27) with mirror 

real • 



transfer matrix 
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and two equations for the running waves inside the cavity: 
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describes the free propagation of light between the mirrors. The solution to these equations is the 
following: 
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The equations set for the quantum fields has the same structure, but with more sophisticated 
boundary conditions, which include mirrors' motion as described in Section 2.2.5[ 



where 
and 
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The solution to these equations reads: 
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A.2 Proof of Eq. (376), (377) and (378) 



In the lossless case [see Eq. (380l], the noises (double-sided) power spectral densities (376), (377) 
and ( 378 ) are equal to 
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Note that the squeezing factor r has gone. 
Taking into account that 
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and 



SxSf — \Sxf\ 



It is easy to show by direct calculation that 
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which proves Eq. ( 379 1 



A. 3 SNR in second-order pole regime 

Substituting Eqs. ( |528| and ( |527[ ) into Eq. ( |526[ ), we obtain that 

^ _2 Qo f°° dy no 1 



TT rj, y_^(V_A2)2 + l + 2s2 72^9 ^(1 + 2s2 + A4)(yr+ 2s2 + A* - A2) 



where 



i/ A f7o . 



rj„ 



r!„ 



1^„ 



If A = (the sub-optimal case) , then 



a2 = -S^(l + 2s2)-3/4. 



/2r2„ 



The maximum of ( 556 1 in A corresponds to 

A 

and is equal to 



2 , 1 + 2s2 



3 ' 



Slo / 1 + 2.s2 



2^/2^J, V 3 



2\-3/4 



In both cases ( 558 ) and ( 560 ) , we should minimize the expression 



(l + 2s2)-3/4 



n„ 



f7,(l + 2g^Lh 



02 

]2 

'-3 



3/4- 



in f2q. It is easy to see that the optimal Qq is equal to 



(556) 



(557) 



(558) 



(559) 



(560) 



(561) 



(562) 



which gives the SNR (5301 for the optimal case (529) and (531) - for the sub-optimal one A = 
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